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: CHAPTER I 
INTRODUCTION 

Diffusion Is a fundamental process and alumina Is both a widely 
occurring and widely used ceramic. Interest In oxygen diffusion In 
a-Al^Og naturally follows. Such processes as sintering, creep, 
electrical conductivity, and oxidation of metals may all Involve the 
mass transport of oxygen. The motivation for this work stems from 
the latter, In particular the oxidation of NICrAl alloys. 

These alloys exhibit good oxidation resistance due to an a-Al,,0 3 
scale which forms on their surface and across which material must be 
transported for further oxidation to take place. It Is believed that 
the growth of the oxide arises primarily due to an Inward diffusion 
of oxygen from the gas across the scale to the metal, with outward 
diffusion of aluminum along short circuit paths playing a secondary 

The oxide scale tends to spall. This has been attributed to 

growth of oxide within the existing oxide occurring when Inward 

oxygen and outward aluminum fluxes coexist. Small additions of 

oxygen-active elements such as zirconium and yttrium reduce or el 1ml- 
( 1 - 4 ) 

nate this problem' . The postulate has been that these materials 
segregate to the short circuit paths and Inhibit outward aluminum 
diffusion. 
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The oxygen which diffuses across the scale sees a complex envi- 
ronment. The scales are polycrystalllne and the grain size varies 
across the scale. One side of the scale Is exposed to an oxidizing 
atmosphere and the other, which Is In contact with the metal. Is In 
a strongly reducing environment. This scale, while primarily A1 2 0 3 , 
will pick up Impurities as It grows. Concentration gradients will 
likely occur. 

In order to get basic knowledge about this problem. It was the 
Intent of this work to study oxygen diffusion In single crystals of 
o-Al^Oj singly doped with Cr, N1, Y, and Zr, In both a 1 atm 0 2 and 
In a C0/C0 ? atmospheres of low P_ . These dopants were chosen be- 

cause they are likely contaminants In the scales. 

Chromium would not be expected to produce defects since It Is 
Isovalent and the difference In Ionic radii between chromium and 
aluminum Is small. Yttrium Is also Isovalent, however It has an 
Ionic radii much larger than that of aluminum and It has been sug- 
gested that this strain can cause yttrium to act as a donor^. 
However, as will be discussed, the current author questions the 
experimental evidence. 

Nickel Is divalent with small misfit and Is expected to behave 
as an acceptor. Since a nickel Ion on an aluminum site Is a negative 
defect and the crystal must retain electroneutrallty, the Incorpora- 
tion of nickel will necessitate an Increase In the concentration of 
positive defects. The positive defect on the anion sublattice Is an 


oxygen vacancy. 
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Zlrconlum Is quadrivalent with large misfit and Is expected to 
behave as a donor. Zirconium on an aluminum site Is a positive 
defect. The concentration of oxygen Interstitials which are nega- 
tive defects would be expected to Increase with Zr-doplng. 

A study of crystals whose behavior Is governed by these dopants, 
then, provides Information about the mechanism of oxygen diffusion 
In, and the defect structure of, a-Al^. 



CHAPTER II 


LITERATURE SURVEY 

In this section the literature on self-diffusion and the defect 
structure In alumina will be reviewed. Properties such as electrical 
conductivity, sintering, creep, and deformation are related to diffu- 
sion through the defect structure. Therefore, It Is appropriate to 
also review the studies of these properties. 

2.1 MATHEMATICS OF SELF-DIFFUSION MEASUREMENTS 

At temperatures above absolute zero, atoms In a solid are In 
constant motion about their lattice sites. When an atom obtains 
enough energy It may jump to a neighboring site. Such a process Is 
termed diffusion. The random motion of atoms In the absence of a 
concentration gradient Is termed self-diffusion. 

The study of self-diffusion Is most often accomplished by Intro- 
ducing a tracer to the system, that Is atoms which are similar to 
those under study yet In some way Identifiable. Most often an Iso- 
tope with a small natural abundance Is used. The tracer Is allowed 
to make contact with a sample of material at a given temperature for 
a given period of time. After which, the tracer concentration 
profile Is determined. A tracer diffusion coefficient may be deter- 
mined knowing the Initial profile, the final profile, and the 
appropriate solution to the diffusion equation. The relationship 
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between the tracer diffusion coefficient and that of the species of 
Interest Is dependent on their similarity. In the case of Isotopes, 
the difference In diffusion coefficients Is related to the Isotope 
mass difference^ by: 


°1 « 2 
oj - 1 ■ f i ‘<s;> 


( 2 . 1 - 1 ) 


Where D 1 Is the Isotope diffusion coefficient, Is the Isotope 
mass, and f^ Is the Isotope correlation coefficient. The value of 
f.j depends on the crystal structure and mechanism of diffusion. 

When the Isotope oxygen-18 Is used to study oxygen diffusion, the 
difference, neglecting correlation, Is roughly 6 percent. 

In a gas exchange experiment, an atmosphere of 0 ? enriched with 
tracer was brought In contact with the sample at a given temperature. 
The Initial profile In the sample Is only the natural abundance, 

0.204 percent, which can be subtracted out. The final profile Is 
then compared with a concentration profile predicted by a solution to 
the diffusion equation. 

Flck's First Law, In one dimension, may be written as: 

J 1 ■ -») • o iT 

2 

where Is the flux density (kg/m s) of the 1th species, 0^ Is 

2 3 

the diffusion coefficient (m /s), and p Is the density (kg/m ). 

With conservation of mass, and constant density. It may be shown^ 7 ) 

that : 

!ii . sl o ii 

at ax 1 ax 


(2.1-3) 
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Wlth the assumption that 0^ Is Independent of position this 
becomes: 



a 2 c, 


' 3 x 2 


( 2 . 1-0 


Solutions to this equation, using the appropriate boundary 
conditions, which give the tracer concentration as a function of 
position, x, and time of anneal, t, may be obtained. 

For diffusion Into a semi-infinite slab from a gaseous source 
of tracer, which Is well-mixed and Infinite, the following three 
cases are considered. 

(I) The surface Is maintained In equilibrium with a gas con- 
taining a constant concentration of tracer. The boundary conditions 
In this case are: 

c(x,t) = c for t<0 0 < x < « 

OP “ 

c(x,t) = c Q for t > 0 x = 0 
c(x,t) = c as x -* ® 

OB 

The solution to equation (2.1-4) for these conditions 1s^: 

c(x,t) - c = (c - c ) erfc(e) (2.1-5) 

co o co 

where e = x/2 Vot. 

(II) The surface does not Instantaneously equlllbrlate with the 
gas containing a constant concentration of tracer. The boundary 
conditions for this case are: 

c(x,t) = c for t < 0 0 < x < ® 

OB — — • 

fr = K[C c(0,t)] for t>0 
dX x=0 9 


-0 
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c(x,t) = c as x -* ® 

00 

where K Is the first order surface exchange coefficient, C g Is 
the concentration of tracer In the gas, and C(0,t) Is the concentra- 
tion of tracer In the solid at the surface. The solution to equation 
(2.1-4) for this case^ Is: 

c(e, 4 >) - c^ ? 

— = [erfc(e) - exp(4> + 264 >)erfc (0 ♦ $)] (2.1-6) 

c g * ® 

t 1/2 

where <t> = K (p) 

It has been shown that for <t> > 6 the rate of surface exchange 

/ o \ 

may be considered Infinite' . 

(Ill) The surface does not Instantaneously equlllbrlate and the 
concentration of tracer In the gas Is allowed to vary with time In a 
manner described by: 

C g (t) = A 2 t + A^ 172 + A q 

The solution to equation (2.1-4) 1$^: 

E 2 A t n/2 - C d 

( -D n+1 7 — • (2.1-7) 

^n 

n=0 

n 

r (1 + h (exp [«*»( 26 + 4 >) ] erfc (0 + 4>) - £ (-2$) 1 erfc (0)} 

r=0 

where r Is the tabulated gamma function, 1 r erfc Is the rth Inte- 
gral of the complementary error function, and 3 Is the Kronecker 

no 

delta (1 for n = o and 0 otherwise). 
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Oefects must be present In the lattice for diffusion to take 
place. Diffusion mechanisms are usually thought of In terms of point 
defects, either vacancies or Interstitials, although defect complexes 
may play a role. 

An atom on a lattice site adjacent to a vacancy may move to the 
vacancy's site. The vacancy may be thought of as moving In the oppo- 
site direction of the mass flow. 

An atom In an Interstitial site may move directly to another 
Interstitial site or It may move to a normally occupied site and 
knock the lattice atom to an Interstitial site. 

The tracer diffusion coefficient can be related to the defect 
diffusion coefficient by: 

D - fc d O d (2.1-8) 

where D Is the tracer coefficient, f Is the correlation factor 
(1 for Interstitials, <1 for Interstltlalcy and vacancy), c d Is 
the fractional concentration of defects, and D d Is the defect 
diffusion coefficient. 

From random walk considerations It can be shown^^ that: 

1 - 2 

°d = 6 a F d (2J - 9) 

_2 

where a Is the mean square of the Jump distance and Is the 


Jump frequency. 

Since Jumping Is an activated process, the defect Jump frequency 
depends on temperature, as described by: 


r 


d 



v exp 


( 2 . 1 - 10 ) 
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where u Is the vibrational frequency and AG^ Is the activation 
energy for defect migration. 

For the case of thermally generated defects, termed Intrinsic, 
the relation for the concentration of defects to temperature Is: 

-AG, 

C(j = expf-^p 1 ) (2.1-11) 

From equations (2.1-8) to (2.1-11): 

, _2 -AG -AG, 

D = i fa v(exp -^JexpJ-pY 1 ) (2.1-12) 

Since AG = AH - TaS, this may be rewritten as: 

_2 (AS + AS ) -(AH + AH I ) 

D = £ fa « exp (-^ )exp ( ^ ) (2.1-13) 

For the case of defects produced by the presence of an allo- 
valent Impurity, termed extrinsic, where there Is no appreciable 
temperature dependence on defect concentration, the expression for 
defect concentration Is: 

c d = bc f (2.1-14) 

where b Is a number dependent on the charge of the Impurity atoms 
and c^ Is the concentration of the Impurity atoms. The value of 
b is Independent of temperature If the defects produced by the 
Impurity are the majority defects contributing to diffusion. However 
when the diffusion Is by a minority species and the Impurity creates 
a majority defect, b will contain a reaction coefficient and have a 
temperature dependence. 

From equations (2.1-8) to (2.1-10) and (2.1-14): 
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1 - 2 AG m 

0 = £ fa vbCf exp (gy ) (2.1-15) 

which may be rewritten as: 

i - 2 AS m AH 

D - £ fa t,bc f exp (-f) exp ('^) (2.1-16) 

Experimental diffusion coefficients can be represented by an 
Arrhenius type equation: 

0 = D o exp(^) (2.1-17) 

2 

where 0 Q Is the pre-exponential term (m /s) and Q Is the activa- 
tion energy (kJ/mol). Comparison of equations (2.1-13), (2.1-16), 
and (2.1-17) shows the different quantities measured when the mate- 
rial behaves Intrinsically and extrlnslcally. 

Intrinsic Extrinsic 

_2 AS + AS _2 AS 

0 Q £ fa »exp ( ^ — -) ^ fa vbc f exp (-^ ) 

0 AH, ♦ AH AH 

t m m 

2.2 SELF-DIFFUSION STUDIES 

There have been four studies conducted, since 1960, on the self- 
diffusion of oxygen In alumina and one study of aluminum self- 
diffusion In alumina. 

The first was that of Olshl and K1ngery^\^. In this work a 
finite amount of a-Al^O^, In the form of either spheres (0.30 In. 
diameter) or crushed Verneull grown rods (200 to 300 mesh) was annealed. 
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1n an atmosphere enriched In oxygen-18 at the temperature range 1450° 

to 1750° C. The powder was then reduced, the oxygen converted to CO,,, 

18 16 

and the ratio of CO^ / CO 2 was determined by mass spectrometry. 

An Arrhenius plot of log^ q( D) versus 1/T showed a break In 
slope at 1650° C. The da1;a at temperatures higher than this were fit 
with: 


. .. ,„-l .-636 kJ/mol. 2. 

0 = 1.9x10 exp ( ) m /s 


(2.2-1) 


The data below this temperature displayed much more scatter and were 
fit by: 


D = 6.3x10 


-12 


exp ( 


-241 kJ/mol 
RT 


p 

mVs 


( 2 . 2 - 2 ) 


This behavior was ascribed to an Intrinsic-extrinsic transition. 
However, when samples In the low temperature region were preannealed 
before the diffusion anneal, the measured diffusion coefficients fell 
on the extrapolation of the high temperature data. 

In Olshl and Klngery's Investigation, It was assumed that the 
exchange between the gas and the solid was sufficiently fast that 
the process was diffusion controlled. Since what was measured was 
the total amount of tracer, either gained by the sample (In most 
cases) or lost by the atmosphere, the presence of a surface exchange 
effect would cause an underestimation of 0. It was not until late 
In the study that the effect of preannealing the specimens prior to 
exchange was Investigated. The fast transport of oxygen to the 
Interior of the material, by dislocations, again undetectable due to 
the measurement of only total amounts of tracer, would lead to over- 
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estlmatlons of the values of 0. This was determined In later work 
to be the cause of 'extrinsic' break. 

The second and third studies of oxygen self-diffusion occurred 

(12 13) 

concurrently. In 1978 one by Reed and Weunsch' ' and the other 
study by Reddy and Cooper^ 15, In 1980. 

In the work of Reed and Weunsch, diffusion of tracer Into a 
semi-infinite slab was Investigated. The source of tracer was an 
0-18 enriched polycrystalllne Al^O^ layer on a single crystal slab. 
This was generated by successive deposition of A1 metal and oxidation 
In an 0-18 enriched atmosphere. 

The diffusion coefficients that were determined could be de- 
scribed by 


D = 6.4x10 


1 OYn / -2.87 kJ/mol . 2. 
®xp ( nj ) m /s 


(2.2-3) 


This Is about a factor of fifteen below that of Olshl and Klngery. 

One-half the single crystal slabs were preannealed In air, prior 
to thick film deposition, at 1650° C for four hours. The other half 
of the slabs were used In the as-polished state. The diffusion 
anneal took place with a pair of samples, one preannealed and one 
not, bound together with the tracer sources In contact. A schematic 
of the configuration Is shown In figure 1. The concentration versus 
depth profile of 0-18 was then determined with an Ion Probe (SIMS) 
and a diffusion coefficient was determined. 

Reed and Weunsch state that those samples which had been pre- 
annealed showed Flcklan profiles and those which had not been pre- 
annealed showed an anomolously high enrichment In the near surface 


-13- 


region. This was attributed to fast transport In the region damaged 
during polishing. 

Examination of Reed's thesis, however, reveals that the thick 
film solution to the diffusion equation, which Reed employed In his 
analysis, does not describe the experimental situation. 

The thick film solution requires that the diffusion coefficient 
not be a function of position. Referring to figure 1 It can be seen 
that this Is not the case. The work of Olshl and Klngery showed that 
diffusion of oxygen In Al^O^ Is enhanced In the presence of grain 
boundaries. Hence the diffusion coefficient In the labelled poly- 
crystal Is different than In the single crystal. Further, Reddy 
showed that the damage produced by polishing Al^ also enchances 
diffusion. Therefore, the diffusion coefficient Is different In the 
as-polished and preannealed single crystals. 

Both of these enhancements will have time dependencies. The 

(141 

work of Monahan and Halloran' suggests that the single crystal 
will consume the polycrystal, and Reed demonstrated that the polish 
Ing damage can be annealed out. 

A clear experimental observation that the thick film solution 
has been violated Is a flux of tracer across x = 0, In this case 
the sample surfaces correspond to x = 0. Such a flux was observed 
In the work of Reed and Weunsch. Equal amounts of tracer were de- 
posited on the samples that were simultaneously annealed. However, 
after annealllng the as-polished samples contained substantially 
greater amounts of tracer than did the preannealed samples, up to 


1 .5 times . 
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Figure 1. - Schematic of the sample configuration employed by Reed. 
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Figures 2 to 5, were made by tracing the figures In Reed's 
thesis and arranging them so that the diffusion couples In this 
experiment are shown. It can be seen that the profiles are not 
symmetric. Other aspects worth noting Include the fact that there 
Is no discontinuity between the samples which Indicates a good 
couple; the degree of skewness seems to Increase with Increasing 
temperature of the experiment; the areas under the Initial and final 
profiles are very nearly the same. 

It Is also puzzling that the as-polished samples which contained 
the most tracer (suggesting higher effective diffusion coefficients) 
also had the steepest concentration gradients, (suggesting lower 
diffusion coefficients). 

In addition, Reed and Weunsch allowed the Initial concentration 
to be an adjustable parameter during the fitting procedure. It Is 
Invalid to assume different values of c q for the as-polished and 
preannealled samples that were bound together during the experiment. 
Figure 6 shows profiles predicted with the reported values of both 
the as-polished and preannealled diffusion coefficients, using the 
thick film solution, compared with Reed's experimental data. The 
lack of fit Is apparent, however. It should be pointed out this pro- 
file was the most skewed. 

Finally, Reddy analyzed Reeds preannealed samples by proton 
activation. He determined diffusion coefficients which were In 
reasonable agreement with those determined by Reed. However In order 
to obtain fits to the profiles, Reddy found It necessary to assume 
that the sample contained more tracer than Reed reported depositing 



Figure 2. - Reed's diffusion couple at 1584° C for 257 hr 





Figure 3. - Reed's diffusion couple at 1676 




Figure 4. - Reed's diffusion couple at 1742° C for 2.0 hr. 
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Figure 5. - Reed's diffusion couple at 1839° C for 



figure 6. - Comparison of the reported 
with Reed's experimental data at 1839 
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on the sample. Note that this Is the opposite of the finding of 
Reed; that tracer left the preannealed samples during the anneals. 

In summary, this determination of diffusion coefficients seems 
to contain major flaws and It would be useful to have a second solid 
state experiment, to see If there Is a true dlscrepency with the gas 
exchange results. 

_5 

In addition, these exchanges were carried out at 10 atm total 

pressure In the furance, both Reed and Reddy pointed out that defect 

chemistry may be a function of P_ and therefore It Is difficult 

U 2 

to compare these results with those carried out at 1 atmosphere. 

In the study of Reddy and Cooper, semi-infinite slabs of A1 0 

c O 

were annealed In an 0-18 enriched atmosphere. The profile of 0-18 
versus depth was determined by the Single Spectrum Proton Activation 
technique (see Chapter 3). After deconvolatlon with a solution to 
Flck's Second Law which took Into account a surface exchange reac- 
tion, this data was fit by: 

0 = 2.66xl0' 2 exp ( ~ 615 R * J — ) m 2 /s (2.2-4) 

This was In good agreement with the high temperature data of Olshl 
and Klngery. 

In addition to Flcklan profiles, when samples of as-received 
material were exchanged, deeply penetrating "tails" of low concentra- 
tion were observed. It was stated that these "tails" could be re- 
duced or eliminated by preannealing the samples In air. 
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As mentioned, Reddy analyzed Reed's five preannealed samples and 
determined diffusion coefficients by Proton Activation which corre 
sponded well with the profiles determined by Ion Probe, Indicating 
that the difference In D did not result from the technique used. 
Reddy also exchanged a slab of Reed's crystal In the apparatus at 
C.W.R.U. and found the resulting profile to be In agreement with the 
rest of his (Reddy's) work. Reddy offered an explanation of the dls 
crepancy under the following assumptions: 

(1) The A1^0 3 In each case was pure 

(11) The P_ In Reed's setup was between 10~ 5 and 10~ 6 , 
u 2 

(111) No clustering, and 

(1v) The mechanism of diffusion Is Interstitial. In pure 
A1 2 0 3 Interstitial oxygen has an oxygen partial pressure dependence 
resulting from the reaction: 

l 0 2 (g) -» oj ♦ 2h (2.2-5) 

then, the equilibrium constant for this reaction Is: 


k 1 " 


[O^fh ]‘ 


, 1/2 


therefore: 


( 2 . 2 - 6 ) 


« l /? 

C 0 1 J • p o (2.2-7) 

1 U 2 

Since the concentration of oxygen Interstitials depends positively on 

P n , any transport by Interstitials would decrease with decreasing P„ . 
°2 ° 2 
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In addition to his Investigations of undoped Al^, Reddy Inves- 
tigated two types of doped samples; one with a divalent Impurity Hg„ 
and one with a quadrivalent Impurity T1. The Mg-doped sample showed 
diffusion very similar to the undoped material. This was attributed 
to the evaporation of Mg from the near surface region, or the forma- 
tion of second phase particles, where the profiles are determined. 

The Tl-doped material yielded values for the diffusion coefficient 
which were lower than those obtained with the undoped material by a 
factor of five. 

Tl-doplng Is expected to enter the material by one of the 
following reactions: 

3T10 2 - 3T1 a1 ♦ v"J ♦ 60* (2.2-8) 

or: 

2T10 2 -» 2T1^i + oj + 30* (2.2-9) 

Either case would cause the concentration of oxygen vacancies to tend 
to decrease and the concentration of oxygen Interstitials to Increase. 
The fact that the diffusion coefficient decreased In the situation 
where the concentration of oxygen Interstitials Increased, was Inter- 
preted as evidence that oxygen diffusion does not take place by 
Interstitial mechanism. This Is Inconsistent with the Reddy's 
explanation of Reed's data. It should be noted that the level of 
doping was high, 800 ppm. The reduction of 0 by a factor of five 
Is surprisingly low. 
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The most recent work on alumina has been done by Olshl, Ando, 
and Matsuhlro^ 17 ^ ; and Olshl, Ando, and Kubota^ 18 ^. This work In- 
volved two procedures . 

The first was a comparison of oxygen diffusion. In a very high 
purity A1,,0 3 * prepared by chemical vapor deposition (CVD), with that 
In the relatively Impure Verneull of Olshl and Klngery. The samples 
were prepared using Olshl and Klngery's method of crushing the mate 
rial and measuring the total amount of tracer absorbed. These re- 
sults correlated very well In both the high and low temperature 
regions with those of Olshl and Klngery. This lack of dependence on 
Impurity content was explained as having the secondary defects pro- 
duced In crushing (e.g., dislocations) as having carried most of the 
tracer. 

The second procedure employed semi-infinite slabs of Al^Og. 

There were preannealed for thirty minutes In natural oxygen at the 
temperature of the diffusion anneal. Two types of samples were used. 
Verneull material, with Impurities: 100T1, 90Fe, 40Ca, 20S1, <10Mn, 

<10Mg, <1 ON 1 , <10Na, <10Cu, and <1 OV ppm; and CVD material with total 
Impurity content reported to be <0.1 ppm. 

It was shown that the behavior had no dependence on composition 
and could be described by: 

0 = 1 . 1 2x1 0 _1 exp ( ~ 649 R * J/mo1 ) m 2 /s (2.2 10) 

This was Interpreted as evidence for Intrinsic behavior. Other 
explanations Include the possibilities that the samples were not as 
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clean as suggested or that the surface was contaminated by furnace 
atmosphere during the exchange. 

Olshl, Ando, and Kubota also Investigated the effect of pre- 
annealing. Samples which were preannealed at 1730° C In air showed 
enhanced diffusion whereas those samples preannealed at 1900° C In 
vacuo do not. This was attributed to the formation of dislocation 
networks or subgrain boundaries, from polygonl zatlon . 

A summary of oxygen self-diffusion data from the four studies 
outlined above Is presented In figure 7. 

While aluminum self-diffusion Is not directly related to oxygen 
self-diffusion, both are dependent on the defect structure of the 
material. Since anion and cation defect populations are related 
through Schottky and Frenkel products, the study of aluminum self- 

( 19 ) 

diffusion by Paladlno and Klngery Is of Interest . 

Diffusion couples were constructed by hot pressing two poly- 
crystals together; one of which was enriched with a radioactive 
tracer, Al 26 , which decays to Mg 26 ♦ B*. The y-rays resulting from 
the annihilation reaction B + + B” — ► 0.511 MeV y-ray were detected. 
The couples were annealed In air and the profile of tracer was sub- 
sequently determined by sectioning. 

The diffusion behavior could be described by: 


2.8xl0' 3 exp ( 


-477 kJ/mol . 2 . 

RT > m /s 


(2.2-11) 


While this was Interpreted as Intrinsic, the Impurity analysis showed 
Hg, Ca, Fe, Mn, and SI all present In concentrations <100 ppm. Also, 



FFICIENT (m /sec) 
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the concentration of Mg Increases with time as the tracer dlslnte 
grates, therefore, It seems doubtful that this Is truly Intrinsic. 

2.3 DEFECT ENERGY CALCULATIONS 

There have been two sets of calculations, that of Dienes 

et al./^ and recently, that of Catlow, James, Mackrodt, and 
( 21 ) 

Stewart' . A summary of the results of the calculations Is pre- 
sented In table 1 . 

Both sets of calculations were made with a shell model-Hott 
Littleton approach. This Included the assumption that the material 
under consideration was completely Ionic. The calculations of Dienes 
et al., and the Hades II calculations of Catlow et al., assume 
dielectric Isotropy, which Is not appropriate for Al^Og because all 
the Ions do not sit In cubic symmetry. The Hades II and Hades III 
calculations of Catlow et al., use a larger number of atoms resulting 
In better convergence properties. Therefore It Is felt that the cal- 
culations of Hades III are the current best estimates. 

The Hades III calculations were made with both an empirical 
Interionic potential derived from fits to available data on cohesive 
energy, lattice parameters, dielectric and elastic constants, and 
with nonemplrlcal Interionic potentials obtained by the method of 
Mackrodt and Stewart. The nonemplrlcal approach allows potentials 
for redox and doping to be calculated In the same manner as those of 


the lattice. 
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TABLE 1. - OEFECT ENERGIES IN 0 -AI 2 O 3 (kJ/defect) 


Defect type Empirical Nonemplrlcal 

Potential Potential 



Hades 

II 

Hades 

III 

Hades 

III 

Dienes 
et al . 

Schottky Quintuplet 

453 

403 

496 

550 

Anion Frenkel Pair 

392 

366 

798 

675 

Cation Frenkel Pair 

644 

504 

679 

956 


The choice of potential greatly affects the energies calculated 
for the defects and hence the type of disorder predicted to dominate. 
The calculation with the empirical potential suggests anion Frenkel 
will predominate, whereas the calculation with the nonemplrlcal 
potential suggests Schottky disorder. 

The empirical potential has the advantage of being derived from 
real crystals. However the Interatomic forces at small distances, 
such as encountered with an Interstitial, are not well determined by 
the experimental values (of cohesive energy, lattice parameters, and 
dielectric and elastic constants) used to determine the potential. 

Calculations on doping were also performed. The results, with 
the nonemplrlcal potential, for Tl-doplng Indicate that It Is 
energetically more favorable for the compensating defect to be a 
cation vacancy, enthalpy of solution = 341 kJ/TI, rather than an 
oxygen Interstitial, 726 kJ/TI. However, when the empirical poten- 
tial was used It was not possible to select one mode of compensation 


over another. 
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In the calculations for Mg-doplng, the results (with both types 
of potentials) show three mechanisms with nearly the same energy 
penalty. These are shown table 2. 

Although Reddy explained the lack of effect of Mg-doplng on the 
basis of evaporation, It Is Interesting to note that the type of d 1 s 
order predicted to have the lowest energy Is Mg substitutional on an 
A1 site compensated for by Mg Interstitials. In this situation It 
Is not necessary to produce defects on the anion sublattice and 
oxygen diffusion could be unaffected. 

Calculations done on the energy reduction gained by clustering 
of defects shows that the defect formation energies will Involve a 
large stabilizing term from clustering and this will probably deter 
mine the type of disorder present. This Is a new result. Virtually 
all previous work on A1 ? 0 3 has been Interpreted on the basis of 
unassociated defects. 

For example. It Is predicted with these calculations, that even 
as little as 1 ppm Tl-doplng, clustering will dominate up to 1500 K. 

The defect clusters Investigated by Catlow et al., show cluster 
Ing tends to favor compensation by vacancies for both divalent and 
quadrivalent Impurities, see table 3, since the amount of energy 
available from clustering Is larger for those clusters Involving 
vacancies, though as mentioned the choice of potential makes a large 
difference In the T1 case. 
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TABLE 2. - CALCULATED ENTHALPIES OF 
SOLUTION OF Mg IN a-Al 2 0 3 


Mode of Mg Incorporation 

AH 

s, 

kJ/mg 

2MgO -> 2Mg‘ 1 ♦ V” ♦ 20* 

303 

3MgO -» 2Mg^ ♦ Mg” ♦ 30 J 

292 

3MgO + Al^ -» 3Mg^ + Alj** + 30* 

346 

3MgO ■» 3Mg'* ♦ 2v“j ♦ 30 * 

558 

MgO -» Mg'* + 0^ 

1348 


TABLE 3. - BINDING ENERGIES 
OF CLUSTER IN a-Al 2 0 3 

Cluster Energy, 

kJ/cluster 

293 
503 
246 


2T ’ai°, 

3T 'ai v ai 

2M 9 A 1 V o 

"■l*Al 


-234 
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2 . 4 DEFECT STRUCTURE OF a-AK Oj 

2.4.1 Defect Diffusion In doped g-Al ^O^ 

The diffusion of point defects has been studies by the movement 
of color boundaries In doped Al^. 

(221 

In 1969, Jones, Coble, and Mogab' ’ Investigated the movement 

of the color boundary In Tl-doped, 5000 ppm, material. The presence 

43 44 

of T1 In Al^Og produced a pink-violet color, whereas oxidized T1 

Is optically Inactive leaving the material colorless. Crystals which 

had been reduced by annealing In hydrogen were heated In air. After 

given Intervals of time the samples were removed and the progression 

of the boundary Into the sample was determined. It was stated that 

the defect with the highest dlffuslvlty would be rate controlling. 

It was found that the movement of the boundary could be described by: 


0 . 4.3X10- 2 exp r 335 n, ? /s 


(2.4-1) 


It was suggested that this dlffuslvlty was that of the point defect 
causing the oxidation. 

( 231 

Hower, Akul'Onok et al., who studied the same system In 
1978 pointed out that the dlffuslvlty of the boundary Is related to 
the defect dlffuslvlty by: 

N 

°boundary 2 C Q °defect (2.4-2) 


where N Is the equilibrium defect concentration and C Is the 
o o 

Initial concentration of Impurity. They suggest the high dopant 
level, 5000 ppm, would make N q /C 0 a small enough number to slgnlfl- 
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cantly change the deduced value of D; however, the activation energy 

term should remain very similar. 

Akul'Onok et al., reduced the T1 by vacuum annealing and found 

good agreement with the value of 0 from Jones' reduced crystals. 

This Indicates that the diffusion of hydrogen was not rate controlling 

Akul'Onok et al., also studied the effect of P n on color bound- 

u 2 

ary motion by annealing crystals In mixtures of He-O.,. A dependence 

of was found which most closely resembled that predicted If 

u 2 

triply Ionized aluminum vacancies were assumed to be the diffusing 
species. Supporting evidence for the rate limiting step Involving a 
negatively charged defect was given In the form of an observed shift 
In color toward a positive electrode when the oxidation anneal was 

II I 

carried out under an electric field. This Indicates the V A1 were 
compensated by holes, but that they are not bound Into a complex. 

The diffusion coefficient for aluminum vacancies, estimated from 
color boundary motion, at 1500° C Is: 

~4xl0~® m 2 /s 

compared with Paladlno and Klngery's self-diffusion coefficient of 
aluminum at 1500° C Is: 
lxl O' 6 m 2 /s 

The difference suggests that the self-diffusion of aluminum takes 
place by a mechanism Involving a considerably less mobile defect. 

If aluminum diffuses by a vacancy mechanism, then: 

Al 


°A1 = [V A1 ]0 V 


(2.4-3) 
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therefore: 


t V Al) * 



using the values of Akul'Onok: 


(2.4-4) 



! ®- 6 - = 2.5x10"® 


4x10 


-9 


(2.4-5) 


This would correspond to a defect concentration of 0.025 ppm which 
Is considerably less than the allovalent Impurity concentration, a 

puzzling result . 

(24) 

Cox' ' observed the motion of a color boundary In Mg-doped 
Al^Og and found It proceeded much faster than the boundary In T1 
material. It was stated that the rate of boundary motion In Mg-doped 
Al^Og at 1000 K was equivalent to the rate In Tl-doped A1 ? 0 3 at 
2000 K. He used this as an argument for an Interstitial cation 
rather than an anion vacancy based on Indications that the mobility 
of the former Is greater than that of the latter. However, Mohapatra 
and Kroger have suggested that what was more likely observed Is the 
diffusion of the hydrogen used In the reduction of the crystal. 

2.4.2 Density Measurements 

(25) 

Rasmussen and Klngery Investigated the effect of the addi- 
tion of dopants on the density and lattice parameters of Al^O^. 

The dopants employed were Ca, Mg, Cr, T1, V, and SI. 


-34- 


The density measurements Indicated that the donor dopants of SI 
and T1 were compensated by aluminum vacancies. There was also evi- 
dence for the solubility limits. In the case of Si-doping, regard- 
less of the concentration of SI In the starting powder from which 
the crystals were grown (250, 500, 750, and 1000 ppm), the resulting 
crystal contained <250 ppm. In the Tl-doped material a second phase 
was observed In the samples with >500 ppm dopant. 

In the crystals doped with the acceptors Ca and Hg, a second 
phase was present In all samples. The effect of the second phase on 
the density measurements was much greater than that arising from the 
production of defects. This made It difficult to draw any conclu- 
sions on the Incorporation mechanism. What It does suggest Is that 
the solubility limit for these Ions Is low, <400 ppm. 

The Cr-doped samples show density change consistent with no 
defect production and the one sample containing vanadium Indicated 
that It was also Incorporated as an Isovalent Ion, V + ^. 

2.4.3 Spectroscopic Investigations 

Discussion In this section will be limited to results on defects 
In undoped Al^O^, defect clustering In acceptor doped Al^Og, and 
A1 ? 0 3 co-doped with acceptors and donors (since In practice all 
samples contain both), and nickel-doped Al 2 0 3 (since this Is our 
acceptor). 

2 . 4 . 3 . 1 Undoped g-Al ^O^ 

Lee and Crawford^) found that when undoped alumina was heated 
to 2000° C In a strongly reducing environment the material became 
colored, with the principle absorption at 6.1 eV. This was attrl- 
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buted to an F-Center, an oxygen vacancy associated with two elec 
trons, on the basis of the observed symmetry. 

The absorption of 6.1 eV. was unaffected by postcoloration 
anneals up to 1300° C and that an anneal of one-half hour at 1400° C 
was insufficient to remove much color. This was contrasted to the 
ease with which coloration produced by neutron irradiation can be 
annealed out . 

The difference in rate of color center removal by annealing may 
be used as an argument for Schottky defects predominating in undoped 
alumina. It is considered that the majority of defects produced by 
neutron irradiation are Frenkel pairs. The neutron knocks an ion 
from its lattice position into an interstitial site. In annealing 
out such Frenkel defects the relevant diffusion distance governing 
the kinetics is that between the interstitial and the vacancy which 
can recombine. The Frenkel reactions for A1.,0 3 are: 


v ai + a 1 T' < 2 - 4 ~ 7 > 

However, in the case of Schottky disorder, the vacancies must travel 
to a sink, et., dislocations or the crystal surface. The Schottky 
reaction for A1 ? 0 3 is: 

2V A i ♦ 3\T ’ -* Null (2.4-8) 

In this case the distance across which the defect must migrate Is 
related to the dislocation density and the sample thickness. This 
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1s presumably a much greater distance than for Frenkels, leading to 
the Increased resistance to the defect removal through annealing. 

It should be noted that there Is some question as to the purity 
of the material since no results of a chemical analysis are given. 

2. 4. 3. 2 Defect Clustering 

The EPR spectra of Al^Og and Al^O^ doped with Iron have been 

121) 

studied by Bauer and Whitmore . The material was doped with 15, 
30, and 230 ppm Fe. In order to account for the observed g-values, 
the symmetry of the centers, and the dependence of the absorption on 
Fe content, the Investigators found It necessary to Invoke clusters 
of defects containing the cation Impurities. 

In a study of Fe-doped, Tl-doped, and Fe-TI co-doped Al^Og, 
Elgenmann and Gunthard^ 28 ^ used Verneull grown material doped with 
Fe at 75 ppm, T1 at 350-600 ppm, Fe-TI at 250 to 1000 ppm Fe, and 50 
to 350 ppm T1 . 

It was reported that the ESR and optical spectra of the Fe-TI 
co-doped material could not be accounted for by the summation of the 
spectra from the singly doped materials. This Is evidence for 
cluster formation. In this case, Fe-TI Ion pairs. 

2. 4. 3. 3 Nickel Doping 

In the spectroscopic studies of Nl-doped Al^Og the valence state 

1291 

of the nickel has been unclear. McClure' ' observed that heating a 
Nl-doped crystal In a vacuum to 1600° C for fifteen minutes (Verneull 
grown with 0.2 mol percent N1 In the starting material) removed the 
coloration associated with N1 +8 but that the coloration associated 
with N1 did not appear. Also observed was that upon heating the 
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+3 

crystal to 1600° C In an oxidizing environment, the N1 color was 
restored. 

In another study of Nl-doped AI^Oj, Marshall, Klkuchl, and 
Relnberg^*^ also observed the color change upon heating the sample 
In a vacuum to 1800° C for two hours; but no change In the EPR spec- 
tra was observed. It was also noted In this study that two Verneull 
grown crystals with quite different amounts of nickel oxide In the 
starting material, 0.4 percent and 4.0 perccent, showing Identical 
EPR spectra with only a slight difference In Intensity. This sug- 
gests that the solubility of nickel In alumina <0.4 percent. 

( 31 1 

Muller and Gunthard' were able to reduce Nl-doped A1 2 0 3 with 
hydrogen at 800° C and found the EPR spectra changed. The amount of 
hydrogen uptake was measured and found to be 1.0+0. 7 H/N1 . It was 
also observed that was present In the crystals stemming from the 
Verneull process. 

In summary, the spectroscopic studies of nickel-doped alumina 
have produced Interesting yet Inconsistent results and It remains to 
be explained why the nickel Is apparently not reduced In vacuum. 

2.4.4 Electrical Conductivity Measurements 

The conductivity In alumina Is very small and Is therefore 

1321 

difficult to measure. In 1965, Peters, Flensteln, and Peltzer' 
found that the conductivity measured by their equipment was the same 
value whether or not the sample was In place. This was attributed 

to conduction through air. 

(331 

Kltazawa and Coble measured the conductivity of undoped 
alumina In the temperature range of 900° to 1700° C, being careful 
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to choose sample configurations so as to block the external circuits 
through the gas phase and along the surface. The resulting conduc- 
tivity versus temperature plot shows two linear regions with differ- 
ent activation energies: 238 kJ below 1300° to 1400° C and 335 kJ 

above 1400° C. 

It was also observed that conductivity was electronic at high 

temperatures (>1600° C) being n-type at low P_ and p-type at high 

u 2 

P Q , and Ionic at lower temperatures. On the basis of the subsequent 
wo?k of Kroger et al. t It seems more likely that the observed con- 
ductivity was mixed Ionic-electric In all temperature ranges. 

Kroger and his co-workers have conducted numerous studies on 
the conductivity of both undoped and doped alumina. It Is appro- 
priate to review these studies In chronological order stating the 
findings of each and offering comments. In none of these studies Is 
a break observed In conductivity versus 1/T, however, these measure- 
ments are confined to the high temperature region (>1300® C) of 
Kltzawa and Coble. 

In 1975, Dutt, Hurrell and Kroger^ 34 ^ studied cobalt-doped 
alumina (Impurity levels In ppm 31 Co, 4.6 Cr, <8 Fe, 52 SI, 1.3 Mg, 
and 5.0 Ca In ppm) using optical absorption, ESR, electrical conduc- 
tivity, and EMF measurements. The presence of cobalt was found to 
Increase the conductivity. When the cobalt Is divalent It Is a 

I 

negative defect Co^ and the concentration of positive defects \T ' , 
Al^’’, and h’ must Increase to compensate. Based on a comparison of 
aluminum and oxygen self-diffusion In alumina It Is postulated that 
the defect responsible for the Ionic conductivity Is Al^”. 
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3Co a1 ♦ Al”' ♦ | 0 2 (g) t 3Co J, ♦ AlJ, ♦ | 0* (2.4-9) 

To reconcile discrepancies between theoretical predictions and 
experimental data It had to be assumed that an unknown donor was 
present, that this donor was fixed valence, and that the donor level 
was above the cobalt level In the bandgap. While this Is possible, 
the presence of the donor was not experimentally verified. The only 
donor which showed up In the spectro-chemlcal analysis was silicon 
and this was thought to be an artifact. 

Dutt and Kroger^^, also In 1975, Interpreted the data of 
Pappls and Klngery on conductivity In alumina and came to the con- 
clusion that aluminum Interstitials were the Ionic charge carries 
and holes were the electronic charge carries. The measurements 
of Pappls and Klngery were made In an Improperly guarded cell. 

A correction factor based on the work of Dutt and Kroger on Co-doped 
alumina. In a different setup, was applied to these measurements, 
which casts some doubt on the reliability of the analysis. 

In 1977, Mohapatra and Kroger^) Investigated Mg-doped Al^Og. 

It was not clear however that the sample was Mg-doped as the spectro- 
chemlcal analysis showed other Impurities of significant concentra- 
tions: 5 Cr, 10 Mg, 5 Fe, 5 SI, 1 Co, 1 Ca, 2 Na ppm. It was also 
stated that between 12 to 90 percent of the total conductivity mea- 
sured was due to surface and gas phase conduction. The authors were 
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unable to conclude whether aluminum Interstitials or oxygen vacancies 

were the defect responsible for Ionic conduction. 

/ 37) 

In 1977, Mohapatra and Kroger also studied Tl-doped alumina. 
Analysis of this material Indicated a well-doped crystal, 430 ppm 
T1 , 6.2 Ca, 1.2 Mg, and 26 SI. The authors suggest that clusters of 
(3T1 ^^V a ^) x may form but they assume, without evidence, that 
the clusters have disassociated at the temperature of measurement, 
1600° C. Further, they conclude that aluminum vacancies are the 
Ionic charge carries and electrons the electronic charge carries, 
although It Is pointed out that a model based on oxygen Interstitials 
Is only slightly less satisfactory. 

/ QQ \ 

In 1978 Mohapatra and Kroger' ' reviewed the data available 

at the time and compared It with the defect energy calculations of 

Dienes et al. The point Is made, by the authors, that to get 

Information on the defect structure In A1 ? Q 3 from conductivity, 

self-diffusion, creep, and sintering Involves testing models against 

P_ dependencies. 
u 2 

X = B ( P ) r [dopant] 5 (2.4-10) 

2 

where X Is the property of Interest and B Is a proportionality 
constant. The exponents r and s are often very similar for 
different models. There Is often a large uncertainty In the experi- 
mentally determined values of these processes In alumina. This Is 
contributed to by the fact that conductivity Is low, self-diffusion 
rates are low, etc. This makes It difficult to distinguish between 
models . 
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The authors then go on to Interpret the available data as Indi- 
cate ng a predominance of Schottky disorder, despite mentioning argu 

ments by others that there Is evidence for Frenkel predominance. 

1 39) 

In 1979, Hohapatra, Tlku, and Kroger': addressed the question 

of finding Al^O^ pure enough to show behavior not dominated by 
Impurities. One sample used, grown on a molybdenum die, displayed 
acceptor dominated behavior despite the fact that the spectrochemlcal 
analysis showed only marginal traces of acceptors 1.3 Hg, 0.7 Cu, 

2.8 ppm Ca, and a sizeable donor concentration, 53 ppm SI, with the 
41 ppm Ho assumed Isovalent with Al. A second sample grown by the 
Czochralskl method Initially showed donor doped behavior which 
changed to acceptor dominated after annealing Is oxygen. The donor 
behavior was attributed to hydrogen Incorporated Into the material 
during growth. In the Interpretation of this data, arguments are, 
again, found for both Schottky and Frenkel defects as the majority 
species . 

(41 ) 

In 1980, Tlku and Kroger' further Investigated the presence 
of hydrogen In alumina and Its action as a donor. Despite the fact 
that data points garnered by extrapolation are plotted In several 
figures and that linear behavior Is sometimes tested with only two 
measurements, there seems to be convincing evidence for hydrogen 
acting as a donor. This Is In accordance with the observation that 
hydrogen may reduce Impurities. Efforts to deduce the dominant type 
of native disorder using this Information were unsuccessful. 

Also In 1982, E 1 -Al at and Kroger^*^ Investigated yttrium as 
a possible donor, due to misfit strain. In alumina. The material 



-42- 




used was obtained from the same source as that used In the present 
work. The crystal we obtained, as will be shown later, had all of 
the yttrium segregated to Its surface and had bulk Impurity levels 
virtually identical to that of our undoped crystal. This raises a 
question as to whether or not the donor doped behavior of the crystal 
In El-Alat's work Is truly due to the presence of yttrium. It Is 
also stated that the crystal used as an undoped standard contains 
"cloudy areas" and behaves as If It were acceptor doped. 

Other points of the paper which tend to diminish Its credibility 
are the fact that the transference numbers are derived from the 
slope of a curve determined by Just three points and the fact that 
the yttrium-doped crystal has a conductivity which Is time dependent 
after a temperature change. It does not seem there Is enough evi- 
dence at present to conclusively Indicate that yttrium does In fact 
act as a donor. 

A conclusion based on the electrical conductivity data In terms 
of unasssoclated point defects seems to be very difficult. In light 
of the new defect calculations It may be that clustering Is of major 
Importance and will have to be taken Into account. Such a review 
and reinterpretation will not be attempted here. 

2.4.5 Sintering Studies 

The sintering behavior of Al^Og has been widely studied. 

In this section only a few points relevent to the present study will 
be brought out. 

The role of both acceptors and donors has been Investigated. 

Rao and Cutler^ 4 ^ studied the effect of Fe and T1 doping. Bagley, 
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/ 441 

Cutler, and Johnson also studied the effect of T1 doping. 

MSI 

Brook' ' further Interpreted their data. The conclusions arrived 
at were based on Frenkel defects with cation diffusion being the rate 
controlling step, with the compensating defects being aluminum Inter- 
stitials for Iron and aluminum vacancies for titanium. 

Mistier and Coble^ 48 ^ using the single crystal and polycrystal 
oxygen diffusion data of Olshl and Klngery and the aluminum diffusion 
of Paladlno and Klngery, under the assumption that aluminum diffusion 
Is not enhanced by grain boundaries, calculated that with a grain 

size below 5 pm, the diffusion of aluminum would be rate controlling. 

(47) 

Delaunay, Huntz, and Lacombe conducted an Investigation 
Into the effect of yttrium on sintering of Al^Og. The particle size 
of the A1 ? 0 3 was 50A, which Is well below the Mistier critical size. 
It was found that the addition of yttrium Increased the observed 
sintering rate, with additions of up to 1 percent. By Inference, 
the addition of yttrium Increases the diffusion of aluminum rather 
than Impedes It, as was suggested by Gollghtly et al. 

2.4.6 Creep Studies 

Creep studies can give Information about the rate controlling 
step. In A1.,0 3 this Is generally believed to be lattice diffusion 
of aluminum, although Lessing and Gordon have shown that with large 
grain size In a heavily acceptor doped material, grain boundary 
diffusion of oxygen may become limiting. 

Studies of acceptor, Fe, and donor, T1, doped material have been 
carried out by Hollenberg and Gordon^ 48 ^; Lessing and Gordon^ 48 ^; 

Hou, Tlku, Wand, and Kroger^ 50 ^; and El-Alat, Hou, Tlku, Wand, and 
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dlffuslon control, with the exception mentioned above. The defect 
In Fe-doped material was felt to be A1 ^ and In Tl-doped material to 
be V-t Gordon and co-workers felt this to be evidence for Frenkel 
disorder while Kroger and co-workers were able to show that this was 
not the only situation which would satisfactorily explain the results. 

Kroger and co-workers^ an< ^ ^ also measured the conductivity 
of the materials used In the creep studies. The results Indicated 
that while there was grain boundary contribution to electronic con- 
duction, no contribution was evident In the Ionic conduction. This 
was Interpreted as evidence for grain boundary diffusion by a neutral 

oxygen species, 0*. 

2.4.7 Deformation Work 

{521 

Pletka, Mitchell, and Heuer' studied the hardening of sap- 
phire with Isovalent and an allovalent Impurity. The hardening 
achieved from the solid solution of Isovalent Impurities Cr 3+ , T1 3+ , 
was due to the development of a strain field resulting from the 

differences In Ionic radii. The much greater effect, an order of 
4+ 

magnitude, of T1 was explained using a model of defect clusters. 
The proposed cluster consisted of a cation vacancy surrounded by 
three quadravalent Impurities. 

Later work by Phillips, Mitchell, and Heuer on the growth of 
dislocation arrays during the precipitation of rutile In Tl-doped 
sapphire (star sapphire) led to the conclusion that the compensating 
defect for the T1 was an oxygen Interstitial rather than an aluminum 


vacancy. 
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( 54) 

Currently, Lagerlof and Heuer' are Investigating the anni- 
hilation of dislocation loops In cylinders of Al^ are de- 

formed, generating a large population of dislocations. The radii of 
small loops are observed by TEM as a function of annealing time. 

Based on measurements of the radii of these loops, and the assumption 
that self-diffusion of oxygen Is rate controlling, an oxygen self- 
diffusion coefficient can be calculated (see fig. 8). The calculated 
diffusion coefficients for undoped Al^Og agree very well with the 
extrapolation of the data of Olshl, Ando, and Kubota. However, In 
Mg-doped Al^ the calculated diffusion coefficients are much higher, 
a factor of 50, than In the undoped case and In Tl-doped Al^ the 
calculated diffusion coefficients are lower, a factor of 5, then In 
the undoped case. This dopant dependence Is In contrast to the Im- 
purity Independence of Olshl, Ando, and Kubota. 

2.4.8 Summary 

The determination of the defect chemistry In terms of un- 
associated point defects has been generally unsuccessful. There Is 
growing evidence that clustering of defects Is Important even with 
low defect concentrations and at high temperature. 
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temperature ( *c > 



RECIPROCAL TEMPERATURE (K*'.|0 4 ) 


Figure 8. - Diffusion coefficients determined from dislocation loop 
annihilation, (figure taken from Reference 54). 
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CHAPTER III 


EXPERIMENTAL PROCEOURE 

3.1 SAMPLE PREPARATION AND CHARACTERIZATION 

3.3.1 Single Crystal Samples 

Single crystals were obtained through NASA Lewis Research 
Center^. These crystals were grown In France by La Pierre 
Synthetlque Blakowskl International. 

The crystals were grown by the Verneull process. It was the 
Intent of the study to have five crystals of different composition, 
an undoped crystal and four singly doped crystals; one with yttrium, 
one with chromium, one with nickel, and one with zirconium. The 
starting material for the doped crystals was reported to have con- 
tained 100 ppm of the desired Impurity. The undoped crystal appeared 
clear and well-formed, that Is, with a smooth surface and an even 
shape. The chromium-doped crystal had a uniform red color and was 
also well formed. The nickel-doped crystal appeared green-yellow, 
however, It had a deeper coloration toward the center of the boule 
than at the edges. It was also a well-formed crystal. Both the 
yttrium-doped and the zirconium-doped crystal showed what appeared 
to be precipitates on the surface of the boule. These precipitates 
were a brown color In the case of the yttrium-doping and a white 
color In the case of the zirconium-doping (see figs. 9 and 10). 

^Courtesy of Or. J. Smlalek, NASA Lewis 
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Whlle the zirconium-doped crystal appeared to be evenly shaped, 
the yttrium-doped crystal had an Irregular shape with three large- 
scale facets having formed along the length of the boule. Both crys- 
tals contained numerous cracks In their Interiors. 

Two types of samples were submitted to Coors Spectrochemlcal 
laboratories for analysis. The results are shown In table 4. 

TABLE 4. - ANALYSIS OF SINGLE CRYSTALS 3 (ppm) 



Sample 

Ca 

Cr b 

Fe 

Mg 

N1 b 

SI 

Na 

yb 

Zr b 

Total 

Undoped 

Cl 

5.6 


16 

23 


27 

4.5 



76.1 


N8/4 

6.0 

— 

14 

30 


34 

5.6 

— 

— 

89.6 


N8/2 

10.0 

— 

14 

26 

-- 

42 

7.0 

— 

— 

99.0 

Y-doped 

C2 

5.8 

r “ — 

14 

27 

_ _ 

37 

5.6 

<10 



89.4 


N10/3 

6.0 

— 

11 

18 

- - 

28 

5.2 

150 

— 

236.2 


N10/1 

4.0 

— 

14 

17 

-- 

52 

6.2 

18 

— 

111.2 

Zr-doped 

C3 

9.4 

» _ 

10 

30 

_ _ 

57 

5.2 



<10 

121.6 


N13/4 

8.5 

— 

16 

17 


28 

7.0 

— 

410 

486.5 


N13/2 

6.0 

r : 

15 

21 


50 

6.3 

— 

86 


Nl-doped 

C4 

4.6 

_ _ _ 

11 

15 

15 

33 

4.2 

— _ - 

— « — 

82.8 


N9/4 

5.4 

— 

9 

16 

30 

32 

5.0 

— 

— 

97.4 


N9/2 

17.0 

— 

18 

17 

22 

37 

5.2 

— 

— 

116.2 

Cr-doped 

C3 

6.0 

118 

24 

20 

_ _ 

33 

5.2 

_ _ _ 



206.2 


N7/2 

17.0 

114 

22 

17 

— 

42 

16.2 

— 

— 

228.2 


N7/3 

8.5 

no 

22 

15 

— 

47 

6.2 

— 



208.7 


3 Coors Spectrochemlcal Laboratory. 
b Dopant . 

The samples with the prefix 'N‘ consisted of a full slice of the 
boule, and those with the prefix ' C ' were slices of the boule from 
which the outer edge, or the "skin," had been ground off with a dia- 
mond lap (45 ym) . In the case of the undoped, Cr-doped, and Nl-doped 
samples, there Is little difference between the two types. The 
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ana lysis of the undoped crystal reveals a level of background Impuri- 
ties around 100 ppm. The Cr-doped crystal had a similar amount of 
background Impurities with the Intended level of dopant, 100 ppm, 
almost precisely being obtained. The Nl-doped crystal, while appear- 
ing colored to the eye, did not have a sufficient amount of nickel 
such that Us presence would dominate the crystal behavior. There 
Is as much magnesium and twice as much silicon. In both the case of 
Y-doped and Zr-doped samples there Is a striking difference between 
the two types. The "surface removed” type had no dopant detected and 
showed background Impurities similar to the undoped crystal. The 
"full slice" types showed the presence of the dopant In 
concentrations which varied dramatically between samples by a factor 
of between five and eight. This Indicates that the dopant Is 
rejected to the surface of the boule during crystal growth. This was 
later confirmed by SIHS analysis, the surface precipitates on each 
boule found rich In the dopant, (figs. 11 and 12). 

A 'single precipitate' was removed from the surface and analyzed 
2 

with an SEM . Prior to analysis the sample was polished with 6, 3, 
and 1 ym diamond paste followed with a Syton etch. The polish was 
done at a slight angle to reveal the boule/preclpltate Interface. 
Figure 13 shows low magnification micrographs of the precipitate. 

It can be seen that the microstructure Is complex but very ordered, 
similar to that of eutectic phase. Also visible are Isolated Islands 
of the eutectic phase and a long "stringer,” which may be responsible 

2 NASA Lewis' AMRAY 1200 
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1H0I3H MV3d 


Figure 11. - Mass scan of precipitates on the Y-doped boule. 



1H0I3H XV3d 
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Figure 13. - Low magnification optical micrographs of precipitate 
from Y-doped boule. a) SEH micrograph of sapphire and precipitate 
showing eutectic type microstructure, b) close-up of sapphire 
precipitate Interface, c) close-up of "stringer" showing facets. 
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for the internal cracking In the boule. It Is Interesting that the 
facets on the “stringer" form an angle which appears to be the same 
as the large scale facets on the boule. Evidently there Is an 
epitaxy problem In the growth of this phase, though no serious 
Investigation Into It was made. 

A region of rather coarse microstructure was Investigated with 
EDAX. An EDAX spectrum Is shown In figure 14. All the peaks could 
be accounted for by Al, Y, Ca, and Au. The Au was deposited to 
dissipate charge In the microscope. The other three were also pre- 
sent In the SIMS. A micrograph and element maps for this area are 
shown In figure 15. 

It Is clear that the phase which appears white In the micro- 
graph Is enriched In yttrium while the dark phase Is enriched In 
aluminum and has a very low yttrium content. The calcium does not 
have a significant correlation to the apparent microstructure and Is 

3 

treated as a trace Impurity . 

The Implication of this Is that during crystal growth the dopant 
Is continuously rejected to the melt until the local concentration 
becomes sufficient to form this eutectic mixture on solidification. 
The fact that spectroscopic samples taken from a crack free Interior 

3 Samples of this boule were given to Bob McCune (at Ford Motor Co.) 
and his analysis found this boule to contain more Ca than the 
boule of Intentionally Ca-doped sapphire purchased from the same 
manufacturer. In addition he found an axial gradient, first grown 
end ~15 ppm last grown >300 ppm. 



000 h 
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SINDOD 


Figure 14. - EOAX spectrum of precipitate from Y-doped boule showing 
the presence of Al, Ca, and Y. The Au was deposited to prevent 
charging. 
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e from Y-doped boule. a)SEM 
ure, b) EDAX map of Yttrium, 
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region (presumably with no “stringers" or Inclusions) did not contain 

detectable amounts of yttrium suggests the solubility of yttrium In 

sapphire Is less than the 10 ppm detectability limit. The eutectic 

phase formed In this system would be expected to be YAG, however It 

has been found that when molten YAG Is quenched a different phase 
(55) 

YAP ' Is formed, which has the perovsklte structure. What Is clear 
Is that samples from this boule contain, at best, levels of yttrium 
which are small compared to the background Impurities. 

While the Zr-doped crystal was not examined In detail, similar 
behavior Is expected. 

The net result Is that we have three crystals with virtually 
Indistinguishable Impurity levels which will be termed undoped no. 1, 
undoped no. 2, and undoped no. 3; a crystal with Its dopant, nickel, 
representing ~20 percent of the total Impurity content which will be 
termed Nl-doped; and a crystal with the Isovalent dopant chromium 
representing ~50 percent of the total Impurity content which will be 
termed Cr-doped. Such a set of samples Is wholly unable to provide 
the Investigator with the Information sought regarding the Influence 
of Cr, N1, Y, and Zr content on the self-diffusion of Al^. ** 

4 

unfortunate then, that the spectrochemlcal analyses which were 
obtained at the beginning of this Investigation proved to be 
Incorrect and that the results presented here were not obtained 
until the single crystal diffusion anneals and proton activation had 
been completed. 

4 

Herron Testing Labs, Cleveland, OH. 
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The boules were oriented by the Laue back reflection and sliced 
perpendicular to the c-axis. These samples were polished on a brass 
lap with successive grades of diamond paste suspended on Beuhler 
Texmet paper. Some rounding of the sample occurred during polishing. 
It was observed, by fortunate coincidence 5 6 , near the end of this 
work that using Buehler AB Selected Silk allowed one to take a rough 
ground sample directly from the watering machine to 3 ym diamond 
paste. This resulted In an enormous savings In time (approx 1 hr vs 
two days using 30 and 15 ym diamond discs and successive 6, 3, and 
1 ym diamond pastes). The sample was also much flatter by the time 
a scratch free surface was obtained. 

All samples were preannealed In air, at the temperature to which 
they would be exposed during the diffusion anneal. This was done In 
order to remove polishing damage and equlllbrlate the sample. The 
samples which were given diffusion anneals at 1300° and 1350* C were 
preannealed at 1600* C for 10 days eradication of polishing damage at 
1350* C would have required a prohibitively long time. 

3.1.2 Polycrystall Ine Samples 

Polycrystal line samples were prepared 6 from Al^Og powders which 
were supplied to us through NASA-lewls from Blakowskl. These con- 
sisted of five basic compositions, undoped, Cr-doped, Nl-doped, 


5 0n the advice of R. J. Schut. 

6 The assistance of F. W. Oynys In the operation of the C.W.R.U. hot 
press Is gratefully acknowledged. 
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Y-doped, and 2r-doped. The doped powders were of three levels, 200, 
500, and 1000 ppm. 

The powders were hot pressed at C.W.R.U. using a GCA/Vacuum 
Industries hot press. The billets were formed In a graphite die, 

3 In. diameter, at a temperature of 1500° C, under a pressure of 
4000 psl, for a duration of 1 hr. 

The spectrochemlcal analyses on the resulting material showed 
that the anticipated level of dopant was very nearly attained, with 
background Impurities ~250 ppm, (see table 5). 

TABLE 5. - IMPURITY LEVELS IN HOT PRESSED SAMPLES 5 (ppm) 


Sample 

Ca 

Cr a 

Cu 

Ga 

Fe 

Mg 

N1 a 

SI 

Zr a 

Na 

Y a 

C 

S 

Undoped 

30 



30 

3 

56 


88 


38 


31 


Cr-doped 

30 

210 


30 

7 

56 

— 

54 

— 

42 

— 

32 

21 

Cr-500 

30 

510 

- 

30 

19 

48 

- • - 

54 

— 

50 

— 

6? 

20 

Cr-1000 

30 

970 


30 

24 

48 

- - - 

56 

— 

38 

— 

16 

11 

N1-200 

30 

— 

20 

30 

3 

45 

200 

58 

— 

44 

— 

13 

21 

N1-500 

30 

— 

30 

30 

3 

56 

460 

65 

— 

42 

— 

7 

18 

N1-1000 

30 

— 

20 

30 

4 

52 

960 

70 

— 

59 

— 

7 

10 

Y-200 

30 

— 

- 

30 

5 

48 

- - 

70 

— 

51 

190 

30 

24 

Y-500 

30 

— 


30 

2 

52 

— 

88 

— 

53 

480 

4 


Y-1000 

30 

— 


30 

5 

56 

— 

54 

— 

53 

1000 

34 

26 

Zr-200 

30 

— 


30 

2 

56 

- - 

74 

195 

39 

— 

14 

18 

Zr-500 

30 

— 


30 

6 

59 

— 

70 

480 

38 

— 

4 

6 

Zr-1000 

30 

— 


30 

4 

64 

- • ~ 

56 

1010 

48 

— 

33 

18 


a Dopant . 

5 Coors Spectrochemlcal Laboratory. 


In the N1-500 ppm material, normally a green color, a sub- 
stantial number of "white spots" ~1 mm In diameter appeared through- 
out the volume. Under an optical microscope these "white spots" 

(see figs. 16 and 17) appeared as a core of very dense material sur- 
rounded by a halo of quite porous material. The origin of these 





Figure 16. - Optical micrograph of "white spot" In N1-500 ppm 
sample (40X). 



Figure 17. - Optical micrograph of "white spot" In N1-500 ppm 
sample ( 1 60X) . 
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spots Is unknown. X-ray diffraction of samples with and without 
spots were the same. EOAX showed an enrichment of Iron In the white 
spot compared to the matrix, however, the Iron content of the N1-500 
ppm was not different from the other samples (see table 1). If Iron 
Is the cause, there Is no reason why "white spots" should form In one 
material and not the others. 

One set of samples of each composition were subjected to a 
grain growth anneal at 1850° C for 1 hr. The samples were burled In 
powder, of the same composition. In a McOanel 998 crucible. It was 
observed that the undoped samples which had been In contact with the 
crucible were discolored, picking up the yellow tinge characteristic 
of McDanel 998. The resulting grain sizes were >0.1 mm for all sam- 
ples except the Nl-doped sample which had grain sizes <0.02 mm. 
Attempts to polish these samples were unsuccessful. This was attri- 
buted to mlcrc-cracklng upon cooling due to thermal stresses gen- 
erated between the large grains. 

A second set was subjected to a more modest anneal of 1630° C 
for 10 hr. These samples polished readily. The near surface den- 
sities and grain sizes for these samples are given In table 6. It 

should be noted that the density reported In table 6 Is not the bulk 

density, all samples registered densities greater than 99 percent by 
Archimedes' principle. Due to the very small penetration distance 
of the tracer, the bulk density Is not the measure of Interest, 

rather the density In the near surface (<1 mm). 


-62- 


T ABLE. 6. - PERCENT THEORETICAL 
OENSITY ANO GRAIN SIZE OF HOT 
PRESSED A1 2 0 3 



Sample 

Percent 

TO 

Grain size, 

pm 


Undoped 

97.00 

5 - 

10 

Cr-200 

98.00 

5 - 

15 

Cr-500 

97.00 

10 - 

20 

Cr-1000 

99.98 

10 - 

30 

N1-200 

96.00 

10 


N 1 - 500 

99.00 

5 


N1-1000 

96.00 

5 


Y-200 

97.00 

’ 50 - 

150 

Y-500 

99.98 

20 - 

60 

Y-1000 

95.00 

10 


Zr-200 

99.95 

10 - 

60 

Zr-500 

95.00 

10 


Zr-1000 

97.00 

10 



The reported densities were determined by comparing optical 


micrographs with 
density near the 
of grains during 
Examples of 


computer generated standards. The reduction of the 
surface Is doubtless contributed to by the pullout 
polishing. 

the microstructures are shown In figures 18 to 22 . 


3.2 EXCHANGE APPARATUS - 

A schematic of the apparatus In which the exchanges were per- 
formed Is shown In figure 23. It Is essentially the apparatus used 
by Reddy with some modlf Icatlons . The geometry of the shelves 
was changed to permit a greater number of samples to be exchanged 
side by side. The addition of a carbon reservoir allowed for a low 
p 0 to be obtained by the reaction: 
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Flgure 18. - Optical micrograph of undoped hot pressed AI 2 O 3 
sample annealed at 1630° C for 10 hr. (160X). 
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200 


500 


1000 


Figure 19. - 
at 1630° C 





Optical micrographs of Cr-doped AI 2 O 3 annealed 
for 10 hr (160X). 




200 ppm 


500 


1000 


Figure 20. 
at 1630° 



ppm 


ppm 




- Optical micrographs of Nl-doped AI 2 O 3 annealed 
C for 10 hr (160X). 
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500 ppm 



1000 ppm 



Figure 21. - Optical micrographs of Y-doped AI 2 O 3 annealed 
at 1630° C for 10 hr (160X). 
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500 ppm 




Figure 22 . - Optical micrographs of Zr-doped AI 2 O 3 annealed 
at 1630° C for 10 hr (160X). 
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1. FURNRCE 

2. SAMPLE SHELF 

3. CARBON RESERVOIR 

4. O-RING SEALS 

5. GRS SRMPLE OUTLET 

G. ZEOLITE 0-18 STORAGE 
7. THERMOCOUPLE VACUUM 
GAUGE 


0. FLOW METER 

9. DRIERITE 

10. PRESSURE/'VRCUUM GAUGE 

11. Zr0 2 CELL 

12. PUMP 

13. MILLIVOLT METER 


Figure 23. - Schematic of the exchange apparatus 
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3C ♦ 20 2 -» C0 2 ♦ 2C0 (3.2-1) 

The pump was added to circulate the gas and a Zr0 2 cell was placed 

In the system to monitor the P n . 

U 2 

3.3 EXCHANGE PROCEDURE 

The exchange for the 1 atm oxygen experiments was conducted In 
a manner similar to that employed by Reddy. The samples were placed 
on the shelf, the chamber was then evacuated, and the temperature 
was raised to ~1000° C, typically within thirty minutes. The tracer 
atmosphere was then allowed In, to a pressure slightly higher than 
that of the ambient In order to reduce contamination by air If a 
leak was present. The temperature was then raised to that desired 
for the diffusion anneal, typically In less than ten minutes. The 
temperature was controlled at ±5° C for the duration of the anneal, 
24 hours, shortest; 80 days, longest. 

In the bulk of the experiments, gas samples were taken periodi- 
cally to record the change In tracer concentration as a function of 
time. At the end of the anneal the furnace was shut off and the 
tracer atmosphere reclaimed to the zeolite storage pump. 

One exchange was carried out at low P Q . In this case the pro- 

2 

cedure was Identical to the one just described, except that a reser- 
voir of carbon was placed In the hot zone of the furnace. The con- 
figuration was such that the pump caused the gas to flow down the 
Inner tube over the carbon and across the samples. It was then re- 
circulated up through the outer tube. 
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Although the Zr0 2 cell suffered a loss of continuity during 
the experiment, by monitoring the pressure change as the gas equlll- 


brlated with the carbon, the calculated P n was verified. 

u 2 


3.4 ANALYSIS OF THE SAMPLES 

After the samples had undergone the diffusion anneal, the con- 
centration of oxygen-18 versus depth Information was obtained by two 
methods. Single Spectrum Proton Activation and Secondary Ion Mass 
Spectrometry (SIMS). Comprehensive treatments of both procedures 
are available^®* . Only the basic concepts and Items 

relating to our particular experiments will be discussed here. 

3.4.1 Proton Activation 

18 15 

This technique Is based on the nuclear reaction 0(p,o) N. 

In this reaction a proton striking an oxygen-18 nucleus causes an 
alpha particle to be ejected, yielding a nltrogen-15 nucleus. 

This work was carried out at the Unlverslte de Paris VII. The 
geometry of the setup Is Illustrated In figure 24. The samples were 
Irradiated with a monoenergetlc beam of 750 KeV protons. The 3 MeV 
alpha particles which resulted were detected with Ortec surface 
barrier detectors and energy discriminated with a Multichannel 
Analyzer, or a dedicated computer. The energy of an alpha particle 
Is related to the depth at which the reaction occurred and the 
Intensity of a alpha particles of a particular energy Is related to 
the concentration of oxygen-18. Ideally, the relationship of alpha 
energy to depth Is a unique function. This relation Involves the 
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Figure 24. - Geometry of the chamber used In Proton Activation. 
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rate at which Incoming protons and outgoing alpha particles lose 

energy. This Is given by the stopping powers which are 

. „ , . J57 to 59) 
tabulated 

The stopping power Is a function of the particle energy. Over 
the energy ranges encountered In single spectrum proton activation 
both the proton and alpha particles stopping powers could be accu- 
rately described by linear functions. 

The stopping power for protons In Al^, S p , over the proton 
energy range, E p , 750 to 500 KeV follows the relation 

S = B ♦ M • E (keV/pm) (3.4-1) 

P P P P 

with 

B . 161.637 
P 

M p = -9.512xl0' 2 

and the stopping power for alphas. So, over the alpha energy range, 

E , 3369 to 2500 keV follows 

a 

S * B ♦ M • E (3.4-2) 

a a a a 

with 

B = 582.710 

a 

M = -7.363xl0" 2 

a 

During a nuclear reaction relativistic energy Is conserved. If 

the assumption Is made that the 0-18 nucleus Is at rest then^ 6 ^ 

2 2 2 2 
m, Q c ♦ ♦ Ep = m., c ♦Etc ♦ m c ♦ E (3.4-3) 

18 0 p K 15 N 15 n o o 

where m Is the rest mass, c Is the speed of light, and E Is the 
kinetic energy. The energy balance, Q, may be written 
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7 2 

Q = (n ia ♦ rn - m., - m ) c = E, c ♦ E - E = Amc « AE (3.4-4) 

18q P 15 n a' 15 n a P 

If the value of Q Is positive, kinetic energy Is gained and the 
reaction Is termed exoerglc whereas when the value of Q Is nega- 
tive, kinetic energy Is lost and the reaction termed endoerglc. 

The ^®0{p,a)^N reaction Is exoerglc with Q = 3.98 MeV. 

Over the energy range used for proton activation the alpha energy Is 
linearly related to the proton energy by^ 1 ^ 

E = B + M Ep (3.4-5) 

a nr nr 

with 

B = 3002.724 
nr 

"nr ' °- 482 


Given the above, as well as: 

(1) A detected alpha particle energy, E^, which Is equal to 

a 

the alpha energy at the sample surface, E (x * 0) . 

Q 

(2) A known Incident proton energy. In our case 750 keV, which 
Is equal to E (x « 0). The relationship between the detected 
alpha energy and the depth at which the reaction took place, x r , 
can be Illustrated as follows. The geometry of the chamber Is such 
that the angle between the Incoming protons and the outgoing alpha 
particles Is 165°. The alphas, therefore, travel a distance 



(3.4-6) 
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Slnce S Is a linear function of energy, the effective 
<1 

stopping power, S a , between x » 0 and x = x r Is 
m 

S . B ♦ 5 s (E (x = 0) ♦ E (x = xl) (3.4-7) 
o o c a a r 

Substitution of equation (3.5-5) Into (3.5-4) and rearranging yields 


2 cos(15) «• M x 2 B x f 

E « (x r } * E o (x * 0) 2 cos(15) - M°x + (2 cos(15) - M x ) 

or or 


(3.4-8) 


The energy of the proton, E p (x = x r ), which under went the 
reaction to produce this detected alpha, Is related to E (x = x ) 

(k I 

by equation (3.5-3). It Is also related to the Incident proton 
energy by 


yv - E p < x * °) - Vr 


(3.4-9) 


By analogus substitution to that for alphas, this may be re- 
written as 


2 - Vr 


2 B x 




E p (x r ) s E p (x * 0) ’ 2 ♦ M_x_ 2 ♦ M.x 


(3.4-10) 


P r 


P r 


By combining equations (3.5-6), (3.5-3), and (3.5-8), and re- 
arranging the following may be obtained 


2 cos(15) - H x 
o r 


2 ' Yt 


E . <* ’ °> • ? cos( 1 5) 7 H , • B nr * £ p <* * °) • "nr TTltt 

or p r 


2 H nr B n X r 2 8 X r 

nr p r o r 

2 ♦ M p x r " 2 cos(15) ♦ 


(3.4-11) 


This equation contains as the only variables E (x = 0) and x . 

o r 

It Is a complicated expression, however, as the plot In figure 25 



E d - 3390.686-442. 1 16*X 
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Figure 25. - Plot of the energy of the detected alpha particle versus 
the depth at which the proton activation took place. Also shown Is 
linear regression of the depth energy relationship. 
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shows, the depth energy relationship varies by only a small degree 
from a linear relationship when the appropriate values of the 
constants are substituted In. A second order approximation was used 
In the spectra fitting routine. 

The number of counts In each channel Is related to the concen- 
tration of tracer by: 

* c(x^) Ltofx^ax^&n (3.4-1?) 

where 

= the number of counts In the 1th channel 
c(x 1 ) = the fractional concentration of tracer at a depth corre- 

sponding to the channel, x^ 

L * volume fraction of oxygen 

t = number of protons 

o(x^) = the cross section of the reaction for a proton at depth, 

X 1 

Ax^ ■ the width In distance of the 1th channel 
Aft a solid angle of observation 

Note L Is a material constant and i, o(x^), ax^, and Aft 
are experimental parameters. 

The number of protons was determined by measuring the charge 
Imparted to the sample chamber. A typical charge was 50 pcb corre- 

4 

spondlng to 3.1x10 protons. 

The cross-section Is dependent on the Incident proton energy 
chosen and Is Illustrated In figure 28. 
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(ajs/quj) -o 


Figure 26. - Differential cross section of 18 0(p t a)' 5 N reaction at 
scattering angle of 165*. 
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The width of the channel was selected electronically to corre- 
spond to ~13 keV which Is ~0.03 pm In distance. 

The solid angle Is dependent on the dimensions of the detector 
and Its distance from the sample. In our case. Aft - 0.4 Steradlans. 

This discussion has been for an Ideal case. In reality the 
monoenergetlc particles decay Into wider distributions as they pass 
through the material. The stopping powers give the average energy 
lost. The Imprecision In the detector and the electronics also give 
rise to uncertainties In the detected energies. Both of these 
effects can be accounted for by a spreading function. This function 
Is experimentally determined and found to be very nearly gausslan. 

It Is this spreading function which determines depth resolution. 

An example of a proton activation spectrum Is shown In figure 27. 

3.4.2 Rutherford Backscatterlnq (RBS) 

Many of the protons which strike the nuclei do not cause a reac- 
tion to take place and simply undergo an elastic collision. These 
Rutherford backscattered protons were detected simultaneously with 
the proton activation spectra. This was done to ensure that the 
proton beam was not being channeled. Channeling affects the proton 
activation spectra In two ways. 

First, channeled particles lose energy at a lower rate than 
those which see the lattice In a random way. This makes the energy 
depth relationship ambiguous. Two alphas with the same detected 
energy may have been generated at different depths, one from a 
channeled proton and from an unchanneled proton. 


COUNTS 
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Figure 27. - Example of a proton activation spectrum. 
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Secondly, channeled protons have a much lower possibility of 
reacting with the nuclei. This changes the relationship between 
alpha Intensity and concentration of tracer. 

The RBS spectra proved to be a good Indicator of channeling and 
were used as a guide to selectively reject spectra. 

The energy relationship of a backscattered particle Is: 

wv - wv <3 - 4 -' 3) 

where E RBS ( x r ) Is the energy of the backscattered particle; 
is the energy of the particle prior to backscatter Ing; and, K m Is 
the backscatterlng kinematic factor dependent on the mass and charge 
of the particle, and, the mass of the struck nucleus and the scatter- 
ing angle. 

The kinematic factor may be calculated using the following 
relation 62 


2 2 21/2 
m 1 cos e ♦ (m ? - m 1 sin e) 

m.j ♦ m ? 


( 3 . 4 - 14 ) 


where m 1 Is the mass of the particle, M 2 Is the mass of the struck 
nucleus and, e Is the laboratory scattering angle. 

The energy of a detected RBS particle, E^g^, also depends on 
the depth at which the RBS event took place. This may be written 


Vr 

E dRBS " V E o ‘ Vr* " cos( 1 5) 


( 3 . 4 - 15 ) 


where Is the depth of the RBS event, E Is the Incident partl- 
r o 

cle energy, and, S^ and S^ are the stopping powers for the parti- 
cle traveling Into and out of the material, respectively. 
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3.4.3 Secondary Ion Hass Spectrometry (SINS) 

In SINS the sample Is bombarded with a primary Ion beam, either 
Ar + or O' was used In our case. As the Ion Is Implanted It transfers 
Its momentum to matrix atoms. Near the surface sputtering may occur. 
The sputtered fragments may be elemental or compound and may be posi- 
tively charged, negatively charged, or neutral. Fragments from 
alumina may Include 0, Al, 0 ? , Al 2 , A10, A10 2 , A1 ? 0, etc. The 
fragments are then separated according to AMU/Charge by mass 
spectrometry. 

The bulk of this work was done at the Pennsylvania State 
University MRL (referred to as PSU) using a Gatan SIMS. In this case 
the primary beam was Ar + , typclally 7 keV. Due to the Insulating 
properties of alumina. It was necessary to compensate the positive 
charge of the primary beam by spraying electrons on the sample. 

These electrons also cause some desorption of light elements such as 
oxygen from the near surface region. This Electron Stimulated 
Desorption (ESO) gives rise to a background Intensity for both 
oxygen-16 and oxygen-18. The background was found to be constant, 
for single crystal samples, and could be subtracted out. 

The primary beam was rastered In order to give a large flat 
bottomed crater (see fig. 28). The detector was electronically gated 
to receive Information from only the center third of the crater. 

A quadrapole mass spectrometer was used on this machine. 

Additional work was done at the University of Illinois at 
Urbana MRL (referred to as UIU). In this case no electron beam was 
available, so the use of a negative primary beam was required. 


8 ? 




Cross section of SIMS crater 


Figure 28. - Profllometer trace of a crater sputtered with the 
PSU Gatan SIMS. 
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Unfortunately , the only Ion available on the machine was 0~, which 
complicates the Interpretation of the spectra. The beam was rastered 
and a double focusing mass spectrometer was used to separate the 
masses. 

The output from both machines was directly transferred to com- 
puters. This allowed monitoring more than one mass charge as a func- 
tion of sputtering time. Typical raw data for the two machines Is 
shown In figures 29 and 30. The conversion of the PSU data concen- 
tration profile Is straightforward. The fractional concentration of 
oxygen-18 tracer, c 1Q , at a given value of sputter time, t , for the 
case of the Ar* primary beam Is: 

c 18 (t s ) LS — 5 (3.4-16) 

VV * *H<V 

where h 18 (t $ ) Is the counting rate of mass-18 at time t $ , and 
h,,(t„) Is the counting rate of mass-16 at time t . 

IDS 5 

In order to convert sputter time to sputter depth, the total 
depth of the crater, measured with a profllometer was divided by the 

♦ 

total time under the beam to give the sputter rate, y; such that 
“s ■ 

In the case of the 0" primary beam, at UIU, the simple ratio 
does not hold because of dilution from the beam. The true concentra- 
tion of oxygen-18 Is given by: 
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Figure 29. - Example of the output of the PSU Gatan SIMS during 
depth profiling. 
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Figure 30. - Example of the output of the UIU SIHS during depth 
profiling. 
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C 18 (t s ) 1L 5 (3.4-1?) 

VV * WV ~ V beam) 

In order to determine the contribution from the beam, a sample 
analyzed with the Ar* primary beam was reanalyzed with the 0' beam. 
There were some differences In the data collection. The PSU machine 
counted for 1 second Intervals and plotted counting rate versus 
sputter time, whereas, the machine at the University of Illinois 
counted for ten second Intervals and plotted the number of counts 
versus channel. 

Determination of the sputtering rates allowed by the PSU data 
to be normalized to that of UIU. A plot of the two spectra Is shown 
In figure 31. Apparent In this plot Is a surface peak In the UIU 
data not found In the PSU data. Also, the lower values of UIU ratio 
arise from the primary oxygen dilution. The magnitude of the dilu- 
tion from the beam can be determined at each point and Is plotted In 
figure 32 As seen In this plot, the dilution from the beam makes 
up 80 percent of the total oxygen-16 counting rate and It Is not 
constant. Therefore, an attempt was made to normalize the UIU data 
to the PSU data by using the aluminum signals as the reference. Two 
assumptions were made (the notation used Is; the symbol h refers 
to the counting rate or counts, the subscript refers to the mass/ 
charge being counted, and the superscript refers to the location of 
the Instrument). 





Figure 32. 
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(1) the oxygen-18 signal divided by the aluminum signal was a 
constant for both machines. These ratios would then be propor- 
tional, that is: 


h uiu /h uiu 

- a constant = 5.3 (see fig. 33) 

h 18 /h 27 


(3.4-18) 


(ii) the total oxygen over aluminum signal for PSU is constant, 
that is: 


.PSU //K PSU .PSU. 

Na /(ll ia 1 h u > 

.PSU ..PSU 
h 18 /h 27 


= a constant * 11.9 (see fig. 34) (3.4-19) 


Notice these are not perfect assumptions. Using these assumptions: 


C 18 <V * 


.UIU .PSU ..PSU 

h 18 h 18 /h 27 

.UIU .UIU..UIU 
"l 3. 5 n 18 '"l3.5 


.PSU... PSU .PSU, 

h 18 /(h 18 * h 16 > 

.PSU,. PSU 
hi 8 /h 2? 


(3.4-20) 


which, when values are substituted in, becomes 


.uiu 

11.9 n 18 


5.3 .UIU 
h 13.5 


C 18 <V 


s 


(3.4-21) 
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Figure 33. - The ratio of tracer signals, normalized to aluminum, between 
PSU and UIU. 
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CHAPTER IV 


RESULTS 


4.1 PROTON ACTIVATION 

4.1.1 Standards 

In order to obtain experimental values for spectrum fitting 
parameters, four types of standards were utilized: 

(I) Electronic test pulses which were generated and run through 
the electronics to the amplifier on the detector and back to the MCA. 

(II) Two oxide films of uniform oxygen-18 concentration. A thin 
film, 10 nm, of Al^ on anodlcally oxidized aluminum; and a thick 
film, 600 nm, of SIO^ on thermally oxidized silicon In a water bear- 
ing atmosphere. 

(III) Metallic samples of niobium and tantalum. 

( 1 v ) A sample of the undoped no. 1 boule containing only the 
natural abundance of oxygen-18. 

As previously mentioned, the proton activation spectra was 
spread due to the Imperfect nature of the electronics and the de- 
tector. In order to assess the magnitude and source of the spread- 

18 

Ing, the proton activation spectrum for the thin film of A1,,0 3 and 
an electronic test pulse were analyzed. The actual spectrum Is shown 
In figure 35. The data points are Indicated with crosses and the 
curves drawn through them are Gaussian Functions. In the case of the 
electronic test pulse, the standard deviation was 12.4 KeV (~30 nm) ; 
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Figure 35. - Proton activation spectrum on a 10 nm thick film of 
AI 2 O 3 18 and an electronic test pulse. 
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and In the case of the Al^O^ 8 * he s ^ an ^ ar ^ deviation was 25 KeV. 

The electronics appear to contribute to nearly half of the spreading. 

It was also necessary to calibrate the MCA In terms of the 

energy of each channel. Three physical calibration points were used. 

RBS spectra were obtained, on the metallic samples of niobium and 

tantalum, using alpha particles with an Incident energy of 1850 KeV. 

This gave points near the low end of the energy range of Interest. 

A point at the high end was obtained through proton activation of the 
1 8 

thick film of SI 0^ • These spectra are shown In figures 36 to 38 
and the surface channels are Indicated. Since these are surfaces 
only the kinematic scattering factor needs to be known for the RBS 
and the reaction kinetics for the proton activation, there Is no 
ambiguity Introduced by stopping powers. 

A set of electronic test pulses, (see fig. 39) which were gen- 
erated by passing a constant pulse through an amplifier set over a 
range of decibels, was calibrated using the physical values of the 
expression 


-db = 20 log 



10 V 


Input 


A set of such pulses was run at the beginning and close of each day 
of analysis, by referring to the calibrated set changes In the 
setting were followed. 

The MCA was found to deviate significantly from linearity and a 
second order polynomial was always used, a typical channel vs energy 
plot Is shown In figure 40. 
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Figure 36. - Physical calibration point no. 1. Rutherford backscatterlng 
spectra on Niobium, E surf = 1561.6 keV. 
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Figure 37. - Physical calibration point no. 2. Rutherford backscatterlng 
spectra on Tantalum, E 5urf = 1695.9 keV. 
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Figure 38. - Physical calibration point no. 3. Proton activation on 
6000 A thick slab of 0-18 enriched S10 2 , E surf = 3369 keV. 
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The points are from the electronic test pulses and the curve Is a second 
order polynomial fit. 
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Uslng these values and the tabulated stopping powers, assuming 
the elemental values could be added to give molecular values, an 
attempt was made to fit the proton activation spectrum for the 
natural abundance. The result Is shown In figure 41. It Is lacking 
In three separate ways. The location of the resonance Is off, the 
resonance Is too sharp, and the overall number of predicted counts 
Is low compared to the experimental points. 

The mlslocatlon of the resonance Is felt to arise from errors 
In the tabulated stopping powers (quoted errors are 10 to 20 percent) 
and additional error Introduced by assuming additivity. Based on 
the fact that the energy-depth relationship Is very nearly linear 
the decision was made to multiply this relationship by a factor 
sufficient to bring the predicted and experimental resonances coin- 
cident. The result was a difference of 13 percent. 

The Increased width of the experimental resonance was attri- 
buted to straggllng^ 63 ^ . This arises from the statistical fluctua- 
tions In the number of nuclear collisions which take place per unit 
distance. This may be approximated by a Gaussian function whose 

standard deviation Is proportional to the square root of the distance 

1 /2 

travelled. The resonance was fit using a value of 4 KeV/^m 

The overall low counting rate was attributed to the use of the 
data of Amsel et al. 1 for the reaction cross section. This data 
and the approximation used In the fitting program are shown In fig- 
ures 42 and 43. In previous work Y1nnon^ 3 ^ found these values 

1 Kindly supplied by L. D. Major. 
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Figure 43. - Close-up of the resonance In the cross section for the nuclear 
reaction ® 0 (p,<*) ':>N t showing the quality of the Gaussian approximation. 
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had to be multiplied by a factor of 1.2 In order to obtain a satis- 
factory fit for a spectrum from a sample of Ta^O^ enriched to 98 per- 
cent oxygen-18. In this work a multiplication of the Amsel data by 
a factor of 1.29 gave a good fit. 

A fit of the natural abundance with these corrections Is shown 
In figure 44. The corrections were used In all subsequent fits of 
concentration profiles. It Is clear that an Internal landmark In the 
form of the resonance Is quite useful. 

4.1.2 Channeling 

The channeling of protons by A1 ? 0 3 was Investigated primarily 
to assure It would not occur because of Its deleterious effect on the 
Information In a proton activation spectrum. 

Two types of channeling experiments were performed. First, the 
axial channeling positions were determined. In this experiment the 
samples were placed In a gonlometrlc chamber. In the chamber It Is 
possible to translate the sample and/or rotate It around 3 axes. 

The alignment was performed using a beam of 1.5 MeV alpha parti- 
cles rather than 750 KeV protons because greater mass separation Is 
obtained. 

Rutherford backscatterlng spectra taken In both a "random" 
direction and when the crystal Is oriented such that the beam travels 
down the c-axis, are shown In figure 45. The term random arises 
because, due to the thermal vibrations of the lattice atoms and the 
probabilistic nature of the electron positions, the location of an 
unchanneled particle, after a small number of collisions, can only 
be described statlstlcly. 
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compared to the experimental points. 
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it can be seen that In the channeling orientation the counting 
rate for collisions with the lattice Is markedly reduced both In the 
A1 and 0 plateaus. The counting rate for the amorphous carbon layer 
evaporated on the sample to dissipate charge Is, of course, not 
affected by orientation. 

A second sample from the same boule which had been polished but 
not subjected to any anneal, was treated In the same way. The RBS 
spectra for random and c-axis orientation are shown In figure 46. A 
comparison of figure 45 with figure 46 shows that the counting rates 
In the channeled spectra for the unannealed sample are considerably 
higher than for that of the annealed sample. Note, also, that the 
counting rate decreases. Indicating an Increase In channeling, until 
channel 210, corresponding to a depth of 0.25 yin. This Is certainly 
due to the damaged layer near the surface of the sample which was 
generated during the mechanical polishing. The absence of this be- 
havior In the annealed sample Is experimental evidence that this 
damage can be removed by annealing. 

A second, polished and annealed, sample was Investigated by 
transmission electron microscopy (TEN) . The sample was back thinned 

to reveal the region which had been subjected to the polishing. An 
-4 2 

area of 1x10 cm was Investigated and no evidence of damage was 

seen, e.g., dislocations, figure 47. A material with a dislocation 
6 2 

density of 10 lines/cm It would be expected to find on the order 
of 100 lines within this area. This serves to confirm the RBS data. 
This result Is not new. Both Hockey^ 63 ^ and Reed^ 3 ^, found that 


0009 



OOSZHtl 


Figure 46. - Rutherford backscatterlng spectra on an unannealed a-Al 203 
sample. In and out of axial channeling conditions. 
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and annealed a-Al 
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the damage produced by abrasion, In A1 2 0 3 , could be removed 
through annealing. 

A third set of spectra were taken from a different location on 
the sample which had been annealed at 1600* C for 300 hr. These are 
shown In figure 48. It can be seen that there are counts at higher 
energies than those arising from aluminum at the surface. There Is 
both a plateau, which Is depressed when the crystal Is aligned In a 
channeling orientation, and a sharp peak, which Is unaffected by 
channeling. The energy of the peak corresponds to 1.5 MeV alphas 
scattered from nuclei ~ 39 to 40 Amu which corresponds to K and Ca. 
The fact that the peak remains unaltered when the sample Is In a 
channeling orientation suggests that these cations are located In an 
amorphous (at least not epitaxial) phase on the sample surface, and 
the plateau Is consistent with diffusion of these cations Into the 
sapphire. The counting rate for the heavy Impurities would be ex- 
pected to behave during channeling like that of A1 when the Impuri- 
ties occupy A1 sites . 

This suggests that the sample was contaminated at some point 
during the exchange procedure. It Is also Indicated that the con- 
tamination Is not uniform across the sample. Later work established 
that the samples were Indeed contaminated In the form of a Ca, K, 
A1, SI glass on the sample surface. 

The relative concentrations of the heavy Impurity may be esti- 
mated from the plateau heights In figure 48. The only variable 
between the leading edges of the two plateaus Is the cross sections 
for the two species. The ratio of cross sections can be calculated 
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Figure 48. - Rutherford backscatterlng spectrum on a-Al203 annealed at 
1600° C, containing a heavy Impurity. 
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from the expressions^) 
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where o Is the differential cross section, E Is the particle 
energy, Z, and M, are the atomic number and mass of the particle and 
l 2 and are the atomic number and mass of the struck nucleus, 
while e Is the laboratory scattering angle. 

When the appropriate values are substituted the ratio of Impu- 
rity to aluminum cross sections, R c> Is very nearly equal to the 
ratio of masses 

, !Hi_ . I? 

R c ■ m AL * 27 • '- 5 


Extrapolating the Impurity plateau to the surface, and measuring the 
heights In figure 48 allows the Impurity concentration, c Jt to be 
estimated by 


C 


I 


H I 4 (arb. units) 

1.5 H A1 ♦ Hj = (1.5) • (5) ♦ 4 


35 percent 


This Is a very high concentration. It Is such a high level that It 
Is surprising that the material still has A1 2 0 3 channeling behavior, 
although McCune^ 66 ^ has observed Ca concentrations of 10 percent In 
the near surface region arising from the segregation of a dopant 
level (~100 ppm) bulk concentrations. This sample was the only one 
In which this very high level of Impurities was observed. Concentre- 
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tlons up to a percent, however, would not produce a counting rate 
sufficient for detection. 

The second type of channeling experiment employed an Incident 
beam of 750 KeV protons and two detectors. One detector was located 
on either side of the sample making the same angle to the beam. 

Each detector was connected to a separate set of electronics, one 
was used to collect a RBS spectra and the other was used to collect 
a proton activation spectra. In this way the effect of channeling 
orientation on each could be studied. 

To check for planar channeling, the sample was oriented to 5* 
away from the c-axis. The sample was rotated about the beam In three 
degree steps. The geometry Is schematically Illustrated In figure 

49. At each step the sample was bombarded with the proton beam 

18 15 

spectra both RBS and 0(p,a) N were taken. The Integral of the 

18 15 

aluminum plateau of the RBS spectra and of the entire 0(p,a) N 
spectra are plotted as a function of the angle of rotation In figure 

50. Two points should be noted here. The first Is that when the 
beam Is nearly parallel to a principle crystallographic axis (11.0), 
(10.0), or (01.0), the counting rate for the proton activation spec- 
trum Is less than half of what It Is In a random position. In be- 
tween these major drops there appear smaller decreases ~5 to 10 per- 
cent. Secondly, the RBS A1 plateau Is much more sensitive to the 
channeling and can be used as a guide for acceptance or rejection of 
a set of spectra. 

For the normal course of spectra accumulation In the chamber 
(shown In fig. 24), both the RBS and the proton activation spectra 



Figure 49. - Schematic of the geometry for a rotational scan. 
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Figure 50. - Results of the rotational scan showing the presence of 
planar channeling In AI 2 O 3 . 
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were accumulated. If the RBS Indicated channeling, the run was 
aborted, the sample rotated a few degrees and rerun. Channeling did 
appear periodically during analyses, and those spectra aborted. The 
RBS on A1 was not a perfect Indicator however, due to the differ- 
ences between the crystal structure of the aluminum and oxygen 
sublattices. 

Orientation effects occurred In the RBS A1 plateau which were 
not carried over Into proton activation spectra. The two sets of 
figures 51, 5 2 and 53, 54, show simultaneously recorded RBS and 
proton activation spectra. Note that although the RBS A1 plateau 
has either a depressed, random, or enhanced counting rate, depending 
on orientation, the proton activation spectra for each set may be 
superimposed. For this reason a number of spectra were taken on each 
sample to ensure that “random" Information was obtained. 

An example of the effect of being In a random orientation, 
planar channeling, and axial channeling Is shown In figure 55. The 
effect on the counting rates and the position of the resonance can 
be seen. 

4.1.3 Nuclear Reaction 

4. 1.3.1 Qualitative Interpretation 

The samples analyzed by proton activation were prepared and ex- 
changed In groups. One sample from each of the five boules and 
usually a duplicate of one were exchanged sitting side by side on the 
same shelf. In some cases several shelves were loaded at different 
points In the thermal gradient, to simultaneously exchange at differ- 
ent temperatures . The spectra for the three undoped boules exchanged 
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The A1 plateau displays an enhanced counting rate at -4° of rotation. 
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Figure 52. - A comparison of proton activation spectra taken on Nl-12 
corresponding to the R8S spectra In figure 51. 
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plateau displays a depressed counting rate at 0° of rotation. 




1000 


- 120 - 


Q 

Q 

OJ 



G> 

in 


CP 

( s 

o 


© 


Q 

o 

G3 


© 

Q 

CD 

U5 

'T 


AJ 



s;uno3 


Figure 54. - A comparison of proton activation spectra taken on UD-6 
corresponding to the R8S spectra In figure 53. 
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ure 55. - The effect of axial and planar channeling on a proton activation 
pectrum. 
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at 1458° C, Is shown In figure 56. Because of the similarity of them 
only one undoped spectrum was plotted on the rest of the plots. 

Several spectra were taken on each sample. An example of two 
spectra taken one after the other on the same sample, with a four 
degree rotation of the sample between collection. Is shown In figure 
57. The agreement between them Is quite good but there Is some data 
scatter at the peak. 

An example of two samples from the same boule, exchanged 
together. Is shown In figure 58. Again, the spectra look very 
similar although In this case the "tall" of N161 appears to be a 
little higher than that of N160. 

figures 59 to 62 show spectra for samples of the three composi- 
tions plotted for four exchanges at temperatures ranging from 1415° 
to 1677° C. The most striking feature of each plot Is the sameness 
of each profile. There Is a small amount of scatter but It Is no 
larger than that between spectra on the same sample or between sam- 
ples from the same boule. 

Another feature of these spectra Is the deeply penetrating 
"tall" of low concentration. "Tails" of this type were observed by 
Reddy and Cooper. They claimed to be able to remove these "tails" 
by preannealing their samples In air. Our samples were subjected to 
such preannealing treatments, however the "tails" persist. There Is 
some difference between the "tall" magnitudes but It does not appear 
to be systematic. For example, In figure 60, the Nl-doped sample 
has a larger "tall" than the undoped sample, yet In figure 59, the 
situation Is reversed. 
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Igure 56. - Proton activation spectra for samples from the three undoped 
Doules, exchanged at 1458° C for 404 hr 
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Figure 57. - Two proton activation spectra taken on the same sample 
exchanged at 1458° C for 404 hr. 
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same boule, exchanged together at 1458° C for 404 hr. 
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Figure 60. - Proton activation spectra taken on an undoped sample, a Cr-doped 
sample, and a Nl-doped sample all exchanged together at 1458° C for 404 hr. 



Figure 61. - Proton activation spectra taken on an undoped sample, a Cr-doped 
sample, and a Nl-doped sample all exchanged together at 1634® C for 98 hr. 
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Figure 62. - Proton activation spectra taken on an undoped sample, a Cr-doped 
sample, and a Nl-doped sample all exchanged together at 1677° C for 23 hr. 
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A third feature is that all the experiments up through 1550° C 
have tracer penetration depths which are small compared to the re- 
solution of the analysis. That Is the concentration profiles may 
contain structure that cannot be seen. 

It was also noted that the spectra of Cr had apparent counts of 
an energy higher than that corresponding to the surface of the 
sample, unless a very low beam current was used. If the charging 
from the proton beam was not being dissipated by the carbon film. It 
would continue to build until It reaches the breakdown strength of 
the sample. Such an arcing condition could affect the spectra^ 67 ^. 
The distortion could be accounted for electrical breakdown If the 
ruby were slgnlf Icantly more resistant than the other samples, 
however the reported values for conductivity of ruby are very near 
that of sapphlre^ 58 ^ , or If carbon layer was different and not 
conductive when deposited on ruby. The latter does not seem likely. 
There was no visually observable difference In the carbon layer 

i 

between samples and all the samples were coated together, and the 
samples remained at room temperature. The detector Is sensitive to 
light however, and It seems possible that the distortions of the 
spectra (see figures 63 and 64) which Involve a shift of counts to 
higher energy, might occur If the detector were to count a large 
number of 2 to 3 eV photons and an alpha simultaneously. While all 
samples fluoresced a violet color, the Cr-doped sample had a much 
more Intense bright red color, typical of ruby. 

The experiment necessary to prove this assumption, would be to 
place an opaque piece of aluminized mylar between the sample and the 
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Figure 63. - Proton activation spectra taken on Cr-doped sample with a beam 
current of 210 nA and 120 nA. 
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detector, however. It was not possible to return to Paris to attempt 
the experiment. 

If Indeed the fluoresence of the ruby causes the spectra to be 
altered, then It Is also possible that the fluorescence of the other 
samples might cause this shift of counts to higher energies. It must 
be emphasized that this Is speculation supported by no experimental 
evidence. 

These concerns and the logistics of traveling to France to 
analyze samples prompted the switch to SIMS for the remaining sample 
analysis. 

4. 1.3. 2 Fits to Diffusion Equations 

Many of our spectra were shallow compared to the depth resolu- 
tion of the analysis technique. However, reasonable fits could be 
obtained. Fits, calculated using equation (2.1-7), for three tem- 
peratures (1457° C. 1550° C, and 1677° C) are shown In figures 65 to 
67. It can be seen In these figures that the magnitude of the tall 
(another way to say this Is; the extent of deviation from the 
predicted profile at large depths) Increases with Increasing 
temperature. 

A fit of the tall for the specimen exchanged at 1677° C Is shown 
Irt figure 68. 

It was also found that the value of the surface exchange co- 
efficient, K, was sensitive to fluctuations In counting rate between 
spectra while the value of the diffusion coefficient was not. The 
spectrum shown In figure 69 was taken consecutive to that shown In 
figure 67. The spectrum In figure 67 was best eye fit assuming the 
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gure 67. - Proton activation spectrum, with a near surface fit, from 
sample UO-69 exchanged at 1677° C for 23.3 hr. 
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gure t>9. - Proton activation spectrum, with a near surface fit, from 
sample UD-69 exchanged at 1677° C for 23.3 hr. The value of the value 
for the surface exchange coefficient Is very sensitive to data scatter 
while the value for the diffusion coefficient Is not. 
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surface reaction was rapid enough to play no role, l.e., K = «, 

-19 2 

with a value of 0 = 4.5x10 m /s while the spectrum In figure 69 

-12 -19 2 

was best fit by assuming K = 8x10 m/s and 0 = 5.0x10 m /s. 

t 1/2 

In the later case K (^) = 3.3, which suggests that the surface 

reaction plays a significant role. Note, however, that the diffusion 
coefficients vary by only ~10 percent. Also, the spectra In gen- 
eral were less different than this example pair, and the variation 
In determined diffusion coefficient accordingly smaller. 

There are three possible causes for the Inability to fit the 
spectra. These are: 

(I) a non-FIcklan profile 

(II) an error In the computer calculated spectrum causing It to 
underpredict the concentration at larger depths, and 

(III) an artifact from the proton activation experiment 

The computer program was extensively checked and found to be 
correct. Also, latter work with SIMS found that profiles from these 
same samples confirmed the exlstance of the tails and their relative 
magnitudes. This observation by two very different methods of analy- 
sis greatly diminishes the probability that they are an artifact. 

A summary of the samples analyzed by proton activation Is given 


In table 7. 
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4.2 SECONDARY ION MASS SPECTROMETRY 

4.2.1 Oxvqen-18 Profiles In Single Crystals at P. =1 atm. 

2 

SIMS has a much better depth resolution, ~10 A, than does proton 
activation. It was warranted to recheck the earlier observation that 
all the samples exchanged together have Identical tracer profiles. 
Shown In figure 70 Is a plot of the concentration of tracer as a 
function of depth for two samples analyzed at Penn. State, which 
were exchanged at the same time, one undoped and other Nl-doped. 

The two curves are Indistinguishable. This held true for every case 
that was examined, although this was a limited number of samples. 

The second finding of the SIMS analysis was that, as with proton 
activation, the tracer profile could be divided Into two regions. 

A near surface profile, which could be fit with a solution to the 
diffusion equation, and a deeply penetrating tall. Data from the 
Penn. State SIMS work and fits, for the temperature range 1300° to 
1677° C, are shown In figures 71 through 76. 

As In proton activation, deep penetrating "tails" appear to 
varying degrees In the profiles. Again, they do not appear to be 
systematic. The samples UO-A and UD-B were preannealed together and 
exchanged, on different shelves, during the same run. Sample UD-A 
shows a pronounced "tall" and UD-B does not. 

The data from the University of Illinois SIMS work, despite 
utilizing the Imperfect normalizing coefficients, were able to be 
fit quite well. This data (over the temperature range of 1330° to 
1458° C) and the fits are shown In figures 77 to 79. The "tails" 
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Figure 70. - PSU SIMS profiles on an undoped sample and a Nl-doped 
sample exhanged at 1435° C for 1920 hr. 
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Figure 72. - PSU SIMS profile on an undoped sample UD-8 exchanged at 1415 
for 1920 hr. 
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PSU SIMS profile on an undoped sample UD-69 exchanged at 1677° 
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whlch appear In these profiles may well arise during the normaliza- 
tion procedure used to eliminate the beam contribution, and should 
not be taken too seriously. Evidence for this Is seen In the fact 
that the sample UO-B shows a "tall" In the UIU profile and not In 
the PSU profile. 

The one exception was the experiment at 1634° C which could not 
be fit with a single pair of 0 and K. The PSU profile on this 
specimen are shown In figure 80 (attempts to fit the proton activa- 
tion data on this sample were also unsuccessful). This experiment 
was run with new furnace tubes and new tracer. The starting con- 
centration of the tracer was 95 percent and It dropped to ~18 per- 
cent at the end of the run. Unfortunatenly the rate of change of 
the tracer concentration Is unknown since fact that all but two of 
the gas samples, both taken near the end of the run, were broken In 
the mall on the way to be analyzed. 

An Arrhenius plot of the diffusion coefficients determined In 
this work Is shown In figure 81, along with a best fit line to the 
PSU SIMS data and a summary of the previous work. Note that the 
data agrees well with the previous gas exchange data, particularly 
with that of Reddy and Cooper and may be described by: 

0 = 1.51xl0' 3 exp m 2 / s (4.2-1) 

Note that the values obtained from proton activation are around a 
factor of two lower than those obtained from SIMS for the same sample. 
The proton activation fitting procedure Involved a number of assump- 
tions (Sec. 4.1.1) for the depth-energy relationship while with SIMS 
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Figure 81 
coefficient 


An Arrhenius plot of the determined oxygen self diffusion 
s In A 1^03 with an edited summary of previous work. 
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the only assumption was that the sputtering rate was linear. It Is 

possible that the profllometer used to measure the crater depth was 

out of adjustment, however comparison of results obtained by Inter- 

2 

ferometry compare well with those from the profllometer. Therefore, 
the SIHS results are felt to be more accurate. 

4.2.2 Oxygen-18 Profiles In Single Crystals at P n « 10~ 15 atm. 

y- 

The tracer profiles In the samples exchanged In a CO/CO^ atmos- 
phere, equivalent to P n = 10~ 15 atm., looked very similar to those 

2 

all four samples; two undoped, one Cr-doped, and one Nl-doped, looked 

the same. The profile of an undoped sample Is shown In figure 82. 

Similar to the P ft = 1 atm. profiles, there Is also a near surface 
°2 

region that appears Flcklan and the presence of a low concentration, 
deeply penetrating "tall". 

The value for D Is plotted on the ^Arrenhlus plot (fig. 81). 

Note that It fits very well on the line determined by the P n = 1 atm 

u 2 

data, Indicating no P Q dependence for oxygen diffusion In Al^Og. 

A summary of all the samples analyzed by Secondary Ion Mass 
Spectrometry Is given In table 8. 

4.2.3 Oxygen-18 Profiles In Polycrvstals at P . =1 atm. 

The determination of oxygen profiles In the polycrystalllne sam- 
ples, by SIMS, proved to have complications. Recall, from Section 
3.5.3, that It was necessary, with the PSU SIMS, to have an electron 

2 

Personal Communication Oavld Roach. 
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beam coincident to the Ar+ beam to avoid charging. In the case of 
single crystals this produced a low level (see fig. 29) constant 
background which could be subtracted out of the data. With the poly- 
crystals the electron sputtered desorption (ESO) could not be dealt 
with In this simple way. 

Shown In figure 83 Is a SIMS profile on a sample of undoped, 
hot pressed, A1 0 . This sample was analyzed consecutive to a single 

c O 

crystal sample which had a reasonable 25 percent ESO, and the 
settings on the SIMS unit were unaltered between the two. The curve 
for sputter time less than -1280 sec represents the count rate under 
normal profiling, that Is when both the argon Ion beam and the elec- 
tron beam are striking the sample. The curve for sputter time 
greater than -1280 sec represents the count rate when only the elec 
tron beam Is allowed to strike the sample. 

There Is a sharp depression when the argon beam Is stopped, 
however the signal rapidly returns to very near the level that was 
obtained during profiling. Apparently the vast majority of the Ions 
detected were generated through grain boundary etching by the elec- 
tron beam. This behavior was observed on every polycrystalllne 
sample which was analyzed. 

By reducing the argon beam current to one-third of Its normal 
value (1.5 yA to 0.5 yA) It was possible to lower the electron beam 
Intensity and collect data, although the ESO still produced up to a 
third of the counts. The resulting profiles, for undoped and singly 
doped (1000 ppm) with Cr, N1 , Y and Zr, are shown In figures 84 to 



cocoot 



SINHOO 


Figure 83. - PSU SIMS profile on a sample of undoped hot pressed Al 
showing the high level of electron sputtered desorption (ESD). 
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88. Some charging problems occurred on the Zr-doped sample and the 
near surface region was extrapolated. 

A qualitative comparison between the profiles can be made using 
the heavy smooth line that appears In each figure. This line was 
obtained by drawing a smooth curve by eye through the un doped data. 
Translation along the concentration axis was allowed when the line 
was reproduced on the other profiles. It can be seen that this line 
Is also a good representation of the data for the Cr-doped, Nl-doped, 
and Zr-doped samples. The profile from the Y-doped sample, however, 
appears to be different. 

It Is not known If this Is a real result or an artifact. This 
Y-doped sample was found to be physically stuck to the sample shelf 
after the exchange. This was the only one, of the thirteen samples 
on the shelf, which stuck. Although no clear connection between this 
and the oxygen profile Is seen It may be that this sample was not a 
good sample. Also, the ESO was not rechecked to be sure It was 
Independent of sputter time. If the rate of etching of the grain 
boundaries changes (which Is dependent on microstructure) or the 
concentration of tracer In the boundaries changes during sputtering 
the shape of profile may be altered. 

In summary, the diffusion of oxygen In polycrystal 1 Ine Al^O^ 
appears Insensitive to 1000 ppm of Cr, N1, or Zr but may be affected 
by 1000 ppm of Y. However, the analysis method was not optimized 
for this type of sample and the results are not conclusive. Oue to 
this uncertainty quantitative fitting of the profiles was not 
attempted . 
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Figure 85. - PSU SIMS profile on a sample of Cr-doped hot pressed AI 2 O 3 . 
The smooth curve Is the reference from the undoped profile. 
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Figure 87. - PSU SIMS profile on a sample of Y-doped hot pressed 
The smooth curve Is the reference from the undoped profile. 
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The smooth curve Is the reference from the undoped profile. 


-166- 


4.2.4 Attempt to Profile Polycrystals at Low P n 

J 2 

Samples of polycrystalline material were subjected to a low P 

2 

preanneal (In CWRU Hot Press) and a low P n diffusion anneal. 

u 2 

However, due to severe surface face roughening no tracer profiles 
were determined. 

The surface roughening Is Illustrated with profllometer traces 

of four samples shown In figure 89. It can be seen that the surface 

roughening In a low P environment Is considerably greater than In 

°2 

1 atm. P_ , for the same time and temperature. 

U 2 

The prospect of determining grain boundary diffusion coeffi- 
cients In low P n , by this method Is dim. In order to get a profile 
°2 

which Is large with respect to the surface roughness, would require 

tracer penetration and subsequent sputtering to depths >20 pm. 

Notice that with the 1 atm 0 2 sample, It would be satisfactory to be 

at a depth of 1 pm. Also, the single crystal sample retained Its 

finish despite the low P_ . 

u 2 

4.2.5 Contamination of Single Crystals 

In tracer diffusion studies the assumption Is made that the 
chemistry of the material which contains the tracer profile Is char- 
acterlstlc of the chemistry of the bulk. In a gas exchange experi- 
ment, particularly when the diffusion coefficient Is as low as In 
alumina, this region Is only a few microns Into the surface. 

In order to check the validity of this assumption, the near 
surface region of some samples were mass analyzed by SIMS. Typical 
mass scans from the UIU SIMS and the PSU SIMS are shown In figures 
90 and 91. Differences appear between mass scans taken with the two 


machines . 
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SURFACE ROUGHNESS OF HOT PRESSED ALUMINA 
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Figure 89. - Surface roughness of hot pressed alumina 
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Figure 90. - PSU mass scan, typical of the PSU Instrument, on an undoped 
as-received sample. 
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The scan from the UIU camera shows a greater number of peaks and 

has better counting statistics (note that log counts Is plotted for 

UIU and counts for PSU) . The Increased counting rate probably 

1 551 

arises from oxygen enhancement.' 

A second, more puzzling difference Is the presence, In the PSU 
scans, of peaks for sodium and potassium which do not appear In any 
of the UIU scans. Since sodium was detected In the spectrochemlcal 
analysis (table 4), It Is surprising that It does not appear In the 
UIU scans. Potassium, on the other hand, was not detected In the 
spectrochemlcal analysis and yet It appears In the mass scans from 
Penn. State. The Isotopes of potassium appear In the correct ratio 
AMU 38.96 (93 percent) and AMU (40.96 (7 percent). The argon primary 
beam Is known to be a source of potassium however the beam contribu- 
tion should be constant and gradients were found. 

Fingerprint contamination Is a possible source of these ele- 
ments. However, If such were the case, the concentration of sodium 
would be expected to be tied to that of the potassium, and this Is 
not the case. Also, the mass scans were taken after sputtering 
through the carbon layer and Into the sample at greater than 0.1 pm 
which should be well past surface contamination. It may be that the 
difference In primary beams Is the source of this discrepancy. 

In order to assess the contamination of our samples, mass scans 
were taken on the as-received boules during and after the annealing 
process. All of the peak heights were normalized to the height of 
the Al* + peak. This allowed a comparison from sample to sample. 
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Mass scans of two undoped, as-received, samples; one from undoped 
no. 3 and one from material used by Reddy, are shown In figure 92. 
The two appear very similar. Note the presence of a little sodium, 
some silicon, and an Indication of potassium. Also, the counting 
rates of the other fragments of AI^O^ show some fluctuation. 

This fluctuation, plus the knowledge that various chemical species 
have different sputtering yields, dependent on their environment, 
cautions against quantitative Interpretation of these scans. 

However, the qualitative comparisons show a distinct trend. 
Figure 93 shows scans for two samples from the undoped no. 3 boule, 
one as-received and one after an anneal at 1550° C. The difference 
Is striking. After the exchange the sample picked up magnesium and 
calcium which were not present before. Also, peaks for silicon and 
potassium are distinctly larger. 

Mass scans on a sample exchanged at 1410° C and a sample of as- 
received material, are shown In figure 94. The level of contamina- 
tion Is significantly less than the sample exchanged at 1550° C, 
despite the fact that that Dt's for oxygen were similar. There Is 
an Increase In the sodium peak and a magnesium peak appears after 
the exchange, while the apparent concentration of silicon has de- 
creased. The Indication that potassium Is off scale stems from the 
fact that the peaks on the mass scan were offscale to begin with. 

The ratio of the peaks may not be off scale but It cannot be 


determined. 
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Figure 94. - PSU mass scans on two Nl-doped samples, one as- received and the 
other after an exchange at 1330° C for 1920 hr. 
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As verification of the contamination, cations were profiled 
along with the tracer on some samples. A plot of the counting rate 
of silicon normalized to that of aluminum, for the sample exchanged 
at 1550° C Is shown In figure 95. The concentration of silicon near 
the surface Is ~8 times that at a depth of 0.5 ym Into the sample. 

If we assume the concentration at 0.5 ym Is the background, then, 
based on the spectrochemlcal analysis, we can estimate the near sur- 
face silicon concentration to be 350 ppm. 

Attempts to subtract out the background, assume a surface con- 
centration, and fit the profile with complementary error function, 
proved unsuccessful, (see fig. 96). This may be due to simultaneous 
contamination by Ca, Mg and K. 

The cations, SI, Na, Ca, and K were profiled on a sample ex- 
changed at 1300° C. The oxygen-18 profile Is shown In figure 9T, 
and the normalized cation profiles In 98 to 101. In this sample the 
cation profiles are much deeper than the tracer profile. Note that 
the profiles all contain a spike at the surface. Also, the SI, Ca, 
and K profiles have negative slope Indicating an Inward flux, while 
the Na profile has a positive slope Indicating an outward flux, or 
evaporation. 


Four samples; two undoped, one Cr-doped, and one Nl-doped, were 

exchanged In the 1634° C low P n experiment. The oxygen concentra- 

u 2 


tlon profiles In all four were Indistinguishable. The silicon pro- 
files In the undoped and Cr-doped samples were also very similar. 
However, the silicon profile In the Nl-doped sample Is distinctly 


3.00 



profile for silicon on an undoped sample exchanged at 
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profile. 
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Figure 97. - PSU SIMS profile of oxygen tracer In sample UO-LT, exchanged at 
1300° C for 1073 hr. 
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Figure 99. - PSU SIMS profile of sodium normalized to aluminum. In sampl 

lift I T /\ u « U . - 1 _ 1 « A AAA M . r 






Figure 101. - PSU SIMS profile of potassium normalized to aluminum. In sample 
UO-LT exchanged at 1300° C for 1073 hr. 
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dlfferent. This Is shown In figure 102. The undoped samples have 
just a spike of mass-28 near the surface, whereas, the penetration 
of silicon In the Nl-doped material Is roughly equal to that of the 
tracer. This dependence of contamination upon dopant Is Interesting 
partlcuarly since a dopant effect upon oxygen diffusion Is absent. 

It would be useful to have a reanalysis of other samples to see If 
this Is an Isolated case or If It Is systematic. 

The high purity Al^Og furnace tube visibly changed as a function 
of time at annealing temperatures. The as-received tubes were smooth 
and a waxy yellow color. After a long time at high temperature, the 
color became distinctly whiter and the surface roughened. 

To check on whether the,mater1e1 evaporating from the tubes 
could act as the source for the contaminants, a piece of a new and a 
piece of well used tube were scanned. The results are shown In fig- 
ure 103. The tubes do have all the contaminants present and their 
concentrations do drop off after the tube Is used. 

Contamination of the surface of single crystal Al^Oj sam P^ es » 
after annealing In commercial "high purity" alumina tubes has been 
observed by both Scott^ 69 ^ and Heuer^ 7 ^. The contamination was re- 
ported to consist of small droplets of a silicate glass containing 
Al, Ca. and K. (69) 

Examination of the surfaces of our samples revealed the same 
type of droplets. The droplets were tenacious, they were not dis- 
turbed during the removal of a sputter deposited Pd coating from the 
sample by physical rubbing with a Klmwlpe saturated with acetone. 
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They were, however, etched completely away by a 10 sec Irrmerslon In 
48 percent HF. 

The population and composition of the droplets varied from 
sample to sample, and from location to location on a single sample. 
Figure 104 shows typical droplet densities for two samples. The one 
exchanged at 1600° C has a much higher droplet density than does the 
sample at 1300° C, due presumably to the difference In vapor pres- 
sures of the glass formers at those two temperatures. Figure 105 
shows the droplet population, both near the edge and towards the 
center of a single specimen. The region of the sample near the edge 

was rounded during polishing and surface steps developed during 

annealing. In this region, the droplets were less numerous and of a 

larger diameter than those away from the edge. 

The radius of curvature may be estimated from the Newton's rings 
which appear In the large droplet shown In figure 106. The sample 
was Illuminated with an Hg source and each ring corresponds to ~l/4 
vm. The droplet Is around 2 ym high and — 1 20 ym In diameter. The 

_4 

radius of curvature Is therefore very small, ~5xl0 ym. This Indi- 
cates that the glass wets A1 2 0 3 very well. 

The chemistry of the droplets also appears to vary. Silicon 
was present In every case examined however the alkalis and alkaline 
earths differ. Increases In Mg, Ca, and K were observed by SIMS, In 
one sample (fig. 91) while Increases In only K and Na were observed 
In another sample, and further one sample +vas analyzed with EOAX, In 
an SEM, found only Ca and SI present In the droplets. It Is reason- 
able to expect the oxides evaporating from the furnace tube to change 
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Figure 104. - Optical micrographs of the typical glass droplet populations on 
samples annealed at a) 1300° C and b) 1600° C. (400X) 
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Figure 105. - Optical micrographs of the typical glass droplet populations 
a sample annealed at 1634° C a) away from the sample edge and b) near the 
sample edge. 
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composltlon as a function of time. Initially, higher vapor pressure 
species (with high alkali content) would be dominant and at later 
time the glass would be depleted of alkali. The area density of the 
droplets varied greatly. Regions of a millimeter square were ob- 
served to have no droplets while other areas contained on the order 
of 10* droplets per square centimeter. 

The existence of the droplets raised a question as to whether 
the cation profiles detected by SIMS are truly located within the 
A1 ? 0 3 or all generated by the sputtering of droplets during 
profiling. 

All of the cation profiles have a very sharp, and Intense count- 
ing rate during the Initial stage, and then exhibit a discontinuity 
followed by a smooth profile (see fig. 95 and figs. 98 to 101). 

This Is consistent with sputtering off of glass droplets and then 
profiling A1 ? 0 3 . The RBS data (see fig. 47) also suggests that 
cations In contact with the surface during the anneal will diffuse In. 

A sample which has been cut In half and the glass beads etched 
off of one half was analyzed. A cation profile was observed In both 
specimens. The magnitude In the etched sample roughly two thirds of 
the unetched. 

4.2.6 Contamination of Polycrystals 

Cation profiles similar to those In the single crystal samples 
were observed In the polycrystalllne material. An example Is shown 
In figure 107. 

Another possible method of contamination Is cross-contamination 
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~ Cation profiles In undoped hot pressed Alpth annealed at 
decade C F ° r clar,ty the profile was translated along the y-axls 
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between samples with different dopants. A set of hot pressed Al^O^ 
samples, singly doped with 500 ppm of Cr, N1 , Y, and Zr, were ex- 
changed together at 1300* C for eight hours. The oxygen-18 penetra- 
tion was too small to be measured. However, a profile for chromium 
was observed In the undoped specimen. This Is shown In figure 108. 
The procedure of simultaneous exchange to ensure that samples will 
see the same diffusion anneal may be unwise, particularly for such 
volatile dopants as chromium. 

It Is also observed that the samples with dopants that produce 
a color, Cr and N1, left the shelf In the annealing furnace colored 
after the anneal was complete. Although traces of the Y and Zr 
dopants were not visible. It Is quite probable that they have similar 
behavior. Therefore, the concentration of dopant In the near surface 
region Is lower than that of the bulk and changes with time. There- 
fore, the procedure for exchange of these types of samples will re- 
quire considerable refinement. 
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CHAPTER V 


DISCUSSION 


5.1 ANALYSIS TECHNIQUES 

Two quite different techniques were used In this study to deter- 
mine tracer concentration profiles, SIMS and proton activation. In 
our proton activation work there was some evidence of misleading data 
due to fluorescence of our samples. Earlier workers, at C.W.R.U., 
would not have seen this since they used a meter long magnetic 
spectrometer between the sample and the detector. If It turns out, 

In fact, that there Is a problem with photons striking the detector. 
It should not be a difficult problem to correct. 

There remain other advantages of SIMS over proton activation. 

(I) Depth Resolution - Even If the electronics and the sample 
were perfect, the best commercial detector gives energies +10 to 13 
KeV (13 KeV corresponds to 0.03 ym In Al^O-j), so the depth resolution 
would be limited. SIMS has routine depth resolution <10 A. 

(II) Penetration - The ranges of protons and alpha particles 
only allow a maximum penetration of ~3 ym. This Is particularly bad 
because of the relatively poor depth resolution. In SIMS the mate- 
rial can be sputtered to any depth of choice. However, concerns 
about machine time and taking data from the the region Into which 
the primary beam has been Implanted, also restricts practical use of 
SIMS to the first few microns. 
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(111) Experimental and Analytical Complexity - Proton activa- 
tion requires a Van de Graaf generator for protons. Also, the 
Interaction of protons on a crystal necessitates care to avoid 
channeling. The primary Ion beam In SINS may also be channeled and 
reduce the sputtering rate. However, since the sputtering rate was 
determined for each sample after the experiment, a channeled beam 
would not lead to errors In apparent depth. 

The conversion of raw data to concentration versus depth Infor- 
mation Involves many more variables with proton activation than with 
SIMS. Each variable carries with It an uncertainty. 

< 1 v) Species Detected - The basis for proton activation Is the 
very specific nature of the nuclear reaction. Only oxygen-18 Is 
detected. The addition of RBS allowed for some chemical Information 
however only species more massive than A1 could be resolved. 

SIMS allows many, up to nine, elements to be concurrently 
monitored. It was this feature which uncovered the contamination. 
Such a feature would also be of great value when studying the oxida- 
tion of metals. This would allow concentrations of cations to be 
monitored along with the tracer. 

(v) Destructive Nature - Proton activation Is essentially a 
nondestructive method. It Is a very small fraction of the tracer 
which reacts. SIMS, on the other hand. Is destructive. The analy- 
sis Is dependent on the sputtering of material. This limits the 
technique to a one time analysis. However, a sample 0.5 cm square 
will accommodate up to twenty craters without overlapping. In prac- 
tice, therefore, this Is not much of a limitation. 
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Overall, SIMS appears a more versatile tool for diffusion work. 


5.2 OXYGEN DIFFUSION 

The results of this Investigation may be summarized as follows: 

(1) The diffusion Is Insensitive to Impurity content. Our five 
boules, though not doped In the manner Intended, did have different 
Impurity levels between them. The near surface region was also 
demonstrated to have picked up Impurities during the anneal. 

Despite this, our data agrees very well with that of Reddy and 
Cooper. Since both studies used the same system with the same 
supplier of 99.8 percent Al^ tubes. It might be argued that both 
made diffusion measurements on equally dirty samples. However, the 
activation energy, 615 kJ/mol, agrees very well with that of the 
other gas exchange studies; 648 kJ/mol for Olshl and Klngery, and 
636 kJ/mol for Olshl, Ando, and Kubota. In these studies the sample 
was heated by Induction on a Pt-Rh susceptor surrounded by a water 
cooled SIO^ tube. Even If contamination occurred It would be 
unlikely that It would be the same as ours. The activation energy 
determined by lagerlof and Heuer In their dislocation loop annihila- 
tion study was also similar, ~600 kJ/mol . The furnace In which they 
annealed their TEM foil was porous firebrick. 

As mentioned In the literature survey, Olshl, Ando, and Kubota 
also found no effect on oxygen diffusion from Impurities. 

It Is Interesting to note that there Is an Indication that the 
rate of contamination of our samples by a cation Impurity, SI, was 
dependent on dopant. In the case of the Nl-doped crystal the appar- 



-197- 


ent penetration by SI was deeper than In the undoped material (see 
fig. 102). However, It should be kept In mind that this transport 
takes place on the cation sublattice. 

(11) The diffusion of oxygen In a-Al^Oj does not depend on the 
partial pressure of oxygen. The diffusion coefficient determined 
from the sample exchange In a CO/CO,, atmosphere, corresponding to 

P n » 10 -15 atm, falls on a line determined by the P n * 1 atm data. 
u 2 2 

Reddy suggested that the difference In P_ and that of Reed 

u 2 

might be responsible for the difference between the data. In light 
of the CO/CO^ experiment, this does not appear valid. 

5.2.1 Thermodynamics of Point Defect Clustering 
It has recently become clear that clusters of associated defects 
are very Important In oxides. The calculations of Catlow et al., 
which will be discussed more thoroughly. Indicate that the ratio of 
associated defects (clusters) to unassociated (free) defects. In 
o-Al^Og, Is very large. Before attempting to Interpret the results 
of self-diffusion experiments It Is necessary to consider the thermo- 
dynamics of clusters and the possibilities for their motion. 

Point defects (vacancies. Interstitials, and Impurities) may be 
treated using the thermodynamics of solutions. The host crystal Is 
treated as a solvent and the defects as solutes. In the following 
discussion of defects, and defect clusters, three assumptions are 
made 

(1) That the only Interactions between defects are those be- 
tween nearest neighbors, giving rise to clusters. 
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(II) The kinetics of the association of defects Into clusters 
Is rapid and therefore equilibrium Is maintained. 

(III) The clusters have a sufficiently long half life against 
breakup Into free point defects that the two may be treated as sepa- 
rate species. 

These assumptions are equivalent to having a dilute solution of 
atomic and molecular species. 

5 . 2 . 1 . 1 General Background 1 

Recall that the enthalpy of a system, H, Is related to the 
Interal energy, E , the pressure, p, and the volume, V, of the system 
by 

H * E ♦ pV (5.2-1) 

and the free energy, F, Is related to the Internal energy and the 
entropy, S. by 

F . E - ST (5.2-2) 

The Gibbs Function, G, then Is 

G ■ F ♦ pV (5.2-3) 

= E - ST ♦ pV 

- H - ST 

In the case of dilute solutions the volume, energy, and enthalpy 
will Increase linearly with small solute concentrations (at constant 
pressure and temperature) as 


1 


This background section Is paraphrased however directly from Flynn. 


(71) 
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(5.2-4) 


V * V ♦ nv 
o 

E « E ♦ ne 
o 

H - H ♦ nh 
o 

where V , U , and H are properties of n mols of pure solvent 
ooo o 

In which n mols of solute are dispersed, the v, e, and h are 
partial molar quantities. In turn, the variables S, f, and 6 
become 

S ■ S ♦ S ♦ ns (5.2-5) 

o c 

F = F - S T + nf 
o c 

G = G - S T ♦ ng 
o c 

where S £ Is the configurational entropy while s, f, and g are 
partial excess quantities. Note for the case of thermally generated 
free point defects, equation 2.1-11, 

AG = G o - (G o - S c T ♦ ng) = S c T - ng (5.2-6) 

From the above a g-functlon may be defined 

g » h - sT ♦ pv (5.2-7) 

The chemical potential, v, of a defect species, a, may be 
written as 



m~ < ScT > 


(5.2-8) 


For crystalline solids the configuration entropy may be derived 
from the formula 


In (*) 


(5.2-10) 
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where the degeneracy, *, Is the number of distinguishable configura- 
tions which the system may adopt. For the case of unassociated, or 
free, defects 


where c Is the concentration of defect species <*. The case of a 
a 

tightly bound cluster must take Into account the multiplicity of con- 
figurations, m , which Is the number of distinguishable differences 
In the shape and orientation of the cluster associated with each 
lattice site. For this case the configurational entropy may be 
written as 


-R n In 

q m q 


) 


(5.2-11) 


where the subscript * q ' refers to a cluster made up of a number, q, 
of defects. 

The chemical potential for equilibrium thermal defects 
(Intrinsic disorder) Is equal to zero, therefore for free defects 

v = 77 * ♦ RT In (c ) ♦ g =0 (5.2-12) 

o on tt o 

o 

or 


c . ■ * ,p 'FT 1 


In an analogus manner for defect cluster 


c q - M , **•> ‘rT 3 ) 


(5.2-13) 


(5.2-14) 
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When the cluster Is looked at as being formed by a reaction of 
free defects 

qX -» <qX) (5.2-15) 

The concentration of clusters may be written In terms of the free 
defect concentration as 

b 

c q * », c « RT (5.2-1M 

with 

b q = q • < g-j - g 2 ) 

where Is the energy to form a free defect and g ? Is the energy 

to form and associate q defects. 

5 . 2 . 1 . 2 Clustering In o-Al ^O^ 

Four types of neutral defect clusters, directly Involving the 
oxygen sublattice, are of Interest In Al^Oj. The first Is the 
association of a stoichiometric set of fully Ionized vacancies (a 
Schottky quintuplet) Into a neutral cluster (a Schottky Cluster). 

The rationale behind considering only the Schottky Cluster and not 
off stoichiometric clusters, Is to give a defect which Is Independent 
of both Impurity level and oxygen partial pressure*(as will be seen) 
and that association will take place to minimize electrostatic 
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2 

misfit. The other Imaginable neutral clusters of Intrinsic 
defects are Frenkel pairs. The existence of these are rejected 
because the Ion In such a bound pair must minimize Its energy by 
relaxing back to the lattice side. For this to be untrue the 
crystal structure would be unstable. 

The second type of cluster to be considered Is two substitu- 
tional acceptor Ions (divalent Ions on aluminum sites) with an 
oxygen vacancy. The third type Is two substitutional donor Ions 
(tetravalent Ions on aluminum sites) with an oxygen Interstitial. 

The fourth type of cluster Is that of a substitutional donor 
with a substitutional acceptor. 

Consider first, Schottky quintuplets and Schottky clusters. 
Schottky quintuplets are formed according to the reaction 

Null - 2 V ;; ♦ 3 Vq‘ (5.2-17) 


with a reaction constant which may be written 



■ ■ , 

where a^ Is the activity of and a^ Is the activity of V Q ’. 

Activity of a species may be written as 


a * yc 


2 

Charged clusters Involving Impurities have been proposed to exist In 

72 1 

for example MgO , and one such cluster Is Cr^ - 0 - , however 

9 9 

the author Is unaware of the proposition of an of f-stolchlometrlc 
cluster of Intrinsic defects In an Ionic material. 
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where c Is the species concentration and y Is the activity co- 
efficient. By assuming Ideal solutions, y = 1. concentrations may be 
substituted Into equation (5.2-18); 

.,2 3 -g 

k S = [V A1 ] lV 0 1 exp ( RT > (5.2-19) 


The association of these defects by the reaction 


2V A i OV-’, (2V a1 3V 0 ) 


yields an association constant, k a , which may be written 


K3V A1 3%n -9, 

k a ' .,2 3 ■ “s exp RT 1 

tv A1 ] (v 0 -) 


(5.2-20) 


where m $ Is the multlpllcty of configurations for Schottky Clusters. 
Solving for the cluster concentration yields 

l(2V A1 3V 0 ) X ] = m $ [V A ’j 2 r V q ' ] 3 exp (^) (5.2-21) 

This Is of the form of equation (5.2-16), with b = -g . It 

q a 

can also be seen that the Schottky constant may be substituted for 
the free defect concentrations In equation (5.2-20). This Implies 
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that, at constant pressure, the concentration of Schottky Clusters Is 
uniquely determined by the temperature of the system, and Is In- 
dependent of concentration of Impurities. 

-9 -9 

[<2V A1 3V 0 ) X J = k a k s = m $ exp (^> exp (^) (5.2-22) 

substituting for the g-functlons gives 

s s -h -h 

1(2V A1 3V 0 )X] = m s exp ( ~ ± T S ) exp ( ~ Rf ~ 1) (5.2-23) 

Secondly, consider Impurity defect clusters. 

For the case of a foreign donor Ion, Fd* + , compensated by oxygen 
Interstitials the Incorporation reaction may be written 

2 FdO ? l 2 Fd A1 ♦ oj ♦ 30 X ' (5.2-24) 

These defects may associate by the reaction 

2 Fd ftl ♦ 0*' l (2 Fd A (jl y )* (5.2-25) 

At this point two variables, a’ and a may be defined, a' Is the 
degree of association (the reaction proceeding to the right) and a 
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1s the degree of dissociation (the reaction proceeding to the left). 
They are related by 

o* = 1 - a 

Under the simplified electroneutrallty condition of 

[Fd] * 3[oJj (5.2-26) 

(which Is only rigorously true when there are no other charged 
species present, but becomes a good approximation when background 

3 

charged species are of low relative concentration ) the concentra- 
tions of the species may be expressed as the product of the dopant 
level, [FdO^], and these variables; 

[(2FdV Q ) X ] = I 1 [FdO ? ] = 1 [Fd0 2 ] (5.2-27) 


[Fd ftl ] = (1 - a 1 ) [Fd0 2 ] * a [Fd0 2 ] 

[ojj - i^ 2 ] = 3 [Fd0 2 3 


The association reaction constant may be expressed in terms of 
the degree of association as 


a'[FdO ) -g 

* 3 3 = m (j ex P ( ot ) 

a (1 - a' )^ [Fd0 2 ) J ° Kl 


^The question of when this Is a valid approximation Is dependent on 
the relative concentrations and also the power of the expression 
In which It Is used. For example, If the donor under considera- 
tion makes up 90 percent of all those present, the 10 percent 
error of this assumption Is magnified to 27 percent when this term 
Is cubed. 
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The degree of association may be directly obtained for any 
values of k and [FdO-J from; 

a c 

[FdO,] 2 (5.2-29) 

(1 - a'r 

and the concentrations of the various species determined. In the 
limiting cases of nearly no association and nearly total association 
defect concentration may be directly obtained. In the limiting case 
of the degree of association Is small 

(1 - o') 3 r 1 (5.2-30) 

and solving for the concentration of clusters obtains 

[(2 FdV 0 ) x ] = 3 [FdO ? ] 3 m d exp (^) (5.2-31) 

It can be seen that there Is a very strong dependence, third power, 
on the dopant concentration. 

For the other limiting case It Is more convenient to use the 

degree of dissociation. The reaction constant then becomes 

( 

(1 - o) [FdO ] -g 

ka = 3 rFrfn ,3 = m d 6XP ( R1 _) (5.2-32) 

a l r aUp J 


In this case the degree of association Is very large, therefore the 
degree of dissociation Is small and 
(i-«):i 

solving for the concentration of free oxygen Interstitials obtains 
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In the case of a large association constant the concentration 
of unassociated points defects depend on dopant concentration only 
to the 1/3 power, and the temperature dependence of the exponent has 
been redued by a factor of three, as compared to the unassociated 
case. 

For a foreign acceptor. Fa, compensated by oxygen vacancies, the 
Incorporation reaction becomes 

2 FaO * 2 Fa^-j ♦ V" ♦ 20j (5.2-34) 

The association reaction becomes 

2 Fa Al * V ’ < 2Fa Al V* <5 - 2 - 35) 

with the electroneutrallty condition 

[Fa Al ] = 3 [V o'l (5.2-36) 

A set of equation analogous to these derived for the donor case may 
be obtained, for the case of very weak association the cluster con- 
centration Is given by 

-g 

[(2 FaV 0 ) x ] = 3[FaO] 3 m fl exp ( (5.2-37) 

and for the case of very strong association the vacancy concentration 
Is given by 

[Vq*] - <[FaO]/m a ) 1/3 exp( ^ ) (5.2-38) 

Finally, consider donor acceptor clustering. In the above the 
assumption was made that the crystal was singly doped with either a 
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donor or acceptor. It Is quite possible for both to be present. In 
this case the free oxygen vacancies and Interstitials which are pro- 
duced are related through the Anion Frenkel product. 

°i : * v < 5? - 39 > 

k F - CO*] [Vj ] • exp (5.2-40) 

The Impurity-defect clusters can react to form donor-acceptor 
clusters according to 

(2 FaVg) x ♦ (2 FdO^) X . 2 (FaFd) X ♦ V" ♦ oj (5.2-41) 

with a reaction constant 



[(faFd) x ] ? [Vq’HO^ 

[(2 faV Q ) x ] [(2 FdO^} = da 


exp ( 



(5.2-42) 


The g da Includes the dissociation energies of the impurity- 
defect clusters as well as the association energy of the donor- 
acceptor cluster. The free defect concentrations may be replaced by 
the Frenkel constant reducing the energy penalty for reaction 


[<2FaV 0 ) x ] [(2Fd0 1 )) 


_(9 da ' g F ) 

* m da exp( — ) (5.2-43) 


In order to predict the outcome of doping the g-functlons of all 
of the above reactions, and the multiplicities of configurations for 
the clusters, must be known. If the g-functlons are very small, so 
that k Is very small, clusters can be minority species and the 

d 

point defect equilibria may be dealt with using the notions of 
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Schottky and Frenkel disorder. This Is a special case, however, and 
recent calculations along with some experimental evidence (which are 
discussed In a following section) Indicate that the g-functlons of 
at least some clusters are substantial, to the extent that clusters 
may be significant, perhaps dominant, contributors to diffusion. 

Prior to carrying out calculations of defect equilibria It Is 
of benefit to assess the sign and magnitude of Individual terms. 
Consider first the geometry of clusters. The structure of a-Al^Og 
consists of a very nearly hexagonal close packed oxygen sublattice 
with two thirds of the octahedral sites occupied with an aluminum 
Ion. Each aluminum has sixfold coordination as Illustrated In fig- 
ure 109. The a-Al^Og structure may be thought of as being made up of 
occupied and unoccupied octahedra. The packing sequence, which maxi- 
mizes the average Al-Al distance. Is shown In figure 110 and a 
stereographic pair of the occupied octahedra Is shown In figure 111. 
It can be seen that each occupied octahedron shares a face with one 
other occupied octahedron, shares an edge with three other occupied 
octahedron, and shares a corner with nine other occupied octahedra. 

In figure 110, E shares a face with A, an edge with F,* G, and H, and 
a corner with B, C, 0, I, J, K, L, M, and N. 

The true structure of Al^Og differs from this Idealized struc- 
72 

ture In two ways, both of which arise from the fact that only two 
thirds of the octahedra are occupied. First, the octahedra are not 
regular. The edges of the faces shared between occupied octahedra 
are shorter than those of the unshared faces. This Is due to the 
component of the Al-0 bond along the edges being twice as great In 



coordinated aluminum 1on.° a ° oxygen octahedron containing a sixfold 














Figure HO. - Schematic of the pack 
correspond to the octahedral site 
those occupied by aluminum Ions. 



Right eye Image 
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Figure 111. - Stereo pair of the packing of occupied octahedra In the Idealized 
A 1 j >03 structure. Place binding away from viewer when viewing with 
crossed eyes and place binding toward viewer when using a stereographic 
viewer. 






octahedron thick (drawn with u = -0.083, w/^ = 0.017). Small black 
circles are A1 In plane x ♦ y = 1; heavy and light open circles are 0 lying 
In front of and behind this plane, as viewed from the lower right-hand 
corner of Fig. ll-5(b). Corners of the unit cell are Indicated by broken 
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72 

the shared face. Secondly, the cation-cation repulsion causes 
the aluminums to be displaced from the octahedra centers, toward the 
unshared face, figure 112. Note that these two relaxations of the 
structure have opposite effects on the Al-0 bond distance. The dis- 
placement of the aluminums has the slightly greater effect and the 
Al-0 bond Is 2.5 percent shorter In the unshared face (1.89 A) than 
In the shared face (1.95 A). Bond strength Is related to bond dis- 
tance, however based on the fact that due to unknowns, this dis- 
cussion will be forced to neglect larger dimension variations, and 
relaxations of defect clusters will not be taken Into account; the 
Idealized structure will be employed. 

If the only constraint on cluster configurations Is that It con- 
sists of an unbroken string of nearest neighbors then the Schottky 
Cluster could Involve as many as four octahedra, see figure 113. 
However It seems reasonable, though arbitrary, to employ the follow- 
ing two rules for vacancy type clusters: 

Rule 1.: The clusters Involve two octahedra. It Is not possi- 

ble to use less than two since there Is only one aluminum site per 
octahedra, however all clusters may be obtained using only two. 

Beyond avoiding the Inconvenience of dealing with the large number 
of configurations possible when three or more octahedra are Involved, 
this effectively minimizes the cluster equivalent diameter, and hence 

maximizes the local defect density, which Is the driving force for 

75 

the association process. Oamask and Dienes make a similar 
argument for maximizing nearest neighbors during the formation of 
01- and Trl- vacancies In metal. They Justify this by noting that 
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the number of broken bonds Is diminished when vacancies become 
nearest neighbors. 

Rule 2.: The coordination of the defects, by oppositely charged 

defects be maximized. In order to achieve local neutrality. In 
vacancy clusters this effectively requires that the shared oxygen 
sites be occupied at the expense of the nonshared. 

Using these rules and neglecting relaxation there are three 

types of Schottky clusters. The first Is between face sharing 

octahedra, figure 114, In which the cation-cation site distance Is a 

minimum (0.8165 times the anion-anion site distance) and the defect 

coordination numbers are maximized (each V., has three nn V_ and each 

A I U 

V_ has two nn V.,). The multiplicity of configurations In this case 
u Al 

4 

(distinct orientations over occupied sites ) Is 1/5 since the 
location of all five vacancies Is fixed by Rule 2. The second type 
of cluster Is between two edge sharing octahedra, figure 115. In 
this the cation-cation distance Is longer (equal to the anion-anion 
distance) and the coordination numbers are lower (average V ftl has 
2.5 nn V Q and average V Q has 5/3 nn V A1 ). The multiplicity of 
configuration Is 8/5 since the nonshared oxygen vacancy may be at 
any of the eight nonoccupled sites without changing the average 
coordination numbers. The third type, between corner sharing 


*The multiplicity of configurations given here Is per octahedra pair, 
the total multiplicity Is this number multiplied be the number of 
distinct octahedra pairs of each type which will be termed the 
octahdra coordination number. Referlng to figure 110, It can be 
seen that this Is equal to one for the face sharing, three for the 
edge sharing and nine for the corner sharing configurations. 
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figure 114. - Illustration of a Schottky Cluster Involving two face 
sharing octahedra. The location of all the vacancies Is fixed. 
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octahedra. One oxygen vacancy has an unspecified location and may be 
located on any of the remaining eight unoccupied oxygen sites. 








vacanc 
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gure 116. - Illustration of a Schottky Cluster Involving two corner sharing 
octahedra. Two oxygen vacancies have unspecified locations and may be 
located on any of the ten remaining unoccupied oxygen sites. 
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octahedra (fig. 116), has a still larger cation-cation 
d1stance(1 .3478 times the anion-anion distance) and still lower 
coordination numbers (average V A1 has 2 nn Vg and average Vg 
has 2/3 nn V Q ) . The multiplicity of configurations Is 90/5 since 
the nonshared oxygen vacancies have ten equivalent sites to choose 
from. The differences In Interdefect spacing and coordination 
numbers Indicates that the energy of the cluster will vary with 
type. 

The decrease In coordination of the corner sharing configuration 
over the face, or edge, sharing which suggests higher energy may 
tend be offset by the Increase In the distance between the cation 
vacancies . 

Since It Is not a simple calculation to obtain the net differ- 
ences, the clusters will be assumed to have the same probability of 
existing. In order to calculate a multiplicity of configurations. 

When this Is true there are 835 Shottky clusters which will Involve 
any given lattice site which gives a total multiplicity of configura- 
tions equal to 47. This enters Into the preexponentlal of equation 
A1 as a positive term tending to Increase the degree of association 
Into Schottky Clusters. 

For the case of a foreign acceptor cluster, Fa-Vg-Fa, there are 
also three types. The cation-cation distances are the same as the 
Schottky Clusters and In all cases the defect coordination numbers 
are the same (each Fd has 1/2 nn V Q and each V Q has 2 nn Fd). For 
this cluster the multiplicity of configurations Is equal to the num- 
ber of shared vacancy sites. For face sharing (fig. 117) m = 3/3, 

O 
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Figure 118. - Illustration of a cluster of two substitutional Acceptor 
Impurities and an oxygen vacancy, Involving two edge sharing octahedra. 
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for edge sharing (fig. 118) m « 2/3, and for corner sharing (fig. 

119) m ■ 1/3. The total number of Fa-V n -Fa clusters Involving 
aglven lattice site Is 15 which gives a total multiplicity of 5. This 
will enter the preexponentlal of equation 5.2-37 favoring association 
elation of Fa-Vg-Fa clusters. The case of foreign donor dusters, 
Fd-O^-Fd, Is more Involved, and some additional freedom will be 
allowed. If a dl-Interstltlal forms and occupies an oxygen lattice 

N 

site, 0, , the geometry Is perfectly analogous to that of the 
acceptor cluster. If the oxygen Is Interstitial In the classic 
sense the geometry Is dependent on the allowed Interstitial sites. 
Three types of oxygen Interstitial sites will be considered, a 
triangle, a tetrahedron, and an octahedron. In the case of two face 
sharing octahedra there Is a triangular site, In the shared face, 
that Is the only nearest neighbor cluster Involving an oxygen 
Interstitial, (by nearest neighbor It Is meant that the Fd-0^ bond 
distance Is less than, or equal to, the Al-0 bond distance In 
oAl^Og). This Is shown In figure 120. The cation-cation 
distance Is 0.8165 and the defect coordination numbers are the same 
as the Fd-Vg-Fd clusters. The size misfit for Interstitial Is large 
(Ideal radius ratio for triangular site Is > 0.155, and Is In this 
case « 1), while the donor-interstitial distance Is minimized, 

0.4032. It has been pointed out by Gourdln and Klngery 71 that 
radius ratios In Ionic materials may be misleading. 

There are six tetragonal sites around each octahedral site. If 
these are availed to the Interstitial a number of possibilities occur. 
In the case of face sharing octahedra the oxygen Interstitial may 
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equlvalently occupy any one of the twelve total tetrahedra, figure 
121. In this face sharing conf Iguratlon, however, the distance 
between the positive and negative defect Increases (from 0.4082 to 
0.5774) and the Interstitial Is a nearest neighbor to only one of 
thedonors, the coordination numbers are cut In half. This suggests 
a higher energy cluster though to what degree Is unknown. 

In the octahedra edge sharing configuration, with tetrahedral 
occupancy, figure 122, the coordination numbers are restored (each 
Fd A1 has 1 nn 0^ and each 0^ has 2 nn Fd^). The cation-cation dis- 
tance Is Increased and the radial misfit Is smaller (Ideal radius 
ratio for tetrahedron Is > 0.225), compared to the face sharing, 
triangular occupancy, configuration. However, the Fd-0^ distance, 
0.5774 Is longer. The Interstitial may occupy either of two equlva 
lent tetrahedra. 

In the octahedral corner sharing configuration, with tetrahedral 
occupancy (fig. 123), the oxygen Interstitial, like face sharing, can 
only be a nearest neighbor to one of the donors. To retain a next 
nearest neighbor relation with the other donor there are four equiva- 
lent tetrahedra. 

Since only two thirds of the octahedral sites are occupied with 
aluminum Ions, one third remains available to the oxygen Intersti- 
tials. This cluster has the least radial misfit (Ideal radius ratio 
> 0.414 for octahedral site) but Involves three octahedra and the 
distance between oppositely charged defects Is maximized. If only 
linear clusters are allowed (In order to maximize the Fd-Fd distance) 
then, referring back to figure 110, each lattice site Involves one 
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Interstitial, Involving three face sharing octahedra. 










Normally unoccupied 
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gure 126. - Illustration of a cluster of two donor Impurities and an oxyg 
Interstitial, Involving three corner sharing octahedra with the normally 
unoccupied octahedral site filled. 
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face sharing configuration along the c-axis (fig. 124), six edge 
sharing configurations In the basal plane (fig. 125), and three 
corner sharing configurations along < 1011 > (fig. 126). 

Thus, with three types of Interstitial sites allowed, there are 
68 clusters which will occupy any given A1 site giving a multiplicity 
of configurations of 34. If the highly asymmetric clusters, with 
tetraganal Interstitials, are unstable the multiplicity drops to 10. 
As with the other two clusters this will enter the pre-exponential 
to favor clustering. 

Donor-acceptor clusters have thirteen total configurations, the 
number of octahedra In contact, and therefore a multiplicity of con- 
figurations equal to 1/2 for each octahedra pair 

These results are summarized In table 9. 

Consider now, the change In entropy on clustering. The partial 
excess entropy term differs for the two types of clusters. In the 
Schottky cluster case It contains a formation entropy and an associa- 
tion entropy term while the T1-0^-T1 cluster contains only an 
association entropy term. In both cases, however, any change In 
entropy arising from the association process will result In a change 
In cluster concentration. 

The partial excess energy term arises from local electronic, 
vibrational, and elastic changes arising the defects. It Is assumed 
that all the electrons are In the ground state and the electronic 
contributions will be Ignored. Flynn points out that significant 
electronic changes occur In materials such as transition metals. 

The vibrational contribution may be written as 
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TA8l£ 9. - UNRELATED CLUSTER PARAMETERS 


Cluster 

Type 

Octahedra* 

Multiol icity** 

Cation- 

Ave 

Ave 

Distance 


( shared 

coordination 

of configura- 

cation 

V A1 

v 0 

between 


feature) 

number 

tion oer 

distance 

coord 

coord 

nn ♦ ve 

Schottky 

Face 

1 

octaneora 

1/5 

0.8165 

3 

2 

defect and 
-ve defect 

0.7071 


Edge 

3 

8/5 

1.000 

2.5 

5/5 

0.7071 


Corner 

9 

90/5 

1.3478 

2 

2/3 

0.7071 

FaVjgFa 

Face 

1 

1/3 

0.8165 

Ave 

^At 

Coord 

Ave 

01 

Coord 

2 

0.7071 

Edge 

3 

2/3 

1.000 

1 

2 

0.7071 


Corner 

9 

1/3 

1.3478 

1 

2 

0.7071 

FdO-jFd 

Face 

(Tri) 

1 

1/2 

0.8165 

Ave 

Fd Al 

Coord 

1 

Ave 

v 0 

Coord 

2 

0.4082 


2-0ct 

face 

(Tet) 

1 

12/1 

0.8165 

1/2*** 

1*** 

0.5774 


Edae 

(Tet) 

3 

2/2 

1.0000 

1 

2 

0.5774 


Corner 

(Tel) 

9 

A/2 

1.3473 

1/2 

1 

0.5774 


3-0ct 

Face 

(Oct) 

1 

1/2‘t 

1.6330 

1*.. 

?*** 

0.8165 


Eoqe 

(Oct) 

3 

1/2 

2.0000 

1 

2 

1.0000 


Corner 

(Oct) 

9 

1/2 

2.6556 

1 

Z 

1.4141 

FdFa 

Face 

1 

1/2 

0.8165 

Ave 

Fd Al 

Coord 

1 

Ave 

Fa Al 

Coord 

1 

0.8165 


Edge 

3 

1/2 

1.0000 

1 

1 

1.0000 


Corner 

9 

1/2 

1.3478 

1 

1 

1.3476 


•Octahedra coordination number is equal to the number of distinct octahedra pairs 
associated with any given lattice point. 

••Multiplicity of configurations is equal to the number of distinguishable configu- 
rations, within the octahedra, divideo by the number of occupied sites. The 
total multiplicity of configurations per cluster type is equal to the multi- 
plicity of configurations per octaheora multiplieo by the octaheora coordination 
numoer. 

•••Note that there are oxygen sites between the cations and interstitials in these 
clusters. 

♦Normal itec to the length of an edge on an octahedron, which is equal to the 
oxygen-oxygen distance in the idealized structure. 

♦tAllowing only linear configurations. 



where u refers to the phonon frequencies for the perfect crystal 

o 

and «' refers to the phonon frequencies after the Introduction 

Ol 

of defects. Strlpp and Howard (quoted In ref. 71) find for an fCC 
lattice a value of S y * 1.74k fa for a vacancy however they also 
found an association S = -0.24k fa when a dl vacancy forms. Mahanty, 
Haradudln, and Weiss (quoted In ref. 71) found that divacancy forma- 
tion In a simple cubic lattice also has a negative S y . 

The elastic contributions consist of a strain term, S s , 

S * - 8wrV f* (5.2-45) 

S o I 

where 

n Is the strain 

v Is the shear modulus 

r Is the distance at which nr 

Is established 

i 

and a dilation term, 

d 

o r^nav (5.2-46) 

a Is thermal expansion coefficient. 

These terms are both positive for vacancy formation In alkali 

halides 71 because the neighboring atoms relax outwards. The effect 

on these terms of associating the free defects Into clusters Is not 

74 

known. However, the results of Pederson , on an FCC lennard- Jones 
system shows that the formation entropy (which Is the sum of all 
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three terms) of two free vacancies, 3.44k^, Is greater than that of 
a divacancy, 3 . 1 8k^ and the formation entropy of three vacancies, 
5.16, Is greater than any of the trivacancy configurations examined, 
4.65 to 4.73. 

In summary It appears that the preexponentlals, In the defect 
concentration equations, contain terms of opposite sign. Defect 
association to clusters Increases the multiplicity while decreasing 
the entropy. The net change for Schottky clusters may be positive 
or negative when the values for Al^O^ are employed. Indications are 
that it will be positive, but In either case It appears to be small. 
The product for the Impurity defect clusters will, also, vary from 1 
by only a small amount. 

Due to the uncertainty of their value, and Its temperature In- 
dependence, the preexponentlals, for all reactions, will be assigned 
unity. This Is felt to give a lower limit on cluster concentration. 

The partial excess enthalpy, at constant pressure, Is equal to 
the change In Internal energy of the system provided the system plus 
a pV term. In condensed systems the pV term Is small and may be 
neglected . 

With the above approximations, and estimations, the reaction 
constants for defect reactions may be calculated only taking Into 
account the change In Internal energy. For example, the formation 
of Schottky defects has a reaction constant 

k s = [VJj] 3 [V 5 ] 3 = exp <j^> = exp (j^> 


(5.2-47) 
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5. 2. 1.3 Cluster Mobilities 

The mobility of defects, which determines the extent of self- 
diffusion, Is dependent on their Jump frequencies. An Impurity 
boundto a defect may have a different jump frequency, with respect 
to the defect, than do the host lattice Ions. In addition, the Jump 
fre- quencles of host Ions themselves, may be altered by the 
presence of the cluster. Clusters are often thought of as 

Inherently less mobile than Isolated point defects. Damask and 
75 

Dienes point out that this Is not necessarily so. For example, 
divacancy clusters are more mobile than Isolated vacancies. In 
metals. They also postulate that an Impurity bound to a vacancy may 
cause the defect to Interact more often with the host Ions as well 
as the Impurity. Also, In a review of fast diffusion of Impurities 
In metal (when Impurity diffusion Is much faster than 
self-diffusion), by Warburton and Turnbull^, the same mechanism, 
enhanced Jumping of vacancy due to a bound Impurity, Is cited and a 
number of cases attributed to this are referenced. 

It Is not clear therefore how Jumping frequencies will be 
altered In the defect clusters In Al^Og. They may be higher or lower 
than In pure Al^Og. 

The correlation coefficient, f, may also be expected to be 
different for clusters and Isolated defects. Complexity reduces the 
correlation coefficient since a larger number of coordinated jumps 
must be made to move a unit distance and hence there are more oppor- 
tunities for backward Jumps (Jump opposing the direction of the 
cluster motion) to be made. For example, the value of . f for self- 
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dlffuslon by a vacancy mechanism In FCC metals Is 0.78146 77 while 

78 

for a divacancy f Is 0.475 , or slightly more than half. In 
Ionic solids the situation Is more complex owing to the differences 
In Jump frequencies of the constituents of the cluster, 77 which as 
Just discussed are unknown. 

The vacancy clusters may move within the following constraints: 

1. No dissociation of clusters. Each Jump Is to a stable 
configuration. 

2. No cooperative Jumps. All defects are allowed to move 
Independently. 

Within these rules the clusters may accomplish mixed direction 
motion, that Is, motion which has a component both parallel and per- 
pendicular to the c-axis. 

Illustration of possible paths for the Schottky Cluster, and the 
Fa-Vg-Fa cluster, are shown In figures 127 and 128. Only forward 
Jumps, those which aid In cluster translation are shown. It can be 
seen that In order to get from one set of face sharing octahedra to 
another requires 11 atomic steps In the case of the Schottky Cluster 
and only 7 steps are necessary for the F -V n -F . (Note, If a more 

dud 

favorable starting location has been chosen, the Fa-V^-Fa cluster 
could have made It In 6). In addition to being long sequences, each 
Jump requires the selection of one direction from a number of possi- 
bilities. For example. In the case of the Schottky Cluster, once 
the Jump from 1 to 2 has been made, In any of three equivalent direc- 
tions, the next Jump depends on the probability of which vacancy 
Jumps next and where It Jumps to. For the step 2 to 3 to take place 






nearest neighbor relations 
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requlres the vacancy to jump (which depends on relative jump fre- 
quencies for V.. and V.) and that It jump In one out of the 

A I 0 

three directions possible under rule 1. It Is obvious that an eleven 
step process will Involve vastly more than eleven jumps. 

In the case of donor Interstitial clusters the ability for 
motion, through point jumping, depends on the available sites for 
the Interstitial. If the asymmetric clusters are not allowed there 
would be no stable corner sharing cluster and the motion described 
above would not be possible. If these asymmetric clusters are 
allowed then a possible jumping sequence (7 steps) Vs shown In 
figure 129, for the donor clusters Involving only two octahedra. 

The three octahedra clusters can only move If a partial dissociation 
Is allowed, since the various stable form cannot be obtained through 
single jumps. An Illustration of a possible motion of this type of 
cluster Is shown In figure 130. 

Based on the number of steps In the point jumping schemes, the 
Impurity clusters may be expected to be more mobile than the Schottky 
Clusters. However, the energy of each jump has to be taken Into 
account. They are, unfortunately, unknown. The Initial Jump for the 
two octahedra donor clusters, figure 128, seems to be. Intuitively, 
a high energy Jump. Likewise, the partial dissociation of the three 
octahedra donor clusters (fig. 124) would be expected to have a 
significant energy penalty. While the potential of being misled by 
Intuition Is recognized. It Is possible that barriers of sufficient 
energy exist that these clusters are effectively trapped In their 
original configurations. 


242 


tA 

o» 

u 



3 

c r 

*o 


O JZ 
o 
c — 
o £ 
— 3 


o - 
E «— 
<9 

o> «- 

£ -M 



t- t/» 

O v. 

<*- o> 

E c 

V/» «r“ 

c c 
(0 0 > 
.xr o> 
o >> 

IS 

a> c 

r— * (Q 

jz 

f X 3 
<A C 
U> flj 

O 

Cl i/> 
<u 

(Q *r i 


o &- 

C Q. 

o E 


«o fc. 
i- o 
+■> c 

V* O 
3 T 3 

O 

>-* i 


<r> 

<v 


O t/* 
V- 3 
3 r~ 

a> W 


the Interstitial to go to a tetrahedral site. 



243 - 




o» 
•0 c 

< 4 - > 

O r- 
O 

c > 
o c 



o * 
E <- 
<0 
o> ■*- 



i- i/i 
O k- 

< 4 . 0 ) 

E c 

I/I •*- 

c tr 

<0 4 ) 
jC D> 
o >» 
ai x 
E o 


<v c 



*- X3 
i /> C. 
i/> ra 
o 

a. «/« 

a> 


«o *- 



o l- 
3 

e a. 
o E 



face sharing octahedra. 


-244- 


Obvlously, the existence of low energy cooperative jumps would 
greatly reduce the correlation problem. It Is also necessary to 
have cooperative jumps for any cluster to move directly from one set 
of a given type of octahedra to another set of the same type, e.g., 
edge sharing to edge sharing. 

Another aspect which cannot be Intuitively predicted Is the role 
of lattice relaxation around the cluster. The difficulty In pre- 
dicting relaxations Is Illustrated In the work of Vineyard and 
79 

Gibson on a tetravacancy In copper (FCC). Figure 131 shows 
their findings, a tetravacancy In a tetrahedral configuration (with 
nn maximized) was unstable, l.e., had a negative binding energy. 
However, If the vacancies were shared between six sites, and two of 
the face centered atoms grossly relaxed Inward, the configuration 
was stable, stability was Increased when the sharing was among eight 
sites. 

It Is clear that a great deal of the Information necessary to 
predict the atomic detail of clustering, and cluster motion Is 
currently unavailable. 

5.2.2 Impurity Insensitivity of Oxygen Self-Diffusion 

For diffusion to be Independent of the presence of allovalent 
Impurities, either the mechanism must Involve a defect whose con- 
centration Is Itself not sensitive to the Impurities, or the effect 
of the Impurities must be buffered out. 

The degree of the latter effect can be Illustrated using the 
calculations of Catlow et al. It may be recalled that In the case 
of Tl-doplng (Oonor) the choice of potential had such an Influence 
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that the selection of a dominant mode of compensation was not 
possible. For the sake of clarity only one example will be 
examined; the case of oxygen Interstitials produced by donor doping 
will be selected. 5 

Using the energies of Catlow, et al., the cluster dissociation 
constants were calculated using the expression 


k d - <rt> 


(5.2-48) 


and the results are shown In table 10. It can be seen that up to 
temperatures near the melting point, the dissociation constants are 
small and the defect equilibria In doped Al^O^ cannot be accurately 
predicted without taking clustering Into account. 

The binding energy of the donor cluster Is expected to depend on 
the particular donor Involved. For example, Gourdln and Klngery cal- 
culated defect energies In NgO and found that for the complex 

(2Fdg<j0j) x was ~ 127 If the donor was Al 3 * and -166 kJ (30 percent 

3t 

more) If the donor was Fe For the purpose of the following dis- 
cussion of donor doping In- Al^ consider the following three 
species, all compensated by oxygen Interstitials: 


5 The calculated binding energies for clusters (503 kJ for (3T1^iV^)) 
and 293 kJ for (2T1 aiO^)) tends to favor domination by 
aluminum vacancy compensation. However, only one, arbitrarily 
chosen, configuration was Investigated In each case. For vacancy 
mode this was three titaniums (on sites F, G, and H In fig. 110) 
surrounding one Al vacancy (site £). For Interstitials this was 
analogous to that In figure 124. It Is certainly possible that 
the favoritism could be reversed when an exhaustive study Is made, 
thus a freedom of choice Is allowed. It should be noted that It 
Is not the binding energy of any particular cluster which lies 
beneath this argument but the consistently high magnitude of all 
that were calculated. 
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TABLE 10. - THE DISSOCIATION CONSTANT, K d , FOR 
DIFFERENT CLUSTERS AS A FUNCTION OF TEMPERATURE 



T°C 

< 2T1 A1°1> <3T1 a 1 V a) > 

< 2M9 AlV (2M9 A1 H »1> 


500 

1 .46854E-20 

9.79684E-35 

2.25522E-17 

1.57956E-16 

600 

2.7461 2E-18 

7.7041 3E-31 

1.81987E-15 

1 . 01 989E- 1 4 

700 

1 .7521 8E— 1 6 

9.58546E-28 

5.95578E-14 

2.79591 E-13 

800 

5 . 1 5258E-1 5 

3.15963E-25 

1 . 01 735E-1 2 

4 . 1 3491 E— 1 2 

900 

8.51326E-14 

3.87551E-23 

1.07115E-11 

3.86304E-1 1 

1000 

9.05309E-13 

2. 23290E-21 

7.79124E-11 

2.54053E-10 

1100 

6.82248E-12 

7 . 1 2782E-20 

4.24457E-10 

1.26988E-09 

1200 

3.90833E-11 

1 .421 75E-18 

1.83695E-09 

5.10173E-09 

1300 

1.79332E-10 

1 . 93829E— 1 7 

6. 59877E-09 

1 . 71744E-08 

1400 

6.85842E-10 

1 . 93365E-1 6 

2.03440E-08 

5.00063E-08 

1500 

2.25463E-09 

1.48816E-15 

5.52400E-08 

1 .29066E-07 

1600 

6.52635E-09 

9 . 21 058E— 1 5 

1 .34818E-07 

3.01042E-07 

1700 

1 .69668E-08 

4.73889E-14 

3.00580E-07 

6.44403E-07 

1800 

4.02197E-08 

2.08177E-13 

6.20267E-07 

1.28174E-06 

1900 

8.80599E-08 

7/98048E-13 

1.19740E-06 

2. 39306E-06 

2000 

1.79958E-07 

2.71817E-12 

2.181 55E-06 

4.22913E-06 

(1) T1 

4+ 

with E b 

rounded to -300 kJ 


(2) A hypothetical 

donor, H2 4+ , 

with E b * -; 

200 kJ (30 percent 


less than T1) 

(3) A hypothetical donor, H4* + , with E b = -400 kJ 6 (30 percent 
more than T1) 

Also as an example, a temperature of 1600° C typical of a 
self-diffusion experiment Is chosen for the system. With the above, 
the following comparisons may be made. 


^Note : -400 kJ Is a reasonable value for a cluster binding energy 
It Is 100 kJ less than that calculated by Catlow, et al. for the 
cluster (3 T1/y] V^i) . 
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- Intrinsic vs. Doping 

It Is of Interest to compare the number of free Interstitials 
produced by doping with the number of thermally generated free Inter- 
stitials at our example temperature, 1600° C. 

The concentration of thermally generated defects may also be 
estimated using the calculations of Catlow, et al. Recall that 
Frenkel disorder was predicted to be dominant when the empirical 
potential was employed while Shottky disorder was predicted with the 
non-emplrlcal potential. 

In the case of Frenkel domination the electroneutrallty condi- 
tion Is 

[V'-] = [oj] (5.2-49) 


and an estimation of the Interstitial population may be obtained by 
simply substituting this Into the reaction constant, which Is 


k AF = fV 0 3 [0 1^ = exp ( RT * 


(5.2-50) 


giving 

N _e ^ 

[Oj] = exp (YJf) (2.5-51) 

In the case of Schottky disorder dominating the electron- 
eutrallty condition Is 


7 Note the defect energies given table 1 are per defect, for 
example, E*p = 2 defects/Frenkel pair = 366 kJ/mole of 
defects = -732 kJ/ mole of Frenkel pair. 
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2 l V ] " 3 [V A1 3 


(5.2-52) 


with this the concentration of oxygen vacancies may be obtained from 
the Schottky reaction constant 


f»o'] 3 tv“’l 2 « exp < 5 ^> 


RT 


(5.2-53) 


giving 


[v;*] 


<i> 


2/5 


exp ^ 5RT 3 


(5.2-54) 


substituting this Into equation (5.2-55) and solving for the concen- 
tration of Interstitials obtains 


[0,3 



exp ( 



RT 



(5.2-55) 


Using equation (5.2-56) with Catlow, et al.'s empirical results 
and equation (5.2-58) with Catlow, et al.'s non-emplrlcal results, 
Intrinsic oxygen Interstitial populations were calculated and the re- 
sults are given In table 11. At 1600* C: 
with Frenkel domination 

[0^] = 6.21xl0 -11 mol fraction 
= 6.21xl0' 5 ppm 
with Schottky domination 

[0^] » 1.79xl0 -3 ^ mol fraction 
= 1.79xl0' 25 ppm 

These calculated energies thus predict extremely low levels of ther- 
mal Interstitials. The number of extrinsic defects for the three 
donors, T1, H2, and H4, were determined as a function of dopant 
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TABLE 11. - ESTIMATES OF INTRINSIC OXYGEN INTERSTITIAL 
CONCENTRATIONS BASED ON THE CALCULATIONS 
OF CATLOW ET AL. 


Temperature Assuming Assuming 

( # C) Frenkel domination Schottky domination 

(empirical potential (non-emplrlcal potential) 


1100 

1.19 

E-14 

1.20 

E-42 

1200 

1.05 

£-13 

8.35 

E-40 

1300 

7.02 

E-l 3 

2.52 

E-37 

1400 

3.74 

E-12 

3.85 

E-35 

1500 

1.65 

£-11 

3.33 

E-33 

1600 

6.21 

E-l 1 

1.79 

E-31 

1700 

2.04 

E-10 

6.42 

E-30 

1800 

5.99 

E-10 

1.63 

E-28 

1900 

1.59 

E-09 

3.07 

E-27 

2000 

3.88 

E-09 

4.48 

E-26 


level, 10 to 1000 ppm, using an Iterative procedure and equation (21). 
The results are given In table 12 and shown In figure 132. 

The dependencies vary considerably depending on the binding 
energy. However at very low dopant level the degree of dissociation 
becomes determined by the concentration and the Interstitial concen- 
tration for all three donors Is within a factor of three. In all 
cases the concentration of free Interstitials produced by doping Is 

5 

vastly greater (factor of 10 for 10 ppm of the most weakly bound 
donor) than that of the thermally generated free defects. The effect 
of dopants Is not buffered away In the perfect singly doped case. 

- Heavy doping vs. Background Donors 

A second valuable comparison Is the difference In free Inter- 
stitials present In otherwise pure cyrstals containing a background 
donor level of 10 or 50 ppm and In a heavily doped, 800 ppm, crystal. 
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It Is clear from table 12 and figure 133 that the result Is 
dependent on the binding energy of dopant donor and, for the 50 ppm 
case, the binding energy of the background donor. Consider doping 
800 ppm of the three example donors. 

(I) A crystal doped 800 ppm of H2. In this case, the heavy 
doping will produce a population of free Interstitials large 
(340 ppm) compared to a background produced by any of the donors 
(at [FdO ? ] s 10 ppm the [0^] = 2 to 5 ppm; and at [FdO^] = 50 ppm 
the [0^ Is 3 to 25 ppm). 

Also, note that when doping with H2, across the whole range 10 
to 800 ppm, the [0^] = 1/2 [FdO^] as Is expected when there Is total 
dissociation. H2, then. Is an example of a donor whose effects can 
be accurately predicted when clustering Is neglected. 

(II) A crystal doped with 800 ppm of T1. In this case the heavy 
doping will produce a population of Interstitials which, while sub- 
stantial 70 ppm. Is much lower than H2. This level of Interstitials 
may be considered large compared to a 10 ppm background. However 
when the background Is 50 ppm this depends on the nature of the 
background. 

The [0^] = 70 ppm produced by the doppnt may be considered 
large with respect to the [0^] = 3 ppm produced by a 50 ppm back- 
ground of H4. However, In the case of a 50 ppm background of either 
T1 or H2, a significant fraction of the free Interstitials present 
(21 percent and 20 percent, respectively) will be produced by the 
background donor. 
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In addition, when doping with T1 near 10 ppm [0^] = 1/2 [Fd 0^] 
however at [FdO^] = 800 ppm this Is reduced to [0^] < 1/10 [FdO^J. 
Based on this, significant errors are to be expected In predictions, 
of the effects of Tl-doplng which Ignore clusters. 

(Ill) A crystal doped with 800 ppm of H4. In this case the 
dopant produced a level of free Interstitials, 9 ppm, which Is only 
slightly greater than that produced by a 10 ppm background and which 
Is more than a factor of 2 smaller than that produced by a 50 ppm of 
either T1 or H2. 

At high dopant levels, ~800 ppm, the [0^] = 1/100 [FdO^] which 
Is dramatically different than the unassociated case. Predictions of 
the effect of H4 doping are critically dependent on the binding 
energy of H4, and the binding energy of any background donor. 

- Donor doping with backgound acceptors and donors 

Additional complications occur when both acceptors and donors 
are present In the background. In the absence of association. When 
total donors exceed total acceptors, all of the acceptors can be com- 
pensated by equation (5.2-41) and a net donor population will exist. 

The percentage of each donor which exist In acceptor-donor 
clusters, donor-interstitial clusters, and as free Ions Is not neces- 
sarily uniform. For example, large acceptor Ions may prefer to asso- 
ciate with small donor Ions In order to have the average radius 
approach that of aluminum. The exlstance of this type of preference 
can lead to the segregation of various donors and acceptors Into 
different types of clusters. A strongly Ionizing donor added as a 
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dopant could conceivably be used up In the neutralization of 
acceptors and liberate a very weakly Ionizing donor. 

Therefore, In addition to knowing the concentrations of all the 
Impurities, It Is necessary to know the g-functlons of all of the 
possible cluster forms In order to make meaningful prediction of the 
defect chemistry. 

In summary, 

(I) Doping, despite clustering, produces free defect concentra- 
tions which are large compared to the population of thermally gen- 
erated defects. 

(II) Clustering diminishes concentration dependencies and there- 
fore the ability of dopants to swamp out the free defects produced 

by background Impurities. Also, the sensitivity of free defect con- 
centrations to cluster binding energies Indicates that Impurity 
levels below detectability limits can have significant contributions. 

(III) Predicting defect equilibria In crystals with 100 ppm 
background Impurities requires a great deal of Information which Is 
currently unavailable. That which Is available suggests much of the 
Impact of doping will be buffered out. Also, changes In the free 
defect population produced through control of the oxygen partial 
pressure with be similarly buffered, bound clusters will dissociate 
to provide defects sufficient to annihilate those produced by the 

P n change. 

2 

The second circumstance by which oxygen diffusion would be 
Independent of Impurity concentration. Is If the defect was both 
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neutral and Involved a stochlometrlc number of sites. Allovalent 
Impurities alter the concentrations of charged point defects In 
order to satisfy the constraints of charge neutrality and site 
balance. There Is no need for the population of neutral species to 
change with allovalent dopants. 

Both neutral Interstitials and stoichiometric neutral clusters 

have been proposed as dominating anion diffusion In a number of mate 

80 

rials. In a study of diffusion under a field In NaCl, Chemla 
found that a certain Impurity, Cs, obeyed the Nernst-Elnstlen rela- 
tion while anion Impurities P and S were affected to a much 

smaller, 1/10, degree than predicted by the relation. Similar work 

81 

on self-diffusion under an electric field by Nelson Indicated 
also that cation diffusion obeyed the relation while anion diffusion 
deviated from this expectation. Based on this, the suggestion was 
made that anion diffusion obtains a significant contribution from 
Schottky divacancies, ( V N a V Q-|) X - 

Studles 8 ^ -8 ^ of anion self-diffusion In III to V and II to VI 
compounds MX also lead to the proposition of neutral species, neutral 
Interstitials, S * In CdS and Te* In GaTe, as well as 
Schottky clusters P ^' 

Results In the oxides of MgO and Al^Og, as we ^« have 1ed to * he 
proposal of neutral species. The Insensitivity of oxygen diffusion 
In MgO has been attributed to the transport being governed by 
Schottky divacancies. 


Kroger and co-workers^ an< * ^ observed that there Is no grain 
boundary contribution to Ionic conductivity despite a significant 
enhancement of oxygen diffusion by grain boundaries. Based on this, 
grain boundary diffusion was proposed to take place by a netural 

atomic Interstlal 0?, or a neutral molecular Interstitial, 0« . 

1 *1 
Schottky clusters are another possible neutral species (recall 

from equation (5.2-20) that the concentration of these clusters Is 
uniquely fixed by the temperature). Although the population of un- 
associated Schottky defects Is small, the population of Schottky 
clusters need not be. It Is dependent on the binding energy of the 
clusters since 

- (E s + E b ) 

[(3V q 2V a1 )] « exp - (5.2-56) 

For the case of A1 2 °3’ th ^ s « ner 9y has not been calculated. 

It would be of great Interest to have such Information. Gourdln and 
Klngery calculated the binding energy for the Schottky divacancy and 
found It to be substantial, -257 kJ. If, In Al^, the pattern of 
clusters from allovalent solutes Is followed. In which k, < 10” 6 , 

O 

It would not be difficult to Imagine the diffusion of oxygen being 
governed by the cluster mechanism even If there Is a competing 
mechanism by a more mobile species, for example, V Q . This follows 
since. In the case where the defects Interact, It Is the product of 
the mobility and concentration of the species determine Its contrlbuton 
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5.2.3 P n Dependence 
u 2 

Let us consider the P n dependence of the three netural species. 

U 2 

0* , 0^. and (3V Q 2V A1 ) X . 

(1) In the case of 0~ , using the reaction: 

*1 

0 ? (g) - 0* (5.2-58) 

£ 1 * 


which will give: 


k « 


[0* ] 


(5.2-59) 


Therefore: 


to; ] a P f 


(5.2-60) 


(11) In the case of 0^ , the reaction will be: 


| 0 2 (g) - 0* 


(5.2-61) 


which will give: 


k = 


[ 0 *] 

„1 n 


(5.2-62) 


Therefore : 


r n x l pi /2 

C 0 1 ] a p o 2 


(5.2-63) 
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(111) In the case of (3V Q 2V A1 ) no P Q dependence will occur 
since neither k or k depends on P n . 

S a Up 

The experimental evidence, though restricted to one data point Indi- 
cates no dependence on P n and Is therefore consistent with the 

u 2 

Schottky Cluster model. 

In summary point defect clustering allows two separate explana- 
tions for the Insensitivity of oxygen diffusion to the presence of 
allovalent Impurities and changes In oxygen partial pressure. In 
the first the clusters are considered to have mobilities which are 
negligible compared to that of the free defects. The clusters serve 
as point defect traps and a buffering action Is responsible for the 
Insensitivities. In the second explanation clusters are considered 
mobile and the exslstence of Intrinsic defect clusters Is postulated. 
These clusters augment the thermally generated defect population and 
oxygen diffusion Is the result of Intrlnlc behavior. Further work Is 
needed to decide between the two. 

5.2.4 Non-FIcklan Tracer Profiles 

While the near surface region (<0.75 to 1.00 ym) could be ade- 
quately described by a solution to Flcks second law, the profiles 
contained, to varying degrees, low level deeply penetrating tails, 
extending past the furthest depth probed, ~2ym In proton activation. 
Typically the concentration In the tall was 2 to 7 times the natural 
abundance. The tails appeared In both the proton activation spectra 
and SIMS profiles Indicating that they are not artifacts of the 
analysis method. 
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Slmllar tails were observed by Reddy 16 . His results may be 
summarized as follows: 

1. All as-received samples showed this deviation from Pick’s 

law. 

2. The magnitude of the tall did not correlate to the average 
native dislocation content. His crystal with the lowest dislocation 
density had the tall of the highest magnitude. 

3. Preannealing the crystal affected the magnitude of the tall 

but not In a consistent manner. Three samples were preannealed In a 
-2 -3 

vacuum (10 to 10 torr), two of these samples had tails, one 

did not. Samples which had been preannealed In an oxidizing 

atmosphere were reported to have no tails, or tails of lower 

magnitude. (Note that In the present work all samples, exchanged at 

P = latm were preannealed In air, despite this tails occurred). 

U 2 

Reddy concluded that the tails did not occur by short circuit 
diffusion alone since the tall did not correlate with variations In 
dislocation density among his samples. An unreasonable pipe radius, 
on the order of a micron, had to be assumed to account the level of 
enrichment observed. He proposed two explanations. The first was 
based on the presence of excess vacancies quenched In during crystal 
growth. In this case the diffusion of tracer Into the sample Is en- 
hanced along a dislocation and then Is enhanced radially away from 
the dislocation due to the wind of excess vacancies annihilating at 
the dislocation. The lack of dependence on dislocation density was 
ascribed to differing concentrations of excess vacancies dominating. 
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Were this effect to be the source of the tails then samples 
which were preannealed and diffusion annealed together would be 
expected to show similar, though nonflcklan, profiles. However, 
tall magnitudes In the current work are seen to vary between such 
samples. Figures 133 and 134 show the proton activation spectra 
from two sets of samples, two from the Nl-doped boule and two from 
the undoped no. 1 boule. All four samples were preannealed 
together. During the exchange anneal the samples were split Into 
two pairs, of one sample from each boule, located on two different 
shelves In the exchange apparatus. In both cases, the samples from 
the same boule had tails of differing degree. The difference In 
temperature between the two shelves, ~40° C, Is not the source of 
the variation since the high temperature sample had the more 
pronounced tall In one case while the low temperature sample did In 
the other. In summary, tall magnitude varied within groups of 
equivalently equilibrated samples. Note also that atmosphere 
effects cannot be responsible (for example, amount of C0 2 or 
N 2 which leak Into the apparatus) since all these were exchanged 
In the same atmosphere. 

The second proposition offered by Reddy Is the coexistence of 
two noninteracting oxygen species, molecular oxygen Interstitials 
and oxygen lattice defects. Such a situation has been proposed to 
explain the oxygen transport across vltrlous silica films (see 
Appendix A). However, there are two reasons for Its rejection In 
Al^Oj. The first Is the expected P Q dependence. The concentration 
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of dissolved molecular oxygen Interstitials Is linearly dependent on 
the oxygen partial pressure. However, the tall In the specimen ex- 


10 


-15 


atm. had a tall comparable to the samples 


changed at a P Q 

c 

exchanged at P. =1 atm. Secondly, experiments on the permeation of 
U 2 


sintered alumina tubes (99.8 percent Al^) by Hayes, Budworthy, and 
Roberts found the material permeable to oxygen yet not to argon or 
nitrogen. If the oxygen was permeating as a noninteracting molecular 
species, argon with Its smaller molecular radius (Ar - 2.88 A, 

- 2.98 A) would also be expected to be transported. 

The coexlstance of contamination and Impurity-defect clustering 
offers another possible source for tails. When the cation Impurities 
are diffusing In they may drag with them an oxygen defect. Donors 
may conceivably drag noninteracting Interstitials. 

In this case the flux of tracer would simply be the sum of the 
flux due to self-diffusion and one-half of the flux of the contami- 
nant, since 

[0^] = 1/2 [(2FdO i ) x ] 


with the assumption that the diffusion coefficients are not functions 
of position, the solution for a semi-infinite slab Is 

c lfl (x,t) = c, erfc ( ~z) * 1/2 c.c,.. erfc ( — ■ - — ) 

2 vDt ru I t 

<T 2 V0(2Fd A1 0 ) ) t 

where c^ Is the concentration of tracer In the lattice at the sur- 
face (equals gas for K ■* ®) and c^ Is the concentration of the 
donor Impurities In the lattice at the surface. 
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In this case the magnitude of the tall Is directly related to 
the concentration of Impurity and an observable tall would require a 
large level of contamination. 

For example, a tall three times the natural abundance, or 0.6 

percent, when using an atmosphere of 50 percent enriched 0-18, would 

require more than 2 percent of the cation sites to be occupied by 

Impurities. This seems unreasonable given the solubility limits of 

25 

~200 ppm given for donor Impurities. It would be possible for 
tracer to accumulate If the donor-interstitial cluster were a very 
much faster diffusing species and Interacted Infrequently with the 
lattice. However, a sample which was demonstrated to have cation 
profiles, VD-LT (exchanged at 1300° C), did not have a tall In Its 
profile. 

It remains to find a satisfactory explanation for the tails. 


CHAPTER VI 


CONCLUSIONS 

On the basis of this work the following conclusions can be made: 

1. ) SIMS has several advantages over proton activation for 
tracer studies. That of greatest Importance Is the superior depth 
resolution and the ability to detect a greater number of species. 

2. ) The contamination of the samples was the result of the use 
of commercial high purity Al^ tubes. 

3. ) The oxygen dlffuslvlty determined was In agreement with that 
determined In previous gas exchange studies In absolute magnitude, 
although the determined activation energy was lower. 

4. ) Two models for oxygen diffusion In alumina are consistent 

with the diffusion coefficient being Independent of Impurities and 

P . The first Involves a buffer action arising from Impurity point 
u 2 

defect clustering and the second Involves the exlstance of Schottky 
Clusters, (3 Vq 2V a ^) x . Both models can be justified by extrapolation 
from the defect energy calculations of Catlow, et al. 
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CHAPTER VII 


Suggestions for Future Work 

1. ) Contamination by furnace tubes was a serious problem In this 
work. The design and construction of an apparatus which will mini- 
mize the level of contamination Is essential. 

2. ) A study of diffusion under a constant applied electric 
field would be of value. Diffusion which occurs by neutral Schottky 
Clusters should be unaffected by this field while Ionized point de- 
fects would be. A gas exchange experiment using porous platinum 
electrodes Is one possible scheme. 

3. ) It would be of Interest to carry out another set of experi- 
ments using a solid state source of tracer. It Is Important to 
establish that diffusion coefficients are not a function of par- 
ticular experiment techniques. 

4. ) SIMS with Its excellent depth resolution should allow direct 
measurement of diffusion coefficients at Tow temperature where only 
Inferred values currently exist. 

5. ) It would be of great Interest to have the calculations of 
Catlow, et al., extended to Include various cluster orientations, 
clusters of Intrinsic defects and the energies of motion for both 
point defects and clusters. 
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APPENDIX A 


The Thermal Oxidation of Silicon 


I. Introduction 

Beyond Its general scientific Interest, the study of the passive 
oxidation of silicon Is of commercial value. The thermal oxidation 
of silicon plays a key role In MOS technology and, since oxygen 
transport through S10 ? Is Involved In the oxidation process, know- 
ledge of the mechanlstlcs Is helpful In the understanding of glasses. 

A number of different mechanisms have been proposed for the 
thermal oxidation of silicon. Also, there are Indications that the 
mechanism of oxygen transport In vitreous silica (v-S10 2 > Is complex. 
This work was undertaken to Identify different contributions to 
oxygen mobility and to assess the contribution of each In thermal 
oxidation. 

It Is found that oxygen Is mobile both within the network and 
as a permeating species. Network diffusion, however, does not 
contribute to scale growth. 

II . Literature Review 

(1) Silicon Oxidation Kinetics 

The thermal oxidation of silicon has been shown to produce a 
scale of v-SIO^ and the kinetics of this oxidation has been shown to 
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n 21 

be of a linear-parabolic type' ’ . With this type of kinetics, the 

thickness of the scale after time, t, Is given by: 
x 2 ♦ Ax = Bt 

where x Is the oxide thickness, B/A Is the linear rate constant 
K, , and B Is the parabolic rate constant K . 

1 p 

Linear oxidation kinetics are expected when the rate of oxida- 
tion Is controlled by the reaction of the oxidant 1 with the 
metal. Parabolic kinetics, on the other hand, are expected to occur 
when the oxidation reaction Is controlled by the diffusion of the 
re- agents across an existing scale. Oxidation which Is Initially 

controlled by reaction rate may undergo a transition to parabolic 

( 2 ) 

kinetics as the scale thickness Increases. Irene has shown 
that the values of both the linear and the parabolic rate constants 
are functions of water content and silicon substrate orientation. 

The nature of the dependence of the parabolic rate constant on sub- 
strate orientation Is unclear. It Is not obvious that diffusion 
through the v-SIO^ should be affected by the parent silicon 
orientation. 

Investigations of the kinetics of the growth of thin scales 

( 3 ) 

where linear kinetics dominate, lead Ghez and Van der Meulei to 
conclude the two simultaneous reactions are responsible for oxlda 
tlon. They propose the following set: 

1 An oxidant Is a chemical reagent that oxidizes. Oxygen Is one 

example but there are many others, e.g., water. 



-277- 


(I) Silicon reaction with molecular oxygen 
SI - SI + 0 ? -» SI - 0 - SI ♦ 0 

(II) Molecular oxygen dissociation and silicon reaction with 
atomic oxygen 

°2 l 20 

SI - SI ♦ 0 -* SI - 0 - SI 

(A) 

Irene' ' analyzed the curvature In an Arrenhlus plot of the linear 

rate constants and also came to the conclusion that two simultaneous 

reactions occur. However, the reactions were not specified. 

( 5) 

As reported by Major' the values for the parabolic rate con- 
stant and Its temperature dependence showed wide variation. The 

( 2 ) 

sensitivity to the ambient water content shown by Irene' ' may con- 
tribute to the scatter. Deal and Grove^ also suggest that the 
apparent discrepancies may be due to force fitting data with a para- 
bolic law when It Is Inappropriate. 

(2) Silicon Oxidation Mechanisms 

Marker experiments by Atalla^ and Jorgenson^ both Indicate 
that the growth of the oxide takes place at the oxide-metal Inter- 
face. The linear dependence of the growth rate. In the parabolic 
regime, on oxygen partial pressure was observed by both Llgenza and 
Spltzer^ and Deal and Grove^ which suggests that the oxidant Is 
molecular oxygen. 

Llgenza and Spltzer also observed that samples successively oxi- 
dized In oxygen-16 and oxygen-18 contained more oxygen-18 than could 
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be accounted for by scale growth. This Indicates that the oxygen 

within the network undergoes exchange with oxygen In the ambient gas. 

Jorgenson^ found that when a battery circuit was used to apply 

a field across a scale during growth. It was possible to enhance or 

retard the growth rate. This was asserted to Indicate that the 

( 2 ) 

transported species was Ionic. Based on this, Irene proposed that 
the oxygen molecules were doubly charged and Deal and Grove proposed 
that they were singly charged. Raleigh' ' however, has reinterpreted 
Jorgenson's work and pointed out that oxidation Involves the trans- 
port of electroneutral species and therefore an electric field cannot 
provide a steady state driving force. Rather, Raleigh suggests that 
the battery circuit provides an external path for charge which allows 
Ionic species which were previously limited by electrical charge flux 

to participate In the oxidation. 

(5) 

In 1978, Major carried out a set of experiments above the 

temperature range of 900° to 1200° C In which silicon coupons were 

oxidized In oxygen-18 for a time to produce 1/2 vim thick SIO^. 

The atmosphere was then switched to oxygen-16 and the oxidation was 

allowed to proceed until a full micron of oxide was grown. The 

location of the tracer was determined using proton actlvlatlon. The 

2 

results are shown In figures 1 to 4. . At low temperatures, 900° 

C, the location of the bulk of the second 1sotope(16) Is located 
near the oxide-metal Interface while at the high temperature, the 

^Note that the Isotope order In these spectra Is the reverse 
(0-18), then 0-16) of that employed In all the other studies 
mentioned . 
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X (pm) 


Figure 1. - Best fit concentration profile of Major' 
oxidized at 900° C. 


1 


sample 


Cone, of 2nd Isotope 
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Figure 2. - Best fit concentration profile of Major's sample 
oxidized at 1020° C. 
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X (pra) 


Figure 3. - Best fit concentratl 
oxidized at 1140° C. 
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Figure 4. - Best fit concentration profile of Major's sample 
oxidized at 1200° C. 
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second 1sotope(16) Is located near the gas-oxide Interface. At 
Intermediate temperatures the second 1sotope(16) Is distributed 
between the two Interfaces. In order to account for this shift a 
model was proposed In which oxygen molecular Inter- stltlals 
permeate the scale to react with the metal. During the permeation 
these Interstitials were thought to exchange with the network 
oxygen. If the rate of lattice exchange was small at 900° C the 
tracer would be deposited at the oxide-metal Interface. If, as the 
temperature Increased, the rate of exchange with the lattice also 
Increased, then the tracer would be found Increasingly toward the 
gas-oxide Interface. Clark^ 10 ^ has developed a numerical solution 
to this model. 

Also In 1978, Rosencher, et al./ 11 ^ carried out a double oxi- 
dation experiment at 930° C. They also used an activation technique, 
however It was a resonance technique with better depth resolution. 
With this technique, the tracer profiling (see fig. 5) was found to 
consist of a large block of tracer near the oxide-metal Interface 
and a narrow spike of tracer at the gas-oxide Interface. This Indi- 
cate that, at least at 930° C, the majority of the tracer was In- 
corporated Into the scale at the oxide-metal Interface, suggesting 
growth by a permeating species. 

{121 

In 1981, Christy and Condon' ' carried out a double oxidation 
experiment at 700° C, analyzing their samples with SIMS. The result- 
ing profiles have features which suggest that artifacts exist. The 
profiles are shown In figure 6= The oxide-metal Interface Is re- 
ported to be around 3 min. of sputter time. The total oxygen signal 
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Figure 5. - Typical 18 0 concentration profile measurements using 
the 2-keV-slde 18 0(p,a) reaction near the 629-keV 
resonance. Channel 0 corresponds to 619.6 keV; the energy step 
Is 384 ev/channel. The solid lines show the best fit calculated 
assuming concentration depth profiles represented In the 
Insets, the depth scale origin being the external surface of 
S102- The concentrations In the extreme regions were taken 
constant and equal to that of the gas, 99.7 percent, whereas In 
the bulk It was taken equal to natural, 0.2 percent. Respective 
equivalent thicknesses of these regions are Indicated (see Table 
II), from ref. 11 . 
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Oepth (nm) in Si Metal Equivalencies 
36 89 142 195 248 301 354 



Sputtering Time (mini 


Figure 6- - Typical IMMA depth profiles of oriented silicon wafer 
which had been oxidized at 700° C with 16 02 for 10 days 
followed by IBOj for 10 days. (Actual oxide depth Is about 
twice the silicon metal equivalence depth. The location of x s , 
the defined transformation boundary. Is at about 3 min sputtering 
time.), from ref. 12. 
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appears to vary across the SIO^- This Is a reminder that SIMS Is 

not a quantitative Instrument without standards. After the oxide- 

metal Interface the profile shows that the concentration of oxygen 

does not decrease rapidly. This might arise If the Interface were 

rough or If there Is an appreciable solubility of oxygen In silicon, 

( 13 ) 

however, as pointed out by Helms In his discussion neither Is 
the case. The shape of the profiles are consistent with the edge 
effects or a poorly focussed primary beam. Though excellent depth 
resolution Is claimed no standards were presented and on the basis 
of the poor profiles, the model presented In their paper Is con- 
sidered to be unproven. 

In 1982, Revesz and Schaeffer ^ 1 ^ proposed that the v-S10 2 
scale contains "channels" of Increased structural order similar to 
trldymlte. Their model considers scale growth to arise from two 
processes, network diffusion and permeation by molecular Inter- 
stitials. The channels are proposed to exist In sufficient number 
that transport phenomena are dominated by them but also at a low 
enough concentration that they are not detectable by electron or 
X-ray diffraction. It Is also stated that permeation Is greater In 
channels while network diffusion Is decreased. Further, It Is stated 
that slow oxidation or heat treatment of the scale will Increase the 
level of ordering and a corresponding Increase In permeation and a 
drop In network diffusion will occur. Note that It Is necessary for 
these effects to balance each other for parabolic kinetics to be 
observed. The credibility of this model would be greatly enhanced 
If there were diffraction (or perhaps charged particle channeling) 
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evldence to confirm the structural changes. 

(3) Diffusion of Oxygen In v-S10 2 

There have been five studies of the diffusion of oxygen In 
v-S10 ? . 

In 1961, Norton^ ^ studied and permeation of gaseous oxygen 
through a membrane of v-S10 2 and found the diffusion coefficient 
could be described by: 


„ _ , „-8 .-113 kJ/mol. 2. 

D = 2.88x10 exp ( ^ ) m /$ 

In 1962, Haul and Dumbgen^ 16 ^ studied the rate of uptake of 
tracer by S10 2 and found that the diffusion coefficient could be 
described by: 


D = 4.3xl0" 6 exp ( ~ ?34 R !^~ o1 ) m 2 /s 


In 1963, Sucov^ 17 ^ monitored the tracer uptake from a finite 
source by tiny silica cylinders and found that the diffusion co- 
efficient could be described by: 


, , , „-2 .-297 kJ/mol. 2. 

0 = 1.5x10 exp ( jpf ) m /s 


In 1965, Williams' ' also monitored the tracer uptake by 
either silica fibers or bulbs and found the diffusion coefficient to 
be described by: 


0 = 2.0x10 


-13 


exp ( 


-113 kJ/mol 
RT 


? 

m /s 


In 1982, Ylnnon^^ studied diffusion of oxygen In v-S10 2 by 
exposing samples to a tracer enriched gas and then determining the 
tracer profile by proton activation. This procedure has the advan- 
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tage that the determined value for the diffusion coefficient must 
satisfy both the quantity of tracer In the sample and the shape of 
the profile. It was found that this was true for samples exchanged 
at temperatures below 800° C. At higher temperatures, however, the 
profile could not be fitted with a single diffusion coefficient; but 
could, be described by two diffusion coefficients (see fig. 7). 

This was suggestive of two Independent diffusion mechanisms. The 
near surface regions of their profile yielded diffusion coefficients 
which could be described by; 


~ „ , ,„-14 .-110 kJ/moK 2. 

D = 2.1x10 exp ( ^ ) m /s 


In summary, the range of reported values for the diffusion of 
oxygen In v-S10 2 varies over a wide range. The work of Ylnnon and 
Cooper Indicates that the presence of more than one operative diffu- 
sion mechanism, which may be sensitive to P u n , surface condition, 

rlpU 

and temperature, may contribute to the scatter. 


III. Experimental Procedure 
(1) Samples 

The single crystal silicon samples used In this study were 
oriented In the (001) direction and were p-type with a resistivity 
of 12.4 ohm-cm. 

Prior to experimentation the silicon was Immersed In an HF solu- 
tion to etch off existing oxide, washed with acetone, then methanol, 
and finally rinsed with distilled water. 
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(2) Thermal Oxidations 

The apparatus described In the main text section for oxygen 
diffusion In AI^O^ was also used to carry out the thermal oxidations. 
The first oxidation was carried out either In tank oxygen which had 
been passed over drlerlte for 24 hours at room temperature or In air 
with the relative humidity measured with the wet bulk/dry bulb tech- 
nique. The second oxidation was carried out using the tracer stored 
In the zeolite flask, which was also passed over drlerlte for 24 
hours. Typically the sample was cooled to room temperature between 
the two oxidations. This was not desired but rather due to unavoid 
able equipment limitations. 

IV. Results 

(1) Qualitative Interpretation 

Two samples were oxidized simultaneously on separate shelves at 
different temperatures, 967° and 1055° C. The oxidation In natural 
oxygen took place for 3 hr. at 25 min. and the tracer oxidation for 
15 hours. 

These samples were analyzed using proton activation and the re- 
sulting spectra are shown In figures 8 and 9. It was found that 
these spectra could only be fit with a concentration profile similar 
to those assumed by Rosencher, et al., which consisted of the bulk 
of the tracer near the oxide-metal Interface and a very narrow pro- 
file near the oxide-gas Interface, (see figs. 10 and 11). The fit 
was somewhat Insensitive to assumed concentration profiles, that Is, 
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Figure 8. - Proton Activation spectrum of silicon doubly oxidized 
at 967* C. 
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the depth resolution was poor enough so that small scale changes In 
the profile were not reflected In the spectrum. 

One of these samples was re-examined with SIMS. The profile of 
oxygen-18 and oxygen-16 Is shown In figure 12. Due Is the much 
better depth resolution of this technique the two regions of tracer 
concentration are distinct. The concentration of tracer In both Is 
very near that of the gas, ~4Q percent. A region corresponding to 
the natural abundance of tracer can be seen toward the middle of the 
scale. 

The numerical solution by Clark based on the model of Major and 
Cooper, does predict near surface enrichments for particular ratios 
of diffusion to exchange however It does not predict the midscale 
region of natural abundance to occur at these ratios. With this 
region of natural abundance the model would predict no tracer near 
the oxide-gas Interface. Note also that In the SIMS spectra the 
0-18 and 0-16 signals are complementary and the S1/S10 2 Interface 
Is sharp. 

Three samples were simultaneously oxidized In air at 1000° C 
for 2 hours. Subsequent to this they were Individually oxidized In 
a tracer enriched atmosphere for times of 3.5 min. 1 , 11.5 min., 
and 37 min. The resulting SIMS profiles are shown In figures 13 to 
15. By Inspecting these profiles It Is clear that tracer Vs 
accumulating at both the gas-oxide and the oxide-metal Interfaces. 

^These values are approximate because at these short times, the 

time required to bleed In and remove the tracer Is not negligible. 


o o 
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Figure 12. - SIMS profile on a sample oxidized for 3.4 hr In 
0 16 and for 15 hr In 0 18 at 967° C. 
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Figure 14. - SIMS profile of sample 2/11. 






GO 


Figure 15. - SIMS profile of sample 2/37. 
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These regions are separated by a large region devoid of tracer. 

This Implies that they occur by distinct mechanisms. Ordinary gas 
exchange and self diffusion can give rise to a tracer profile of the 
type which Is observed near the gas-oxide Interface while the 
profile which occurs oxide-metal Interface may arise from permeating 
Interstitials which do not frequently Interact with the network. 

Thus from this qualitative Inspection there appears to be two 
Independent mechanisms for oxygen transport, network diffusion and 
permeation. 

Secondly, the peak height of the deep profile In the 11.5 min. 
and the 30 min. experiments (figs. 2 and 3) Is greater than that of 
the near surface profiles. When permeation Is the dominating con- 
tributor to scale growth and the concentration of tracer at the 
oxide-metal Interface would be expected to accurately reflect the 
concentration of tracer In the source of gas, since all new oxide Is 
formed with this composition. Any network contribution to growth at 
this stage would dilute the concentration of tracer at the oxide- 
metal Interface since the network oxygen adjacent to this Interface 
Is of lower tracer concentration. The fact that the profile at the 
gas-oxide Interface does not attain the value of tracer In the gas 
Indicates that the kinetics of the reaction between the gas and 
solid play, at least, a limited role. Unfortunately the concentra- 
tion of tracer In the gas was not known precisely enough to allow a 
quantitative determination of the surface reaction coefficient. 

Note when the tracer exchances take place for longer times, for 
example 8 hours as In figure 16, the gas-oxide and the oxide-metal 



00 


SIMS profile of sample 
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interface profiles have very near the same tracer concentrations 
Indicating that the surface reaction coefficient Is most Important 
at short times. 

(2) Mathematics of Double Oxidation Experiments 

The thermal oxidation of silicon has been shown to follow 

f 21 

linear-parabolic kinetics' , that Is, at short oxidation times the 
thickness of the oxide scale Is directly proportional to the oxida- 
tion time while at long oxidation times the thickness of the scale 
Is proportional to the square root of the oxidation time. A basic 
analysis of linear-parabolic kinetics when the oxidation reaction 
occurs at the metal oxide Interface has been given by Deal and 
Grove^ . What follows Is essentially a restatement of their work. 

It should be noted that this treatment assumes there Is only one 
transport mechanism across the growing scale. 

Consider a semi-infinite slab of silicon covered with an Initial 
film of S10 2 sufficiently thick that steady state has been reached 
(see schematic shown In fig. 17). In this case the oxidant must go 
through the following stages: 

(1) It Is transported from the bulk gas to the gas-oxide 
Interface where It reacts or Is absorbed. 

(11) It Is transported from the gas-oxide Interface towards 
the oxide-metal Interface. 

(Ill) It reacts at the oxide-metal Interface to form SIO^. 

The assumption of steady state means these three fluxes are equal. 

The steady state fluxes are mathematically modeled as follows. 
The flux, J.j , of the oxidant from the gas to the surface of the 
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oxide Is: 

J 1 " K(? e ' ^o ) (A_1 

where K Is a surface reaction coefficient, c g Is the equilibrium 

concentration of oxidant In the scale and cf Is the concentration 

o 

of oxidant In the scale at X - 0. Note: c should not be con- 

o 

sldered to be normal 0 e In network. 

The flux of oxidant across the oxide layer, given by Pick's 
First Law: 


J 


2 


-0 


3c 

eff dx 


( A-2) 


where 3c/3x Is the concentration gradient of oxidant across the 

scale and 0 ,, Is the effective diffusion coefficient which Is not 
erf 

necessarily equal to the self diffusion coefficient, as will be dis- 
cussed. It follows from the assumption of steady state that the 
flux, J ? , Is not a function of position within the scale. When 

Is also Independent of position the consequence Is that ac/3x Is 


also a constant and equation (A-2) may be rewritten as: 


< c o ‘ c 1> 


'2 eff (X 2 - X^ 


( A- 3 ) 


The flux corresponding to the oxidation reaction Is given by a 
first reaction: 

J 3 ■ k s c, (A-4) 


Also from steady state It may be stated that: 


J = J ] = = J 3 


( A- 5 ) 


-303- 


From these equalities It Is possible to solve for and C 0 : 


C 1 " 1 ♦ k $ /K ♦ k s (X 2 - X 1 )/D eff 


(A-6) 


. ? e (l ♦ [k $ (X 2 - X,)^]) 
: 0 " 1 ♦ k $ /K ♦ k $ (X 2 - X-j )/0 e ^ f 


C A— 7 ) 


It Is easy to see by back substitution of A-6 Into A-4 that: 


k s c e 


J " J 3 “ 1 ♦ k $ /K ♦ k $ (X 2 - Xp/D^ 


( A— 8 ) 


Now, the rate at which the oxide grows Is given by: 


N 1 dt (X 2 " V 


k s c e 


1 ♦ k s /K ♦ k s (X 2 - X 1 )/0 eff 


( A— 9 ) 


where N. Is the molar density of oxygen In the scale. With the 

I i 

Initial condition that: 

(X (t) - X (t)) « X 

£ t*0 

where X^ Is an arbitrary Initial thickness of oxide and Integra- 
tion, the following expression Is obtained: 


(X 2 - ♦ A(X 2 - X ] ) = B(t ♦ x) 


(A-10) 


2 D 


eff 



where 
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2 0 


B = 


eff 


and 




+ Ax^ 
B 


Note t corresponds to a shift In the time coordinate to account for 
the Initial oxide thickness, x^, solving for (X ? - gives 


(X 2 - X,) - f { 1 ♦ ( 

for short times such that 
kinetics where: 


t ♦ t 

A 2 /4B 


) } 


(t ♦ t) 


1/2 A 
2 

« A 2 /4B 


( A- 1 1 ) 

this reduces to linear 


(X ? - X 1 ) - | (t ♦ t) ( A - 1 lb) 

2 

And for long times such that (t ♦ i) » A /4B this reduces to para 
bollc kinetics where: 


(X 2 - X^ = B 1/2 (t ♦ t) 


1/2 


( A - 1 I c ) 


The above analysis of Deal and Grove may be applied to the 
analysis of a "double oxidation* experiment. In this type of experl 
ment a semi-infinite slab of silicon Is thermally oxidized for a 
given time t to give an Initial oxide thickness X^. At time i 
the Isotopic makeup of the atmosphere Is altered. In this study by 
Increasing the Isotopic fraction of oxygen-18, and the oxidation Is 
allowed to continue. After the oxidation Is terminated, the location 
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of the tracer Is determined by such techniques as proton activation 
or SINS In order to give Insight Into the underlying mechanisms. 

It has been shown that tracer enrichments develop Independently 
at both the gas-oxide Interface and the oxide-metal Interface. The 
Implication of this Is that oxygen may be transported across the 
scale by two Independent mechanisms. In the following discussion 
oxygen will be considered to be mobile through the network by a de- 
fect mechanism and also mobile through noninteracting Interstitials. 

The expression for the Incorporation of tracer when convection 
plays a role Is: 

J 18 ■ -° £ * vc <*-'*> 

where c Is now the fractional concentration of tracer and v Is 
the convective velocity. 

There are two separate types of contributions to the velocity 
term. One occurs due to a moving boundary or flow of the system, 
while the other one results from the oxygen potential gradient across 
the scale. The latter Is analogous to vacancy "wind* under a field 
considered by Nanning. Equation (A-12) may be rewritten as: 


J 18 “ -° S d S ‘ 


c - 


°eff c ajj 
RT ax * v mb c 


< A- 13) 


where 0 S(J Is the self diffusion of oxygen, Is the amblpolar 

diffusion coefficient responsible for oxidation, c Is the concentra- 
tion of oxygen-18, c Is the concentration of the defect responsible 
for network diffusion, dy/dx Is the chemical potential gradient In 
oxygen, and v^ Is the velocity of the moving boundary. 
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It Is convenient to consider the two profiles within each sample 
Independently and use two sets of coordinates. For the profile near 
the gas-oxide Interface the observer A In figure 18, which Is fixed 
to the gas-oxide Interface Is appropriate and the coordinate system 
x Is used. In this coordinate system X, Is defined as zero and 

d I 

x Is positive towards the oxide-metal Interface. 

The convective velocity due to a moving boundary seen by ob- 
server A Is zero since all growth takes place at the oxide-metal 
Interface and A- 1 3 reduces to: 


ac °eff c a u 

— _ n + r • -- — — 

18 so ax c RT dx 


( A - 1 4 } 


Now let us examine the near surface profiles developed In two 
samples, one of which was oxidized In air for 16 hours and the other 
for 48 hours. Both samples were then oxidized In a tracer enriched 
atmosphere for 8 hours. 

Figure 19 shows a superposition of the SIMS profiles from these 
samples. The near surface profiles are seen to be coincident. The 
tracer concentration profiles are Identical for both. Since the self 
diffusion term Is uninfluenced by oxide thickness, the convective 
terms from the chemical potential gradient must also be equivalent 
for the two samples to have the same tracer fluxes. 


c 



dp. 


ax, 


c 




With the assumptions of Ideal solution, steady state and parabolic 
kinetics this equality may be rewritten: 



-308- 



Figure 19. - Superposition of SIMS profiles from 2 samples, one 
oxidized for 48 hr In air the other for 16 hr. Both then were 
oxidized for 8 hr In 18 0o, at 1000° C. 
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- Ideal solution 

v « RT In c 

3 jj rt ac 
ax “ ~ ax 

- steady state 

ac (? o - Ii 
ax “ x 2 - x 1 


- parabolic kinetics 



. 1/2 


(t. 


T 1 


, 1/2 


c • 0 


< c o - c 1> 


eff B 1/2 (t, ♦ V 1/2 


« c • 0 


< c o - c 1> 


•" b' / 2 <t 2 . , 2 ) 1/2 


1/2 

The only experimental variable In this term Is the time, (t ♦ t) 

When t 1 ♦ ^ t 2 ♦ t 2 the only way these may be equal Is If a con 

stant has a value of zero or Infinity. Infinite value of any of 
these may be rejected on physical grounds. Because the value of 

c Q and e i are fixed at the edges of the scale and the parabolic 


rate constant Is non zero It follows that the effective diffusion 
network coefficient Is zero. This Is a bit surprising since oxygen 
Is mobile In the network. It would be expected to travel down the 
oxygen activity gradient. However, as previously mentioned. It Is 
pointed out by Raleigh that the motion of Ionic species across the 
network must Involve the transport of electrical charge. Raleigh 
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cites experimental work of Schmalzrled which Indicated that at 
1000° C S10^ Is a pure Ionic conductor. Further the work of 
Jorgenson, Interpreted by Raleigh, Indicates that when an externa) 
path for electrons Is provided, an Increase In the oxidation rate Is 
observed. In sum this Indicates that at temperature 1000° C the 
electrical conductivity of S 1 0^ Is sufficiently low that no network 
diffusion of an Ionic species Is permitted to relieve the activity 
gradient . 

The superposition of the profiles In figure 7 also provide 
an opportunity to test a prediction of the model of Revesz and 
Schaef f er^ 1 . Recall that this model predicts that the network 
dlffuslvlty will decrease as annealing time Increases, and that the 
permeation rate will Increase. The coincidence of the self diffu- 
sion profile Indicates that the network diffusion coefficient Is 
unaffected by the 32 hour difference In annealing time prior to 
exposure to tracer. 

The quantity of tracer which was Incorporated Into the scale 
may be accurately estimated by taking the distance from the oxide- 
metal Interface to the point at which the tracer concentration drops 
to 1/2 of Its plateau value. This Is Illustrated In figure 19. The 
experimental value of 1.68 compares very well with a value of 1.62 
calculated on the basis of parabolic kinetics Indicating that the 
rate of permeation Is also unchanged. The prediction of Revesz and 
Schaeffer Is not borne out. 

Now consider the flux of tracer Into the oxide near the oxide 
metal Interface. The fact that the tracer Is being carried by a 
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specles which does not Interact with the network until reacting at 
the oxide-metal Interface allows the specification of a boundary 
condition at X^. The background produced due to the presence of 
a permeating species Is In the ppm range and may be neglected. In 
this case It Is convenient to use a coordinate system fixed to the 
oxide-metal Interface with X^ and X, seen as positive and moving. 

From the point of view of observer B, x ■ 0 at X ? , there Is a 
convective flux of oxygen that Is equal to the growth rate of the 
scale given by: 




it <?' > * < 2 


+ x 1/2 A 

- “) } - f) 


AV4B 


*75 (i ♦ (V -1 )) 

*' D AV4B 


- 1/2 


(A-16) 

For long times when the oxidation may be approximated by parabolic 
kinetics: 



V 


3 _ r 1/2 * J /2 

at 8 (t * t) 


B 1/2 
2<t ♦ t) 


1/2 


(A-l 7) 


Therefore the flux of tracer Into the oxide near observer B Is: 


J 18 * ~°sD dx * c 


B 


1/2 


2<t ♦ t) 


1/2 


( A-l 8) 


In partial differentiation form the diffusion equation (A-12) takes 
the form: 


dc 3J n a 2 c w ac 
9t ■ - ii ■ 0 ^2 - y ii 


(A-19) 
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Thls equation has been solved (see Appendix C) for both the case of 


constant v and 

the case 

where 

V 

has a time dependancy of 

under the simple 

boundary 

conditions 


c(x,t) = 1 

X = 0 

and 

t > 

0 

c ( x , t ) = 0 

X -* OO 

and 

t > 

0 

and the Initial 

condition 




c(x,t) * 0 

0 < X 

< oo 

and 

t = 0 


However, the double oxidation experiment requires a moving boundary 

- 1/2 

with a time dependancy of (t + const) . The solution to the par- 
tial differential equation with this type of coefficient remains ellu 

slve. When t Is small compared to t It may be neglected and v 

- 1 17 

has the character of t . When -t Is large compared to t the 
value of v may be taken as a constant. 

Three samples were fit using solutions to the diffusion equation. 
One sample was oxidized at 960° C under 0-16 for 3.5 hours and the 0 18 
for 15 hours. Two samples were oxidized at 1000° C. One was oxidized 
under 0-16 for 16 hours and the other one was oxidized for 48 hours. 
Then both were oxidized together under tracer for 8 hours. 

Other samples such as those shown In the qualitative discussion 
were not In the regime of the parabolic kinetics and quantitative fits 
were not attempted. 

The predicted profiles were obtained by superimposing a near sur- 
face and a deep profile. In all cases the near surface or gas-oxide 
Interface profile was generated using assumed values for the diffusion 
coefficient In the familiar solution for a semi-infinite slab 
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c(x,t) = c erfc (— *—) 

0 2Vot 

For the deep or oxide-metal Interface profiles, two different 
solutions were used. It was found that for the samples exchanged at 
1000° C the use of average velocity, by treating t ♦ t as a con- 
stant, produced small, 10 percent, errors (see figs. 20 and 21) and 
the solution for constant v was employed: 

C(e,4>) ■ ^ (erfc(e - 4 >) ♦ exp(464>)erfc(6 ♦ ♦)} 

where 

,.-!z 

2 Ot 




2/ W(i + 




(t+T)M/2 vs 





Z/ 



( t + T ) 1 /2 vs 
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2/W(l + l ) 


Figure 22 . - Illustration of velocity approximations of the 
S 1 0p/s 1 Interface for the sample 3.5/15. 
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The SIMS profiles for the samples as well as the predicted fits 
are shown In figures 23 to 25. Good fits could be obtained provided 
different values of the diffusion coefficients were allowed for near 
surface and deep profiles. The fact that differing values were re- 
quired may be due to differences In composition as a function of 
position. The difference may also be contributed to by the known 
approximations In the model, which Include: 

(I) The neglectlon of reactions rate constants are both the gas- 
oxide Interface (where there Is qualitative evidence that It Is re- 
quired) and oxide-metal Interface, 

(II) The assumption that the boundaries move at a rate which Is 
known to be different than reality, and 

(III) The assumption that the growth kinetics are pure 
parabolic. 

This knowledge, coupled with the fact that the apparent dlffuslvltles 

differ only by a factor less than five leads to the conclusion that 

these are fairly reliable estimates. The values are only the low 

side of the reported values for oxygen dlffuslvlty In v-S10 2 . This 

may be due to the fact that tracer uptake studies cannot separate out 

the components of the two mechanisms. 

Values for the parabolic rate constant may be extracted from 

the boundary velocity term. The value obtained at 960* C, 

2 

K p « 9700 A /min compares well with a linear Interpolation of the 

2 

data of Irene for a 25 ppm H^O Kp ® 8000 A /min. The values obtained 
at 1000° C k p « 30,900 A 2 /m1n and K p = 28,000 A 2 /m1n for the 18 and 
48 hour samples respectively, also compare well with Irene's values 
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Figure 24 . 
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2 

K p • 12,000 A /rain. The differences could easily arise from low 
levels of water contamination. 

The results are summarized In table A-I. 

V. Summary 

It has been shown that SIMS Is an excellent tool for the study 
of silicon oxidation through the double oxidation technique. This Is 
primarily due to the very good depth resolution. 

It has also been shown that In a double oxidation experiment 
there are two Independent fluxes of tracer Into the scale, a self 
diffusive flux within the network and a permeative flux through the 
scale to the metal, further, the flux of oxygen through the network 
Is pure self diffusion and does not contribute to scale growth. The 
motion of the network Ionic species Is prevented due to the very low 
electrical conductivity. It Is plausible that scales grown on 
heavily doped silicon may contain Impurities sufficient to Increase 
the conductivity and allow network diffusion to contribute to growth. 
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It has been observed that heavily doped silicon displays Increased 
rates of oxidation^ 21 ^ . 

Finally, a model Is given which, though refinements are needed, 
predicts very well the tracer profiles In the double oxidation 
experiments . 

VI. Suggestions for Future Work 

1. The apparatus employed In this work was not well suited for 
either low P u _ experiments or for switching the Isotopic atmosphere 

MpU 

at temperature. All further work should be carried out In a thought- 
fully designed setup. 

2. The temperature range Investigated was very limited. It 
would be an advantage to extend It. It Is also Important to carry 
out experiments at several time ratios at each temperature. 

3. The development of the near surface profile under a battery 
circuit such as Jorgenson's holds the possibility of determining the 
effective charge on the Ionic species diffusing In the network. 

4. It would be of Interest to experiment with heavily doped 
silicon. If there Is a network component to oxidation It should pro- 
duce a detectable convection effect on the near-surface profile. 
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APPCNDIX B 

The Thermal Oxidation of Silicon Carbide 
I . Introduction 

Silicon carbide, SIC, In Increasing use as a high temperature 

structural material. Under oxidizing conditions. Its usefulness In 

this capacity Is often limited by degradation due to oxidation. 

During passive oxidation silicon carbide, as well as other silicon 

based materials (e.g., S1.N. and MoSI,.) are protected through 

J a i 

the formation of a v-SIQ^ film analogus to that formed on silicon 
metal . 

Another potential application of SIC Is as a semiconductor mate 
rial 1 . In this case thermal oxidation may be used to form a protec- 
tive film for the SIC surface during device fabrication. 

A number of studies of the passive oxidation of SIC have been 
carried out*, however there remain points which are not well under 
stood. During passive oxidation the rate of oxidation Is controlled 
by mass transport across the v-SIO^, the thicker the scale becomes 
the slower the oxidation rate. 


*Thomas K. Glasgow (NASA Lewis Research Center) had undertaken a 
thorough, unpublished literature review on this topic, which was 
drawn upon for this discussion. 
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2 

In 1964, Moltzfeldt proposed that the passive oxidation. In 
dry air, of both silicon and silicon carbide was controlled by the 
Inward transport of oxygen from the gas across the scale to the 
metal or carbide. This was based on the similarity of the activation 
energies for the oxidation of both materials and that for the 
permeation of oxygen through a v-SIO,, membrane. As discussed In 
Appendix A, this appears to be correct for the case of silicon, 
however for SIC the situation appears more complex. In addition to 
the production of S10 and SIO^, the oxidation of SIC can produce 
other products, either CO, C0 2 or ternary S1-C-0 compounds. The 
rate of the oxidation reaction Is dependent both on the rate at 
which the oxidants arrive at the Interface and the rate at the 
products are removed. Since at least 1959^, It has been 
recognized that It Is possible for outward diffusion of secondary 
oxides to be responsible for the observation of parabolic kinetics. 
It has been proposed by several workers^' 6 ^ that the controlling 
step Is the outward diffusion of CO. However, Pampush, et al. 7 
examined oxidized powders using Infared spectroscopy and found no 
absorption bands which could be attributed to the double bonds found 
In CO and C0 2 . Instead they found bands corresponding to single 
C-0 bonds and proposed that the reaction product to be a ternary 
S1-0-C phase. Further complexity Is found In the work of Fltzer and 
Ebl, this work reports three different activation energies for the 
passive oxidation of SIC. They propose that above 1400° C the 
outward diffusion of CO Is controlling; between 1200° to 1400° C 
short circuit transport along crlstoballte Interfaces con- 
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trols; and below 1200° C the Inward transport of oxygen through bulk 
v-SIO^ controls . 

In addition, as a result of the polar form of the crystal, the 
oxidation kinetics are very sensitive to crystallographic orienta- 
tion. Harris^ reports that the (0001) carbon face exhibits oxlda 

tlon rates seven times greater than that of the (0001) silicon face. 

( 9 ) 

This has been confirmed by Hunch and Pfaf fender as well as 

Costello^ . 

The current work Involved the double oxidation technique 
(described In Appendix A) using single crystal silicon carbide, over 
the temperature range 1300° to 1400° C, In an effort to gain some 
understanding of the mechanlstlcs of SIC oxidation. The results 
Indicate that the mechanism for oxidation of SIC Is different from 
the non-interacting Interstitial mechanism for the oxidation of 
silicon. 


1 1 . Experimental Procedure 

Samples The experiments done In this work were collaborative 
with either Diane Mleskowskl* (CWRU) or John Costello** (PSU), 
details of the samples may be found In references 10 and 11. 
Essentially the samples were high purity, light green, single crys- 
tals. All experiments were carried out exposing the (0001) carbon 
face, with the more rapid kinetics, to the oxidizing atmosphere. 


♦Currently at NASA Lewis Research Center, Cleveland, OH. 
♦‘Currently at Army Materials Research Lab., Ft. Monmouth, NJ. 
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Experlmental Procedure The apparatus described In the main 
text and also used for the silicon oxidation, was used to carry out 
these thermal oxidations by the procedure described In Appendix A. 
Analysis of tracer profiles was by both proton activation and 

SIMS. 

III. Results and Discussion 

The double oxidation experiments were carried out at two tem- 
peratures, 1300° and 1400° C. Several time ratios of Isotopic 
exposure were used, which are outlined In table B1 , note sample 1 
received three separate oxidations. In all cases a uniform concen- 
tration of tracer throughout the scale was found. This result may 
be contrasted with the results of double oxidation results on sili- 
con metal . 

TABLE B1 

Sample Temperature Isotopic exposure 

number 

1st 2nd 3rd 


1* 

1400° C 

20 

hr 

0-16 

20 

hr 0-18 

2* 

1400° C 

20 

hr 

0-18 

60 

hr 0-16 

3** 

1300° C 

2 

hr 

0-16 

3 

min 0-18 

4** 

1300° C 

2 

hr 

0-16 

11 

min 0-18 

5** 

1300° C 

2 

hr 

0-16 

35 

min 0-18 


•with D. Melskowskl. 
••with J. Costello. 
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As discussed In Appendix A the profiles from silicon have four 
main features; 

(1) A major accumulation of tracer near the oxide metal Inter- 
face Indicating the oxidation reaction took place at this Interface. 

(2) A region of very low tracer concentration In the center of 

the scale Indicating the transport of oxidant Is by a non-interacting 
mechanism. , 

(3) A region of tracer enrichment at the gas-oxide Interface 
developed through pure self diffusion. 

(4) A sharp Interface, ~10 A, between the oxide and the metal. 

A profile from a silicon sample oxidized for 3.5 hrs In 0-16 

and then 15 hrs In 0-18, at 960° C, Is shown In figure 1. A profile 
from a silicon carbide sample oxidized for 2 hrs In 0-16 and then 
35 min In 0-18 Is shown In figure 2. The counting rate for 0-16 and 
0-18 are parallel throughout the profile. Also the interface between 
the oxide and the carbide Is wide, ~500 A. The counting rate for 
carbon gradually Increases across the Interface. The counting rate 
for silicon apparently decreases however this Is an artifact due to 
oxygen enhancement, which Is well known In SIMS. This type of 
Interface was observed by Suzuki, et al.^. 

The uniform concentration of tracer throughout the scale Is In 
consistent with control by the transport of a non-interacting defect. 
In the case of control by Inward diffusion of non-interacting oxygen 
a profile similar to those observed on silicon would be expected. 

If, on the other hand, control Is by outward diffusion of non- 
interacting carbon monoxide the tracer profile should reflect gas 
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Figure 1. - SIMS profile of the v-S10? scale formed on silicon 
metal thermally oxidized at 960° C for 3.5 hr under 0-16 and then 
for 15 hr under 0-18. 



COUNTS 
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Flgure 2 . - SIMS profile of the V-SIO 2 scale formed on silicon 
carbide thermally oxidized at 1300° C for 2 hr under 0-lb and 
then for 35 min under 0-18. 
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wlth, and self diffusion Into the silica film. If v1)t Is many 
times the thickness of the scale, self diffusion could produce a 
uniform tracer concentration. However VBt In all experiment was 
less than the scale thickness. In some cases as small as a tenth of 
the scale thickness. In these cases surface exchange and self diffu- 
sion would be expected to produce a profile with a maximum at the 
gas-oxide Interface which decreases to zero toward the metal-oxide 
Interface. 

The profile would tend toward uniformity If the Inwardly 

diffusing oxidant Interacted frequently with the network, though the 
H21 

work of Clark' 7 Indicates at least a small gradient would always 
exist. Also the silicon oxidations, though at 1000° C, Indicate 
that Inward diffusing oxygen Is non-interacting. A second circum- 
stance by which a uniform profile could be produced Is If exchange 
and diffusion through the dissolved, outwardly diffusing, reaction 
product phase was much faster than the exchange and diffusion Into 
the bulk silica. The result of this would be a uniform source of 
tracer across the scale from which the tracer could diffuse. 

IV. Conclusions 

1. Double oxidation experiments Indicate that there are mech- 
anistic differences between the thermal oxidation of SIC at 1300* to 
1400° C and the thermal oxidation of SI at 960° to 1000° C. 

2. The uniform tracer profiles which occur during double oxida- 
tion experiments Indicate that transport of either the oxidant of the 
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reaction product proceeds by a mechanism which Is Interactive with 
the network. 

3. The Interface between the scale and the carbide Is much wider 
than In the case of silicon and Is enriched In carbon. 

V. Suggestions for Future Work 

1. In order to examine the Interface In more detail It would be 
an advantage to employ an oxygen primary beam In order to avoid dis- 
continuous oxygen enhancement. 

2. It would be useful to extend the temperature ranges of Inves- 
tigation for both SIC and SI to allow more direct comparison. 

3. A systematic Investigation of the carbon distribution, over 
a range of temperature. Is needed. 
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APPENDIX C 

Solutions to the Diffusion Equation Including Convection 


Introduction 

Many diffusion problem of Interest contain moving boundaries or 
occur under an external driving force. Both of these situations 
Introduce a convection term Into the diffusion equations. Pick's 
First law may be written 

J ■ -D ~ ♦ vc (C-l) 


where 

J Is the flux of the species of Interest 
D Is the diffusion coefficient of the species 

~ Is the concentration gradient of the species 

v Is the convective velocity 

c Is the concentration of the species 

When v Is not a function of position. Flck's Second law 

3c 3J n a 2 c „ 3c 

it * ' ii ’ D ^5 • v ii 


(C-2) 


The time dependence of v, however, varies depending on the 
problem of Interest. For a double oxidation experiment (See 
Appendix A) the velocity may be written In the form 


v ■ 


(C-3) 
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where 8 1 and t are both constants. In the case of 
dissolution, or growth, of crystals limited by the rate of diffusion 
through the metl (see Chapter^) the velocity Is of the form 


v 



( C-4) 


where 6^ Is a constant. When diffusion of a charged species 
takes place under a constant applied fteld a velocity Is Introduced 
which Is constant 

v = B 3 ( C-5) 


Note that (C-4) and (C-5) are the limiting forms of (C-3) 


v = 




for t » 


and 


(t ♦ t) 


B 1 B 1 

r-rr = — -* const. For t « t 

I /<? T 


It Is also the case that In many diffusion profiles are sufficiently 
small with respect to the experimental sample that the problem may 
be considered one dimensional and semi-infinite. Boundary 
conditions for this type of problem are 

c(x,t) = 1 For x * 0 t > 0 ( C-6 ) 

c(x,t) = 0 For x -* ® t > 0 

and the Initial condition Is 

c(x,t) = 0 For 0 < x < ® t = 0 
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In the following the solution to the partial differential equation 

(C-2) under the boundary conditions (C-6) will be obtained for two 

- 1/2 

cases. The first will consider the velocity to depend on t 
and the second will consider a constant velocity. 

CASE I 


t 1/2 ax 


I 


PDE 

0 


at 

BC s 

c(o,t) 


c(o,t) 

IC 

c(x,o) 


It proves convenient to first transform the Independent variable x. 
Define z * x - 2 ^2 ' t « x - 2B 2 t 1/2 (C-l 

and t* t 

transforming the terms of the PDE gives 

.. 3c . az ac 3r ac. 

' ax * az ax + 3 t ax “az 

2 2 

a c a .ac.az a c 

' 2 0 az'az'ax “ 2 

a x az 


«i n C. ,ac. az ac az 
<at>x * <at>a « * at at 


at ‘ t l /2 az 
combining this term finds 


ac B 2 ac n a 2 c 
* ~ t 1/2 32 " az 2 


2 ac 

*1/2 az 
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canceling terms gives 


i£ n a_£ 

at = 0 2 

az 


A known solution to this' 1 Is 

c(z,t) = A erfc ( 

where A Is a yet to be determined constant. Transforming the 
solution back to x-space gives 

x fi 2 

c ( x , t ) = A erfc ( — -r£r) ( 

2 Vot vu 

By application of the boundary condition 
C(o.t) = 1 

It can be seen that the constant A Is given by 

B -1 

A = [erfc (-^)] I 

and the final solution Is 


c(x,t) = 


erfc '2-s/Pt 


" \/0 ) 


erfc (- — ) 

Vo 


This may be written In terms of the dimensionless parameters 


and 9 * 

2 Ot VO 


c{©,9) = 


erfc (6 - 
erfc (-«p) 


(C-14) 


Plots of this function versus Theta for various values Phi are shown 
In figures 1 and 2. A transition from a pure erfc, at q>=0, to a 


Cone 
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C<e,$)-erfc<e-$)'erfc<-$> 



Theta -> 


Figure 1. - Plots of the solution to the diffusion equation for 
v*B/t 1/2 for various 
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CCG, $)«erf c <©-$>) /erf c(-$) 



M v U f l0tS ° f the solut1on t0 the dlffusl 
v=B/t 1/2 for various 


equation for 
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normallzed and displaced erfc, <p>3 can be seen. This behavior 
arises because the value of the denominator In C-14 becomes 
essentially constant for <p>3, see figure 3. 

CASE II 


BCS c(o,t) = 1 
c(-.t) = 

IC c(x,o) « 0 

In this case It Is convenient to transform the dependent variable 
c.< 3 > 


Define c(x,t) » u(x,t) • w(x,t) 
where u(x,t) = exp (B 3 (x-B 3 t/2) /2D) 


In order to slmlpllfy the transformation of the equation the chain 


rule Is applied to give 

ac 8c au 3c a« . ac ac au ac aw 

at “ au at aw at ax 3 au ax aw ax 


Since 


ac 

au 


w 


and 



u 


These become 


ac 

at 



♦ u 


aw 

at 


and 


ac 

ax 



♦ u 


aw 

ax 


Now transforming the Individual terms of the PDE gives 


D 


ac 

ax 



♦ u 


aw 

ax 


Since 


All A 

- fj exp(B 3 (x-B 3 t/2)/20) = ^ exp(B 3 (x-B 3 t/2)/2D) 


C($)-erfc(-$) 
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Figure 3. - Plot of elfc (-$) versus 
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ff . exp<B 3 <x-8 3 t/2)/2D> ‘ tg‘^1 

2 

”» ■ fc <8> ■ 8 (exp(e 3 {*-B 3 t/2)/?0) • tf« ♦ 


A^w Aw Aw 

i-f - .,p( 63 (x- 6 3 t/ 2 )/ 2 0 ) ' ^ 4 _| 


"'>8 

Since 


du aw 

w • it * u • it 


_g c 

ff = ft (exp(B 3 (x-B 3 t/2)/20)) = ^ exp(B 3 ( x-B 3 t/2)/20) 
8 . e*p(6 3 (x-B 3 t/2)/20) tf? - « - f|] 


Since It occurs In all of the terms the exponential may be divided 
prior to the combination which yields 


aw D f a!w ?3 ay 1^3, 6 f aw 1^3, 

= 0 { ax 2 *0 ax 4 40 2 ] - 0 3 { 3x 4 20 } 


at 40 

Hultlpllng through and canceling terms gives 


( C— 1 5) 


aw n a 2 w 


transforming the boundary conditions with 
w( x ,t) = c(x,t)/exp(B 3 (x - (B 3 t/2)/2D) gives 


c(o,t) - 1 w(o,t) « 


exp( -B 3 t/4D) 

c(®,t) * 0 w(®,t) * o/® = 0 


exp(B 3 /40) 
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the Initial condition becomes 

c{ x ,o) = 0 w(x,o) = o/exp(Bx/2D) = 0 
This new problem may be solved using Laplace transforms. 
Transforming the equation' ' gives 




where w Is L[w] and p Is the Laplace variable, 
boundary condition gives 


Transforming the 


3 - 1 

w(o,t) = exp jjp -> w(o) = ^ 

P~ B 3 

4D 


w(«,t) = 0 -♦ w(oo) = 0 

By Inspection the solution satisfying these conditions Is 
w(x) = — ^ expC-l/jjx) 

p-b£ u 

<TT> 

( 2 ) 

The Inverse transformation of this solution, from tables , Is 


B ? t 

w(x.t) = l [exp (-—) exp (—) erfc {■ X 


x ®3 x 

♦ exp (^p) exp { — — 
2 Vot 


-xB, 
' 20~ 


2 Ot 


vf 


♦Vsf-u 


(C-16) 


finally transforming back to the dependent variable c(x,t) gives 

xB 


1 


c ( x , t ) = j [ erfc ( 


2 VfTt 


_J¥. 

▼ 4D 


exp (pp) erfc ( — - — 
0 2 Vot 


T 40 


(C-17) 
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Thls may be written In the dimensionless variables 


e 


X 

2 Ot 


and 9 = 


eit . b t 

40 2 0 


as 


c(©,v) ■ 2 C®rfc(© 9) ♦ exp(40 <p)erfc(6 ♦ 9)] (C-18) 

Plots of this function versus theta for various values of phi are 
shown In figures 4 and 5. The behavior of this function Is similar 
to that derived previously. It undergoes a transition from pure 
erfc, when phi * 0 to very nearly a normalized and displaced erfe as 
phi gets big. Figure C-6 show a plot of the Individual terms of the 
solution and their sum. It can seem that the term erfc (e - 9) 
begins to dominate rapidly with Increasing 9 while the term 
exp (4e 9) erfc (6 ♦ 9) gives a gausslan shaped wave which 
travels on the leading edge of the convective wave. 

Conclusions 

Solutions to the diffusion equation have been obtained for a 

semi-infinite slab when constant convection or convection with a 

1 /2 

time dependency of t exists. 


cc$+e)3* j »<<freto<* x «+<$-e) 3 i j «:i#s'0-c$'0)D 
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Figure 4. - Plots 
with constant co 




C<(>+0) o i jO (^©^) dxo+ <^-©> o ^ J ®^S*0“C^'e)D 
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Figure 5. - Plots of the solution to the diffusion equation for 
constant convection, for various ♦. 





C(0,$)=0.5#Cerfc( 0-<f) ) +exp ( 40$ )e rf c (0+(J> ) 3 
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Figure 6. - Plot of the components, and sum, of the solution to 
the diffusion equation for constant convections. 
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